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Within the past 20 years or so, there has occurred an explosion of interest in the magnetic 
behavior of pyrochlore oxides of the type ^2^52^07 where A is a rare-earth ion and B is usually 
a transition metal. Both the A and B sites form a network of corner-sharing tetrahedra which is 
the quintessential framework for a geometrically frustrated magnet. In these systems the natural 
tendency to form long range ordered ground states in accord with the Third Law is frustrated, 
resulting in some novel short range ordered alternatives such as spin glasses, spin ices and spin 
liquids and much new physics. This article attempts to review the myriad of properties found 
in pyrochlore oxides, mainly from a materials perspective, but with an appropriate theoretical 
context. 
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Part I 

Introduction 

I. OVERVIEW 

Competing interactions, or frustration, are common in 
Systems of interacting degrees of freedom. The frustra- 
tion arises from the fact that each of the interactions in 
competition tends to favor its own characteristic spatial 
correlations. A more operational definition classifies a 
systém as frustrated whenever it cannot minimize its to- 
tal classical energy by minimizing the interaction energy 
between each pair of interacting degrees of freedom, pair 
by pair. Frustration is ubiquitous in condensed matter 
physics. It arises in hquid crystals, in superconducting 
Josephson junction arrays in a magnetic field, in molecu- 
lar crystals such as sohd N2, and in magnetic thin films. 
Frustration is important outside the realm of condensed 
matter physics in, for example, the competition between 
nuclear forces and long-range electrostatic Coulomb in- 
teractions between protons which is beheved to lead to 
the formation of a so-cahed "nuclear pasta" state of spa- 
tially modulated nuclear density in stellar interiors. This 
enhances the scattering cross-section for neutrinos with 
nuclear matter and may be a factor in the mechanism of 
super novae explosions. Yet, perhaps, the most popular 
context for frustration is in magnetic systems. This re- 
view is concerned with the many interesting and often ex- 
otic magnetic and thermodynamic phenomena that have 
been observed over the past twenty years in a broad fam- 
ily of geometrically frustrated magnetic mater ials, the 
pyrochlore oxides. Before launching into the review per 
se, we first briefly outline what is meant by geometrie 
magnetic frustration and comment on the scientific mo- 
tivation for investigating geometrically frustrated mag- 
netic systems. 

Although the details will be provided later, here we 
give the reader a brief introduction to the key issues in 
frustrated magnetic systems. One can distinguish two 
classes of frustrated magnetic systems: those in which 
the frustration is geometrie and those in which it is ran- 
dom. The simplest example of geometrie frustration is 
that of Ising spins which can only point in two possible 
(up or down) directions, interacting via nearest-neighbor 
antiferromagnetic exchange and which lie at the vertices 
of an equilateral triangle as shown in Fig. 1. Under those 
conditions, it is impossible for the three spins to align mu- 
tually antiparallel to each other. When many triangles 



are condensed to form an edge-sharing triangular lattice, 
a massive level of configurational spin disorder results, 
giving rise to an extensive residual zero temperature en- 
tropy and no phase transition down to absolute zero tem- 
perature (Houtappel, 1950; Wannier, 1950). This exam- 
ple is one of geometrie frustration since it is the regular 
periodic structure of the (here triangular) space lattice 
that imposes the constraint inhibiting the development 
of a long range ordered state given antiferromagnetic in- 
teractions between the spins. The same systém with Ising 
spins on a triangular lattice, but with ferromagnetic in- 
teractions, is not frustrated and displays a second order 
phase transition to long range order at nonzero temper- 
ature. 

One can distinguish two types of random frustration. 
The first one arises due to competing interactions among 
degrees of freedom in full thermodynamic equilibrium, 
as in the case of liquid crystals and the nuclear pasta 
mentioned above. Although it is not commonly referred 
to as such, one could speak of "dynamical frustration" 
or "annealed frustration" , by analogy to annealed disor- 
der. Since there is rarely only one type of interaction 
setting a sole characteristic length scale in any realistic 
systém, annealed frustration is really the norm. In this 
context, the most dramatic signatúre of dynamical frus- 
tration occurs when multiple length scales develop, as in 
modulated phases of condensed matter, for example the 
case of stripes in the copper oxide superconducting ma- 
terials, or in the case of reentrant phase transitions, as in 
some liquid crystal mater ials. In other words, a systém 
with competing interactions will often attempt to resolve 
the underlying frustration by developing non-trivial spa- 
tial correlations, such as modulated structures, whenever 
full dynamical equilibrium is maintained. We will return 
below to the topic of multiple dynamical degrees of free- 
dom and annealed frustration when we discuss the role 
that lattice degrees of freedom may play in geometrically 
frustrated magnetic systems. 

The other case of random frustration arises when the 
randomness is "quenched", i. e. frozen in. This occurs 
when a subset of the degrees of freedom evolve on a 
time scale that is, for all practical purposes, infinitely 
long compared to the predominant degrees of freedom 
under focus. Consider the example above of the tri- 
angular lattice Ising antiferromagnet. There, the posi- 
tion of the atoms is quenched and only the magnetic 
(spin) degrees of freedom are dynamical. This is an 
example of quenched frustration. However, the frus- 
tration in that case is said to be geometrie^ not ran- 
dom. Quenched random frustration arises in systems 
where the frozen degrees of freedom (e. g. positions of 
the magnetic atoms) are not related by periodic trans- 
lational invariance. One example of random frustra- 
tion is that of magnetic iron (Fe) atoms in a gold (Au) 
matrix. In AuFe, the interaction between the Fe mo- 
ments is mediated by conduction electrons via the RKKY 
(Rudermann-Kittel-Kasuya-Yoshida) Jrkky(^) interac- 
tion, which can be either ferromagnetic or antiferromag- 
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netic depending on the dištance r between two magnetic 
Fe atoms. Depending on their relative separát ion diš- 
tance, some trios of Fe spins will have either one or three 
of their Jrkky interactions antiferromagnetic, and will 
therefore be frustrated. Because the atomic positions are 
frozen-in, this is an example of quenched random frustra- 
tion. Quenched random frustration is a crucial ingredi- 
ent in the physics of spin glass materials (Binder and 
Young, 1986). These systems exhibit a transition be- 
tween a paramagnetic state of thermally fluctuating spins 
to a glass-like state of spins frozen in time, but random 
in direction. Random frustration can also arise when the 
magnetic moments in an otherwise disorder-free geomet- 
rically frustrated systém are diluted, or substituted, with 
non-magnetic atoms (Binder and Young, 1986; Villain, 
1979). 

While it is common to invoke exchange interactions 
between spins as the source of their magnetic coupling, 
there are cases, discussed in some detail below, where 
magnetostatic dipole-dipole interactions play a crucial 
role. Because of their long-range náture and angular 
dependence, dipolar interactions are intrinsically frus- 
trated, irrespective of the lattice dimensionality or topol- 
ogy. Indeed, two dipóle moments /l^i and /j.2 at posi- 
tions r2 and ri {r 12 = r2 — ri) reach their minimum 
energy for \\ /J^2 \\ ^12- Obviously, this condition can- 
not be met for any systém of more than two dipóle mo- 
ments that are not all located on a perfectly straight line. 
As a result, even systems where long-range collinear fer- 
romagnetism is due to dipolar interactions, are in fact 
frustrated, and display quantum mechanical zero point 
íiuctuations (Corruccini and White, 1993; White et al.^ 
1993). This example illustrates that even systems which 
achieve a seemingly simple, conventional long range or der 
can be subject to frustration of their underlying micro- 
scopic interactions. Being intrinsically frustrated, the in- 
troduction of randomness via the substitution of dipole- 
carrying atoms by non-magnetic (diamagnetic) species, 
can lead for sufRciently large dilution to a spin glass 
phase (Binder and Young, 1986; Villain, 1979). 



II. MOTIVATION FOR THE STUDY OF FRUSTRATION 

Magnetism and magnetic materials are pervasive in ev- 
eryday life, from electric motors to hard disk data stor- 
age. From a fundamental perspective, magnetic mate- 
rials and theoretical models of magnetic systems have, 
since the originál works of Ising and Potts, offered physi- 
cists perhaps the best test bench to investigate the broad 
fundamental concepts, even at times universal ones, un- 
derlying collective phenomena in náture. The reasons 
for this are threefold. Firstly, from an experimental per- 
spective, magnetic materials present themselves in var- 
ious aspects. They can be metallic, insulating or semi- 
conducting. As well, the magnetic species may reside on 
crystalline lattices which are spatially anisotropic, pro- 
viding examples of quasi-one or quasi-two dimensional 



systems and permitting an exploration of the role that 
spatial dimensionality plays in phase transitions. Sec- 
ondly, and perhaps most importantly, magnetic mate- 
rials can be investigated via a multitude of experimen- 
tal techniques that can próbe many aspects of magnetic 
and thermodynamic phenomena by exploring spatial and 
temporal correlations over a range of several decades of 
length or time scales. Finally, from a theoretical view 
point, magnetic materials can often be described by well- 
defined microscopic Hamiltonian models which, notwith- 
standing the mathematical complexity commonly asso- 
ciated in solving them, allow in principle to develop a 
theoretical framework that can be used to interpret ex- 
perimental phenomena. This close relationship between 
theory and experiment has been a key characteristic of 
the systematic investigation of magnetism since its incip- 
iency. Such a symbiotic relationship between experiment 
and theory has been particularly strong in the context 
of investigations of frustrated magnetism systems. In- 
deed, a number of experimental and theoretical studies in 
frustrated magnetic systems were originally prompted by 
theoretical proposals. We briefly mention two. The first 
pertains to the proposal that frustrated antiferromagnets 
with a non-collinear ordered state may display a tran- 
sition to a paramagnetic state belonging to a different 
("chiral") universality class from conventional O (N) uni- 
versality for collinear magnets (Kawamura, 1988). The 
second, perhaps having stimulated the largest effort, is 
that some frustrated quantum spin systems may lack con- 
ventional semi-classical long-range order altogether and 
possess instead a quantum disordered "spin liquid" state 
breaking no global spin or lattice symmetries (Anderson, 
1973, 1987). 

Systems where magnetic moments reside on lattices of 
corner-sharing triangles or tetrahedra are subject to geo- 
metrie magnetic frustration, as illustrated in Fig. 1, and 
are expected to be ideál candidates for exhibiting large 
quantum mechanical spin fluctuations and the emergence 
of novel, exotic, magnetic ground states. 

The cubic pyrochlore oxides, ^2^207, have attracted 
significant attention over the past twenty years because 
the A and B ions reside on two distinct interpenetrating 
lattices of corner-sharing tetrahedra as shown in Fig. 2. 
Henceforth, we shall refer to the lattice of corner-sharing 
tetrahedra simply as the pyrochlore lattice, as has be- 
come customary since the mid 1980's. If either A, B or 
both are magnetic, and the nearest-neighbor exchange 
interaction is antiferromagnetic, the systém is highly ge- 
ometrically frustrated. As a result, antiferromagnetically 
coupled classical Heisenberg spins on the pyrochlore lat- 
tice do not develop long range order at any nonzero tem- 
perature, opening up new avenues for novel, intrinsically 
quantum mechanical, effects to emerge at low tempera- 
tures. 

In all real systems, interactions beyond isotropic 
nearest-neighbor exchange are at play. Exchange be- 
yond nearest-neighbors, symmetric and antisymmetric 
(Dzyaloshinskii-Moriya) exchange, dipolar interactions 
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FIG. 1 Antiferromagnetically coupled spins arranged on a 
triangle or tetrahedron are geometrically frustrated. 



and magnetoelastic couplings are all concrete examples 
of corrections to the isotropic nearest-neighbor exchange 
part of the Hamiltonian, Tíu- These additional inter- 
actions, are typically weaker than Hu and, often, 
two or more terms in H' compete to dictate the de- 
velopment of the spin-spin correlations as the systém 
is cooled from the paramagnetic phase. The specific 
details pertaining to the materiál under consideration 
determine whether the perturbations H' are able, ulti- 
mately, to drive the systém into a semi-classical long- 
range ordered state, or whether quantum fluctuations 
end up prevailing. It is for this reason that many mag- 
netic pyrochlore oxides have been found to display a 
gamut of interesting and unconventional magnetic and 
thermodynamic behaviors. Examples of phenomena ex- 
hibited by A2B2O7 materials include spin glass freezing 
in Y2M02O7, spin liquid behavior in Tb2TÍ207, disor- 
dered spin ice (not to be confused with spin glass) phe- 
nomenology in H02TÍ2O7 and Dy2TÍ2 07, ordered spin 
ice in Tb2Sn207, order-by-disorder in Er2TÍ2 07, uncon- 
ventional anomalous Hall effect in metallic Nd2Mo207, 
superconductivity in Cd2Re207 and the Kondo-like ef- 
fect in Pr2lr207, to name a few. 

The above examples illustrate the broad diversity of 
phenomena exhibited by pyrochore oxides. There exist 
a number of generál reviews concerned with pyrochlore 
structure materials such as those by Subramanian et al 
(1983), Greedan (1992a), and Subramanian and Sleight 
(1993) and more focused reviews for example on the rare 
earth titanates, (Greedan, 2006). As well, pyrochlore 
structure materials are often featured in more generál 
reviews of geometrically frustrated magnetic materials, 
such as those by Ramirez (1994), Schiffer and Ramirez 
(1996), Greedan (2001) and others (Diep, 1994, 2004; 
Gauhn, 1994). The purpose of this article is to pro- 
vide, given the space limit ations, a comprehensive re- 
view of both experimental and theoretical investigations 
of the A2B2O7 materials over the past twenty years or 
so. Whenever possible, we outline what we perceive as 
promising avenues for future work. Thus, space limita- 
tions preclude detailed discussion of other very interest- 
ing highly frustrated magnetic materials based on lattices 



of corner-shared triangles such as kagome systems and 
garnets, or the spinels and Laves phases, which are also 
based on the pyrochlore lattice. However, we will refer 
to those materials in specific contexts whenever appropri- 
ate in order to either contrast or relate to the behavior 
observed in the pyrochlore oxides. 

This review is organized as follows: after a brief his- 
torical overview and some generalities on geometrically 
frustrated magnets we review and discuss experimental 
data from the most státie (i. e. long range ordered phases) 
to the most dynamic (i. e. spin liquid phases). Then fol- 
lows a brief section on the superconducting compounds 
found among pyrochlore oxides, although frustration has 
not been demonstrated to play a role here. Finally, we re- 
port on salient results from experiments performed under 
the influence of applied magnetic field or pressure. We 
conclude the review by a brief discussion of perceived 
avenues for future research. 



III. HISTORICAL PERSPECTIVE 

The term frustration, meaning the inability of a pla- 
quette of spins with an odd number of antiferromagnetic 
bonds to adopt a minimum energy collinear spin con- 
figuration, first appeared in a paper on spin glasses by 
Toulouse (1977) and soon after, in one by Villain (1977). 
However, it is said that the concept of frustration as a 
key ingredient in the physics of spin glasses was first in- 
troduced by Anderson who apparently wrote on a black- 
board at a summer school in Aspen in 1976 ''Frustration 
is the name of the game" (Chandra and Coleman, 1995; 
Harrison, 2004). It is interesting that the notion of frus- 
tration, random frustration to be more precise, had to 
await the discovery of spin glasses before it was intro- 
duced in magnetism. 

Indeed, geometrically frustrated antiferromagnets, 
without being labeled as such, had been studied much 
earlier. In 1950, Wannier (1950) and Houtappel (1950) 
investigated the exact solution of the Ising antiferromag- 
net on the triangular lattice, finding that no transition 
to long range order occurs at positive temperature and 
that the systém possesses an extensive ground state en- 
tropy. Soon after, following on a previous work by Néel 
(1948), Yafet and Kittel (1952) investigated the problém 
of antiferromagnetic spin arrangements in spinel ferrite 
materials and the role of magnetic interactions between 
the octahedral and the tetrahedral sites in determining 
the ground state spin structure, where the octahedral 
sites form a pyrochlore lattice. Motivated by the prob- 
lém of magnetic ordering for spins on the octahedral sites 
of normál spinels and the related one of ionic ordering in 
inverse spinels, Anderson (1956) investigated the physics 
of antiferromagnetically coupled Ising spins on what we 
now recognize as the pyrochlore lattice. This work ap- 
pears to have been the first to identify clearly the pecu- 
liar magnetic properties of this systém and, in particular, 
noted the existence of a macroscopic number of ground 
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States, and the connection with the problém of the zero 
point protón entropy in hexagonal water ice previously 
investigated by Pauling (1935). We will return to the 
connection between geometrically frustrated Ising spins 
on the pyrochlore lattice and the water ice problém later. 

Quantum mechanics made its first noticeable entry 
in the realm of geometrically frustrated magnetism in 
1973 when Anderson (1973) proposed that a regular two- 
dimensional triangular lattice of spins S = 1/2 coupled 
via nearest-neighbor antiferromagnetic exchange could, 
instead of displaying a three sublattice 120 degree Néel 
ordered state, possess an exotic resonating valence bond 
liquid of spin singlets. This proposal was further exam- 
ined and its conclusion supported by more detailed calcu- 
lations by Fazekas and Anderson (1974). While Anderson 
(1973) had written that there were no experimental re- 
alizations of an 6* = 1/2 antiferromagnet with isotropic 
exchange, the situation would change dramatically al- 
most twenty-five years later with the 1987 discovery of 
the copper-oxide based high-temperature superconduc- 
tors, some of which possess a parent insulating state with 
antiferromagnetically coupled 6* = 1/2 spins (Anderson, 
1987). 

The discovery of spin glasses by Cannella and My- 
dosh (1972) prompted a flurry of experimental and the- 
oretical activity aimed at understanding the náture of 
the observed spin freezing and determining whether a 
true thermodynamic spin glass phase transition occurs 
in these systems. Particularly relevant to this review, 
Villain (1979) investigated the generál problém of insu- 
lating spin glasses and the role of weak dilution of the B- 
site pyrochlore lattice in spinels. The extensive ground 
state degeneracy characterizing this systém for classical 
Heisenberg spins interacting via nearest-neighbor anti- 
ferromagnetic exchange J was noted. He speculated that 
such a systém would remain in a paramagnetic-like state 
down to zero temperature, giving rise to a state at tem- 
peratures T <C J that he called a collective paramagnet 
— a classical spin liquid of sorts. Villain (1979) also con- 
sidered the case of tetragonal spinels where there are two 
exchange couplings, J and J', perpendicular and parallel 
to the tetragonal axis. More recently, theoretical (Tch- 
ernyshyov et al.^ 2002; Yamashita and Ueda, 2000) and 
experimental (Lee et al.^ 2000) studies have investigated 
the cubic to tetragonal distortion in spinels, which leads 
to a reliéf of the magnetic frustration and allows the sys- 
tém to develop long range magnetic or der. Villain dis- 
cussed how the effect of a weak dilution of the magnetic 
moments may be dramatically different in the case of cu- 
bic spinels compared to those with a slight tetragonal 
distortion. In particular, he argued that in the case of 
the cubic spinel, weak dilution would leave the systém a 
collective paramagnet, while for the tetragonal spinel the 
ground state would be glassy. 

While there had been some experimental studies on 
non-collinear magnets between the work of Villain and 
the mid-to-late 1980's (Coey, 1987), the high level of re- 
search activity in the field of frustrated magnetism re- 



ally only started following the discovery of the unconven- 
tional magnetic properties of SrOrg-a^Gas+a^Oig (SCGO) 
by Obradors et al. (1988). In the context of oxide py- 
rochlores, the topic of this review, experimental and theo- 
retical interests developed in the mid-1980's with reports 
of spin glass behavior in the apparently well-ordered ma- 
teriál Y2M02O7 and accelerated rapidly after 1990. 



IV. THEORETICAL BACKGROUND 

The amount of theoretical work dedicated to the topic 
of highly frustrated systems, particularly quantum spin 
systems, has exploded over the past ten years. It is be- 
yond our scope here to discuss more than a fragment and 
in fact the field is so active that a dedicated review would 
at this time be warranted. Short of such being available, 
we refer the reader to a somewhat recent book edited by 
Diep (2004) and a forthcoming one by edited by Mila 
(2008). In this section, we merely present the simplest 
and most generic theoretical arguments of relevance to 
the pyrochlore systems. The focus of the discussion is on 
the effects of perturbations on a nearest-neighbor Heisen- 
berg pyrochlore antiferromagnet. 

One simple and highly convenient measure of the level 
of frustration in magnetic systém is the so-called frustra- 
tion index f (Ramirez, 1994), defined as 



Here ^cw is the Curie or Curie-Weiss temperature, de- 
fined from the high temperature paramagnetic response 
of the systém, i. e. the high temperature, linear part of 
the inverse DC susceptibility. The temperature T* is 
the critical temperature Tc (or Néel temperature, Tn) at 
which the systém ultimately develops long range spin or- 
der. In the case of freezing into a glassy state, which may 
correspond to a genuine thermodynamic spin glass transi- 
tion, T* would be the freezing temperature Tf. The more 
frustrated a systém is, the lower T* is compared to ^cw, 
hence an ideál spin liquid systém would have / = 00. 
In most real materials there are, however, perturbative 
interactions. Tí' that ultimately intervene at low temper- 
ature and lead to the development of order out of the 
spin liquid state. The rest of this section discusses some 
of the most common perturbations. 

Anderson (1956) and Villain (1979) were the first to 
anticipate the absence of long range order at nonzero 
temperature in the Ising (Anderson, 1956) and Heisen- 
berg (Villain, 1979) pyrochlore antiferromagnet. The 
Hamiltonian for the Heisenberg antiferromagnet model 
is 

nn = -jY.s,-Sj , (2) 

where we take for convention that J < O is antiferromag- 
netic while J > O is ferromagnetic. The spins Si can be. 
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at least within the model, taken either as classical and of 
fixed length I^S^I = 1 or quantum. A model with n = 1, 
n = 2 or n = 3 components of S corresponds to the 
Ising, XY and Heisenberg model, respectively. Since the 
pyrochlore lattice has cubic symmetry, the Ising and XY 
models with a uniform global Ising easy axis direction 
z or XY global piane are not physical. One can how- 
ever, define local Zi Ising directions or local XY planes 
with normál Zi that are parallel to the four cubic (111) 
directions. 

A mean field theory calculation for the Heisenberg 
Hamiltonian by Reimers et al. (1991a) provided indirect 
confirmation of Anderson and Villain's result by find- 
ing that a nearest-neighbor model with classical Heisen- 
berg spins on the pyrochlore lattice has two (out of 
four) branches of exactly degenerate soft (critical) modes 
throughout the Brillouin zóne. In other words, there 
is an extensive number of soft modes in this systém 
and no unique state develops long range order at the 
mean field critical temperature. In early Monte Carlo 
simulations, Reimers (1992) found no sign of a transi- 
tion in this systém, either to long range order or to a 
more complex type, such as nematic order, down to a 
temperature of J/20, where J is the nearest-neighbor 
antiferromagnetic exchange constant. This conclusion 
was confirmed in numerical studies by Moessner and 
Chalker (1998, a). This behavior of the nearest-neighbor 
pyrochlore antiferromagnet differs markedly from that 
of the nearest-neighbor classical Heisenberg antiferro- 
magnet on the two-dimensional kagome lattice where 
co-planar order-by-disorder develops as the temperature 
drops below ~ J/ 100 (Chalker et al.^ 1992; Reimers 
and Berlínsky, 1993). However, the Mermin-Wagner- 
Hohenberg theorem forbids a phase transition at nonzero 
temperature in the kagome Heisenberg antiferromagnet. 
While no transition is seen in the pyrochlore Heisenberg 
systém (Moessner and Chalker, 1998, a; Reimers, 1992), 
the pyrochlore lattice with easy-plane (XY) spins, where 
the easy planes are perpendicular to the local (111) di- 
rections, displays a phase transition to long range order 
(Bramwell et al.^ 1994; Champion et al, 2003). Finally, in 
agreement with Anderson's prediction (Anderson, 1956), 
Monte Carlo simulations find that the Ising antiferromag- 
net, or equivalently the Ising ferromagnet with local (111) 
directions, which is the relevant model to describe spin 
ice materials, does not order down to the lowest temper- 
ature (Harris et al.^ 1998). 

Moessner and Chalker (1998, a) investigated the condi- 
tions under which an n-component classical spin Si with 
\Si\ = 1 arranged on a lattice of corner-sharing units of 
q sites and interacting with q — 1 neighbors in each unit 
can exhibit thermally-induced order-by-disorder. Their 
analysis allows one to rationalize why the Heisenberg 
n = 3 kagome antiferromagnet develops co-planar order 
and why the XY pyrochlore antiferromagnet (with spins 
in a global XY piane) develops collinear order (Moess- 
ner and Chalker, 1998). From their analysis, one would 
predict a phase transition to long range nematic order 



for g' = 3, n = 3 in three spatial dimensions. Interest- 
ingly, such a situation arises in the S = 1/2 Na4lr308 
hyperkagome (Okamoto et al.^ 2007) and Gd3Ga50i2 
garnet (Dunsiger et a/., 2000; Petrenko et a/., 1998; Schif- 
fer et al, 1994, 1995a; Yavors'kii et al, 2006). Moving 
away from the strictly nearest-neighbor model, one would 
generally expect interactions beyond nearest-neighbors 
to drive a transition to an ordered state, at least for clas- 
sical spins (Reimers, 1992). While this has been stud- 
ied via Monte Carlo simulations (Mailhot and Plumer, 
1993; Reimers et al, 1992), it has not been the subject 
of very many investigations. In that context, we note 
that a paper by Pinettes et al (2002) studies a classi- 
cal Heisenberg pyrochlore antiferromagnet with an addi- 
tional exchange interaction J' that interpolates between 
the pyrochlore lattice ( J' = 0) and the face-centered cu- 
bic lattice {J' = J). It is found that for J' / J as low 
as J' / J ^ 0.01, the systém orders into a collinear state 
whose energy is greater than the classical incommensu- 
rate ground state predicted by mean field theory. A 
Monte Carlo study of another classical model of a py- 
rochlore lattice with antiferromagnetic nearest-neighbor 
exchange and ferromagnetic next nearest-neighbor ex- 
change finds an ordered state with mixed ordered and 
dynamical character (Tsuneishi et al, 2007). Because 
the pyrochlore lattice lacks bond inversion symmetry, 
Dzyaloshinsky-Moriya (DM) spin-spin interactions aris- 
ing from spin-orbit coupling are allowed by symmetry 
and can give rise to various ordered states depending on 
the orientation of the DM vector (Elhajal et al, 2005). 

As we will discuss below, a number of pyrochlore ma- 
terials exhibit glassy behavior where the spins freeze at 
low temperature into a state where their orientation is 
random in space, but frozen in time (Binder and Young, 
1986). It is not yet fully understood whether this glassy 
behavior is due to the intrinsic náture of the systems 
considered, as is the case for the spin ices, or if it orig- 
inates from random disorder (e. g. antisite disorder, ran- 
dom bonds or other) as in conventional random dis- 
order spin glasses (Binder and Young, 1986). Only a 
very few theoretical studies have explored the effect of 
random disorder in antiferromagnetic pyrochlores with 
classical Heisenberg spins and with small and homogo- 
neously random deviations of exchange Jij from the aver- 
age antiferromagnetic exchange J (Bellier-Castella et al, 
2001; Saunders and Chalker, 2007). Interestingly, ran- 
dom bonds on a single tetrahedron lift the degeneracy 
and drive an order-by-(random)-disorder to a collinear 
spin arrangement. However, on the lattice, the sys- 
tém remains frustrated and there is no global order-by- 
disorder (Bellier-Castella et al, 2001). Considering this 
model in more detail, Saunders and Chalker (2007) find 
compelling evidence of a conventional spin glass transi- 
tion. Sagi et al (2005) have investigated how coupling 
with the spins to the lattice may lead to a freezing be- 
havior and trapping of the spins into metastable states 
in an otherwise disorder-free nearest-neighbor classical 
pyrochlore antiferromagnet. Villain (1979) discussed the 
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role of non-magnetic impurities on the pyrochlore lat- 
tice antiferromagnet. He argued that the collective para- 
magnetic behavior of the otherwise disorder-free systém 
should survive up to a finite concentration of impurities. 
A similar conclusion was reached on the basis of heuris- 
tic arguments and simulations for the site-diluted two- 
dimensional kagome Heisenberg antiferromagnet (Shen- 
der et al.^ 1993). It would be valuable if the problems of 
site-disorder and dilute random bonds in the antiferro- 
magnetic pyrochlore lattice were explored in more detail. 
We note in passing that the problém of disorder in quan- 
tum variants of the kagome Heisenberg antiferromagnet 
has attracted some attention (Läuchli et al.^ 2007). 

In insulating A2B2O7, where A is a 4f rare-earth triva- 
lent ion, the inter-rare-earth exchange interactions are 
often small, contributing only 10^ — 10^ K to the Curie- 
Weiss 6cw temperature. This is because the unfilled 
magnetically active 4f orbitals are strongly shielded by 
the 5s, 5p and 5d orbitals and their direct exchange 
over lap or superexchange over lap is small. At the same 
time, trivalent rare earths often possess large magnetic 
moments and hence the magnetostatic dipole-dipole en- 
ergy scale is significant compared to 6cw and needs to 
be incorporated into Eq. (2), 

^int = ~ 2 ^ "^^^ ' ^ 

/MO \ (^LMb)^ {Jj • Jj - ^Jj • UjTij ■ Jj) 

Here we have taken the opportunity to generalize 
Eq. (2) to the case where the total angular momen- 
tum Ji = Li -\- Si is included. is the Landé factor 
and r j — Vi = \rij\ŕij^ where r i is the position of mag- 
netic ion i. In generál, one needs to include the role of 
the single-ion crystal field, TYcf, to the full Hamiltonian, 
^ = Hcí H- Hint • 

A mean field study of TYint, with TYcf = O and with 
only nearest-neighbor antiferromagnetic J^n and dipo- 
lar interactions of approximately 20% the strength of 
J7nn, was done by Raju et al. (1999) to investigate the 
type of order expected in Gd2TÍ207. The assumption 
TLci = O is a reasonable íirst approximation for Gd2TÍ2 07 
since the electronic ground state of Gd^+ is ^87/2 with 
L = 0. Raju et al. (1999) argued that this model has 
an infinite number of soft modes with arbitrary order- 
ing wave vector along [111] and that the transition at 
~ 1 K could possibly proceed via order-by-disorder. 
However, Palmer and Chalker (2000) argued that quartic 
terms in the mean field Ginzburg-Landau theory would 
select a /Cord = 000 ordering wavevector at the criti- 
cal temperature Tc that also corresponds to the zero 
temperature classical ground state. Subsequent mean 
field calculations (Cépas and Shastry, 2004; Enjalran and 
Gingras, 2003) showed that there are in fact two very 
closely spaced transitions and that the transition at the 
highest temperature is at a unique ordering wave vector 



A- Site B- Site 




fíg. 2 The A and B sublattices of the cubic pyrochlore ma- 
terials. Either or both sublattices can be magnetic. 



(^ord — \\\) that had been missed in previous calcula- 
tions (Palmer and Chalker, 2000; Raju et al, 1999). A 
Monte Carlo simulation (Cépas et al., 2005) found a sin- 
gle first order transition. It appears likely that the small 
temperature separation between the two transitions that 
is predicted by mean field theory is not resolved in the 
Monte Carlo simulation. Irrespective of the subtleties 
associated with the transition, it is generally accepted 
that the classical ground state of the nearest-neighbor 
pyrochlore Heisenberg antiferromagnet with weak dipo- 
lar interactions is the so-called Palmer- Chalker ground 
state (Cépas and Shastry, 2004; Del Maestro and Gin- 
gras, 2004, 2007; Enjalran and Gingras, 2003; Palmer 
and Chalker, 2000; Wills et al, 2006). Starting from 
the Palmer- Chalker ground state, a 1/5 expansion cal- 
culation shows that all conventional magnetic spin-wave- 
like excitations of this systém are gapped and that the 
l/S quantum fluctuations are negligible (Del Maestro 
and Gingras, 2004, 2007). For Ho and Dy-based py- 
rochlore oxides, the crystal field Hamiltonian TÍcí pro- 
duces a single-ion ground state doublet and the systém 
maps onto a classical dipolar Ising model. We shall re- 
turn to this model when reviewing the problém of spin 
ice in Section X. 

Another important perturbation of relevance to real 
materials is the coupling between spin and lattice degrees 
of freedom. Spin-lattice coupling can lead to a redistri- 
bution of the frustration among the magnetic bonds such 
that the energy cost for distorting the lattice is more than 
compensated by the energy difference between that of 
the more strongly coupled antiparallel spins on the com- 
pressed bonds compared to the higher energy (frustrated, 
parallel) spins on the dilated bonds (Tchernyshyov et al, 
2002; Yamashita and Ueda, 2000). It is anticipated that 
such coupling can lead to a global lifting of magnetic de- 
generacy and drive a transition to long range magnetic 
order that is accompanied by a cubic to tetragonal lattice 
deformation. Villain (1979) had also discussed a number 
of significant differences in the thermodynamic proper- 
ties and the response to random disorder of a tetragonal 
spinel compared to that of a cubic spinel. The problems 
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of magnetic order, lattice distortion, coupling to an ex- 
ternal magnetic field and magnetization plateaus in oxide 
spinels are currently very topical ones (Bergman et al.^ 
2006; Lee et al, 2000; Pene et al, 2004; Tchernyshyov 
et al.^ 2002). However, in pyrochlore oxides, interest in 
the role of the lattice degrees of freedom is just starting 
to attract attention (Ruff et al.^ 2007; Sagi et al.^ 2005). 

We close with a brief discussion of the quantum S = 
1/2 pyrochlore Heisenberg antiferromagnet. This is ob- 
viously an extremely difíicult problém and there is no 
concensus yet about the náture of the ground state. 
It appears to have been first studied by Harris et al. 
(1991) who found indicators that the systém may form 
a dimerized state at low temperature. By applying per- 
turbation theory to the density matrix operátor, Canals 
and Lacroix (1998, 2000) found evidence for a quan- 
tum spin liquid in the S = 1/2 pyrochlore antiferro- 
magnet. From exact diagonalizations on small clusters, 
they also found that the singlet-triplet gap is filled with 
a large number of singlets, similar to what is obser ved 
in the S = 1/2 kagome lattice antiferromagnet (Lechem- 
inant et al, 1997). Calculations using the contractor- 
renormalization (CORE) method (Berg et al.^ 2003) find 
a ground state that breaks lattice symmetry, analogous 
to that suggested by Harris et al. (1991), and with a 
singlet-triplet gap filled with singlets, similar to Canals 
and Lacroix (1998). Moving away from strictly numeri- 
cal approaches, analytical calculations based on large- 
(Moessner et al, 2006; Tchernyshyov et al, 2006) and 
large-^* expansions (Hizi and Henley, 2007) as well as 
variational calculations away from an exact dimer cover- 
ing ground state solution (Nussinov et al.^ 2007) have all 
been very recently carried out. 



Atom 


Wyckoff 


Point 


Minimal 




Position 


Symmetry 


Coordinates 


A 


16d 


3m (Dsd) 


111 

2' 2' 2 


B 


16c 


3m (Dsd) 


0, 0, 


O 


48/ 


mm {C2v) 


x i i 


O' 


8b 


43m (Td) 


3 3 3 
8' 8' 8 



TABLE I The crystallographic positions for the space group, 
FdSm (No. 227) suitable for the cubic pyrochlore ^2^2 060^ 
with origin at 16c. 



cently sparked considerable interest (Lee et al.^ 2000, 
2002; Matsuda et al, 2007; Ueda et al, 2005) in this 
area of magnetic frustration. Another common family of 
materials that contain a network of corner-sharing tetra- 
hedra are the cubic Laves phases, such as YMn2 (Ballou, 
2001; Shiga et al, 1993). 

As already emphasized, in this review, we will restrict 
ourselves to the cubic pyrochlore oxides with the gen- 
erál formula A2B2O7, where A is a trivalent rare earth 
which includes the lanthanides, Y and sometimes Sc, and 
B is either a transition metal or a p-block metal ion. 
It should be mentioned here that it is possible to form 
742^^2^07 (2+,5+) pyrochlores, however these are not 
common and have not been studied in great detail. Most 
of this manuscript will deal with the (3+,4+) variety of 
pyrochlores, but when appropriate the (2+,5+) will be 
discussed. 



V. CRYSTAL STRUCTURE 



Part II 

Materials 

While this review is focused on the three dimensional 
pyrochlore lattice, which is shown in Fig. 2, there are of 
course other lattice types which provide the conditions 
for geometrie frustration. Among these, in two spatial 
dimensions are the edge-sharing and corner-sharing tri- 
angular lattices - the latter best known as the Kagome 
lattice - and in three spatial dimensions one has garnets, 
langasites and face centered cubic lattices among others. 
There exists a considerable and important literatúre de- 
voted to these materials which lies, unfortunately, outside 
the scope of this work. 

A network of corner-sharing tetrahedra is found in 
several common minerál types in náture. The octahe- 
drally coordinated B-site of the spinel class of minerals 
{AB2O4) forms a network of corner-sharing tetrahedra. 
The A- and the B-site can be filled by many ions, but 
when the B-site is occupied by a trivalent magnetic tran- 
sition metal ion, magnetic frustration can play a role in 
its bulk properties. The trivalent chromates have re- 



A. Space Group and Atomic Positions 

Pyrochlore materials take their name from the minerál 
NaCaNb206F pyrochlore, the structure of which was first 
reported by Gaertner (1930). The name, literally "green 
fire" , alludes to the fact that the minerál shows a green 
color upon ignition. Most synthetic pyrochlores of in- 
terest here are oxides that crystallize in the space group 
FdŠm (No. 227). As pointed out by Subramanian and 
Sleight (1993), confusion can exist due to the fact that 
there exist four possible choices for the origin. Currently, 
Standard practice is to formulate oxide pyrochlores as 
^2^2 060' and to plače the B ion at 16c (the origin of 
the second setting for FdSm in the International Tables 
of Crystallography), A at 16(i, O at 48/ and O' at 8b as 
shown in Table L Note that there is only one adjustable 
positional parameter x for the O atom in 48/. There 
are at least two important implications of these struc- 
tural details. One concerns the topology of the 16c and 
16d sites and the other the coordination geometry of the 
O ligands about the two metal sites. Both the 16c and 
16d sites form a three dimensional array of corner-sharing 
tetrahedra as shown in Fig. 2, thus giving rise to one of 
the canonical geometrically frustrated lattices. 



9 




FIG. 3 The coordination geometry of the A-site by O (shaded 
spheres) and 0^(open spheres) atoms. The O' — A — O' unit 
is oriented normál to the average piane of a puckered six- 
membered ring. This bond is one of the shortest seen between 
a rare earth and oxygen ions. 



B. Local Environment of the A and S-Sites 

The coordination geometry about the two metal sites 
is controlled by the value of x for the O atom in the 4^1 
site. For x = 0.3125 one has a perfect octahedron about 
16c and x = 0.375 gives a perfect cube about 16d. In fact 
x is usually found in the range 0.320 to 0.345 and the two 
geometries are distorted from the ideál polyhedra. For a 
typical pyrochlore the distortion about the B ion or 16c 
site is relatively minor. The 3m(Dsd) point symmetry 
requires that all six B - O bonds must be of equal length. 
The O - B - O angles are distorted only slightly from the 
ideál octahedral values of 90° ranging between 81° and 
100°. The distortion of the A site geometry from an 
ideál cube is on the other hand very large. This site is 
depicted in Fig. 3 and is best described as consisting of 
a puckered six-membered ring of O atoms with two O' 
atoms forming a linear - A - O' stick oriented normál 
to the average piane of the six-membered ring. The A - 
O and A - bond distances are very different. While 
typical A - O values are 2.4 -2.5 Ä, in accord with the 
sum of the ionic radii, the A - bonds are amongst 
the shortest known for any rare earth oxide, ~ 2.2 Ä. 
Thus, the A-site has very pronounced axial symmetry, 
the unique axis of which is along a local (111) direction. 
This in turn has profound implications for the crystal 
field at the A-site which determines much of the physics 
found in the pyrochlore materials. 

Recall that the spin-orbit interaction is large in rare 
earth ions and the total angular momentum, J = L -\- S 
is a good quantum number. For a given ion one can ap- 
ply Hunďs rules to determine the isolated (usually) elec- 
tronic ground state. Electrostatic and covalent bonding 
effects originating from the local crystalline environment, 
the so-called crystal field (CF), lift the 2 J + 1 degener- 
acy of the ground state. A discussion of modern meth- 
ods for calculation of the CF for f-element ions is be- 
yond the scope of this review but is described in several 
monographs for example, one by Húfner (1978). For our 
purposes we will assume that the single ion energy levels 
and wave functions, the eigenvectors and eigenenergies of 



the CF Hamiltonian Hcí have been suitably determined 
either through ab initio calculations or from optical or 
neutrón spectroscopy (see for example Rosenkranz et al. 
(2000) and more recently Mirebeau ct al (2007)). 

Hcí can be expressed either in terms of the so-called 
tensor operators or the "operátor equivalents" due to 
Stevens (1952). The two approaches are contrasted by 
Húfner (1978). While the tensor operators are more con- 
venient for ab-initio calculations, the latter are better 
suited to our purposes here. In this formalism, TYcf is ex- 
pressed in terms of polynomial functions of the Jiz and 
J^_^, with Ji± = Jix ± Jiy, which are components of the 
Ji angular momentum operátor. The most generál ex- 
pression for Hcí is: 

^^f = EE^rorw, (4) 

i l,m 

where, for example, the operátor equivalents are O2 = 
3J/-J(J+1) andO^ = J+^ + J_^ The full CF Hamil- 
tonian for -3m (D^d) point symmetry involves a total of 
six terms for / = 2, 4 and 6 (Greedan, 1992a). In fact, due 
to the strong axial symmetry of the A-site, described pre- 
viously, it can be argued that the single / = 2 term, B2, 
plays a major role in the determination of the magnetic 
anisotropy of the ground state. In the Stevens formalism, 
B2 = A2(r^)aj(l-cr2), where A2 is a point chargé lattice 
sum representing the CF strength, (r^) is the expecta- 
tion value of r^ for the 4f electrons, a2 is an electron 
shielding factor and aj is the Stevens factor (Stevens, 
1952). This factor changes sign in a systematic pattern 
throughout the lanthanide serieš, being positive for A= 
Sm, Er, Tm and Yb and negatíve for all others. So, 
the sign of B2 depends on the product ajA2 and A2 is 
known from measurements of the electric field gradient 
from for example ^^^Gd Môssbauer studies to be posi- 
tive for pyrochlore oxides (Barton and Cashion, 1979). 
Thus, B2 should be positive for A = Sm, Er, Tm and 
Yb and negatíve for all others. From the for m of B 2 
above, it is clear that states of different \Mj) do not mix 
and that the energy spectrum will consist of a ladder of 
States with either \Mj^^.^^) {B^ > 0) or \Mj^^.^^) {B^ < 
0) as the ground state. Note that the former constitutes 
an easy piane and the latter an easy axis with respect 
to the quantization axis which is <111> for pyrochlores. 
A comparison of the known anisotropy for the A2TÍ2O7 
and 742Sn2 07 materials, that is, easy axis for A = Pr, 
Nd, Tb, Dy and Ho and easy piane for A = Ei and Yb 
with the sign of B2 shows a remarkable agreement with 
this very simple argument. [Note that only the Stevens 
formalism works here. The tensor operátor definition of 
B2 is B2 = A2(r^), so this quantity is always positive 
for pyrochlore oxides, independent of the rare earth A. 
We do not maintain that the actual ground state wave 
function can be obtained within such a simple model (al- 
though the agreement for A = Dy and Ho is remarkable), 
rather, that the overall anisotropy can be predicted with- 
out a detailed calculation.] 
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FIG. 4 Alternating kagome and triangular planar layers 
stacked along the [111] direction of the pyrochlore lattice. 
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FIG. 5 A structure-field or stability- field map for A2B2O7 
C. Alternatíve Views of the Pyrochlore Structure materials. The B = Pi serieš are also high pressure phases. 

Adapted from Subramanian and Sleight (1993). 



An important feature of the pyrochlore structure rele- 
vant to the centrál theme of this review is the fact that 
the 16c or 16d sites form layers stacked along the [111] 
direction as shown Fig. 4. From this perspective, the py- 
rochlore lattice is seen to consist of alternating kagome 
and triangular planar layers. As well, from the view- 
point of chemical bonding, the pyrochlore structure can 
be described as either an ordered defect fluorite (CaF2) 
or as two interpenetrating networks, one of composition 
^206, which is a network of corner-sharing metal-oxygen 
octahedra (topologically equivalent to that found in the 
perovskite structure) and the other of composition A20\ 
which forms a zig-zag chain through the large channels 
formed by the ^206 network. These models have been 
described in detail in several of the previous reviews and 
the discussion here focuses on a few important issues. 
While there is a topological relationship with the per- 
ovskites, the pyrochlore lattice is considerably more rigid. 
For example the B - O - B angle in which the O atom 
is shared between two octahedra is restricted to a very 
narrow range in pyrochlores, usually from about 127° to 
134°, with only little influence from the size of the A ion. 
On the other hand this angle can range from 180° to 140° 
in perovskites and there is a strong variation in the rádius 
of the A ion. The interpenetrating network description 
allows an understanding of the formation of so-called de- 
fect or "rattling" pyrochlores, such as KOS2O6, in which 
the entire A2O' network is missing and where the large 
K+ ion occupies the Sb site which is normally occupied by 
in the štandard pyrochlore structure. The disordered 
fluorite model is useful in cases where a disordered, de- 
fect fluorite structure phase competes with the ordered 
pyrochlore, such as for the A2ľí^2^i serieš where Aj B 
site mixing is obser ved. 



D. Phase Stability 

There are over 15 tetra-valent ions that can reside on 
the 5-site of a cubic pyrochlore at room temperature. 
Pyrochlore oxides of most interest here are formed with 
a rare earth ion in the A-site and transition or p-block 
elements in the 5-site. The rare earth elements form 
a serieš in which the trivalent state predominates and 
through which the ionic rádius decreases systematically 
with increasing atomic number (the so-called lanthanide 
contraction). Nonetheless, it is uncommon to flnd that 
the pyrochlore phase is stable for all A for a given B. 

The most efficient way to display the range of stability 
of pyrochlores in the space of the ionic radii is through a 
structure-field map or structure-stability map. Figúre 5 
is an abridged version of a comprehensive map published 
by Subramanian and Sleight (1993). First it is worth 
noting that the 5 = Sn serieš is the only one to form for 
all rare earths. There would also appear to be minimal 
and maximal ionic rádius ratios, A?^ jB^^ ^ which define 
the stability limits for pyrochlore phases of this type. 
For B ions with very small radii, such as Mn^+, the py- 
rochlore phase can be prepared only using high pressure 
methods. This is in strong contrast to perovskites of the 
type A^+Mn^+Os which can be prepared easily under 
ambient oxygen pressures. The first serieš stable under 
ambient pressure is that for 5 = V, but only for the 
three smallest members, Lu, Yb and Tm. The serieš can 
be extended using high pressure synthesis (Troyanchuk, 
1990). For the largest B ion, Pb^+, the smallest rare 
earth for which a stable pyrochlore is found appears to 
be Gd. Thus, an ambient pressure stability range defined 
in terms of rádius ratios, (RR) would extend 
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FIG. 6 Elements known to produce the (3+,4+) cubic py- 
rochlore oxide phase. 



FIG. 7 Variation of unit celí constant and physical properties 
of the serieš A2M02O7 with the rádius (Ah et al, 1989) 



between 1.36 and 1.71 with marginal stability expected 
as those limits are approached. There are exceptions, for 
example Pr2Ru207 is clearly stable with a RR of 1.82, 
yet Pr2Mo207 does not exist with a RR = 1.73. 

Also of relevance to this review are issues of defect for- 
mation in pyrochlores. This is an active area of research 
in materials science where defect pyrochlores which show 
high ionic conductivities are of interest. While the defect 
levels are often small and difficult to quantify via exper- 
iment, some computational guidelines have appeared re- 
cently (Minervini et al.^ 2002). The generál trend is that 
defect concentrations are predicted to increase as the RR 
approaches its lower limit, for example as the B rádius 
increases within any given serieš, and that the predicted 
levels are generally of the or der of 1% or less. 



VI. SAMPLE PREPARATION AND 
CHARACTERISATION 

Polycrystalline samples of pyrochlore oxides readily 
form by a štandard solid state process, where stoichio- 
metric proportions of high purity (> 99.99%) rare earth 
and transition metal oxides are calcined at 1350° C over 
sever al days in air with intermittent grindings. Elements 
known to exist in the cubic (3+,4H-) pyrochlore phase 
are highlighted in the periodic table shown in Fig. 6. 
Several serieš of compounds with particularly volatile el- 
ements, like t in, require excess materiál to account for 
losses during the heating processes and others require 
heating under special atmospheres to form the desired 
oxidation state. For example the molybdates require a 
vacuum furnace or a reducing atmosphere to form the 
Mo^+ state and the rare earth manganese (IV) oxide py- 
rochlores are not stable at ambient pressure and require 
a high pressure synthesis technique (Subramanian et a/., 
1988). Specialized growth conditions will be discussed 
below, as relevant. When high chemical homogeneity is 
required, the sol-gel method (Kido et al.^ 1991) has been 
successful. 



Crystals of a several oxide pyrochlores were grown in 
the 60's and 70's by flux methods (Wanklyn, 1968). The 
size and quality of the crystals grown by these methods 
vary, and it is often difficult to obtain voluminous crys- 
tals, especially for neutrón scattering experiments. The 
'floating zóne' method of crystal growth was first used 
by Gardner et al (1998) to produce large single crystals 
of Tb2TÍ207. Later the entire titanate serieš was pro- 
duced by others (Balakrishnan et a/., 1998; Tang et a/., 
2006; Zhou et al.^ 2007) and within the last few years, 
groups have started to grow many other pyrochlores in 
image furnaces including molybdates, ruthenates and 
stannates (Miyoshi et al.^ 2002; Taguchi et al.^ 2002; 
Wiebe et a/., 2004). Single crystals have also been pro- 
duced by several other methods including vapour trans- 
port (He et al.^ 2006), hydrothermal techniques (Chen 
and Xu, 1998) and KF flux (Millican et al, 2007). 

Powder samples of the ^-pyrochlores ("defect" py- 
rochlores) are formed at lower temperature 500° C) in 
an oxidizing atmosphere (Yonezawa et al.^ 2004). Small 
single crystals are also available (Schuck et al.^ 2006). 

Most oxide pyrochlores are cubic at room temperature 
with a lattice parameter between 9.8 Ä and 10.96 Ä (Sub- 
ramanian et a/., 1983). Of these, the titanates are prob- 
ably the most studied, followed by the zirconates, but 
many are still poorly represented in the literatúre. Both 
the titanate and zirconate pyrochlores were first reported 
by Roth (1956). Since then, over 1000 papers have been 
published on the two families. 



VII. METAL INSULATOR TRANSITIONS IN THE OXIDE 
PYROCHLORES 

Two types of metal-insulator (MI) transitions occur in 
the oxide pyrochlore family. First, are those that, as a 
function of a thermodynamic variable (e. g. temperature, 
magnetic field, pressure) change their transport proper- 
ties, for example Tl2Mn207 (Fujinaka et al.^ 1979). Sec- 
ondly, there are the serieš of compounds where the room 
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temperature character changes from metal to insulator 
as the rare earth ion changes, for example the molybdate 
(Greedan et al.^ 1986) or iridate serieš (Yanagishima and 
Maeno, 2001). There appears to be very little contro- 
versy over the first class of MI transition, however for 
the second type, the exact position of this transition is 
a topic of debate. For example, in the molybdenum py- 
rochlore serieš, studies of their bulk properties have indi- 
cated a strong correlation with the magnetism and elec- 
trical transport properties, i. e., the ferromagnets were 
metallic while the paramagnets were insulating (Ali et al, 
1989; Greedan et al.^ 1987). Indeed, the dependence on 
the lattice constant, ao, on the A^^ rádius showed a dis- 
tinct break at the MI boundary; see Fig. 7. In some 
subsequent studies however, the Gd phase is clearly in- 
sulating (Cao et al.^ 1995; Kézsmárki et al.^ 2004). The 
initial studies were carried out on polycrystalline sam- 
ples, prepared between 1300 - 1400° C and in at least one 
case, in a CO/CO2 "buffer gas" mixture which fixes the 
oxygen partial pressure during synthesis (Greedan et al.^ 
1986) Several subsequent studies have used single crystals 
grown by various methods above 1800° C, including melt 
and floating zóne growths (Kézsmárki et al.^ 2004; Morit- 
omo et al.^ 2001; Raju and Gougeon, 1995). While the 
polycrystalline samples have been fairly well character- 
ized, including elemental analysis, thermal gravimetric 
weight gain and measurement of the cubic lattice con- 
stant, ao, this is less true of the single crystals. The dif- 
ferences between poly and single crystalline samples can 
be monitored most simply using the unit celí constant 
as illustrated in Fig. 8 in which unpublished data (Raju 
and Gougeon, 1995) for a selection of single crystals of 
Gd2Mo207 are plotted. Note that as ao increases the 
samples become more insulating. From accurate struc- 
tural data for the powders and single crystals, it has been 
determined that the increase in ao correlates with an in- 
crease in the Mo - O dištance as shown Fig. 8. The most 
likely origin of this systematic increase is the substitution 
of the larger Mo^+ (0.69 Ä) for the smaller Mo'^+(0.64 Ä) 
which can arise from oxygen deficiency in the crystals re- 
sulting in the formula Gd2Mo2Í2a,Mo2+07-a;. Note that 
other defect mechanisms, such as vacancies on either the 
A or Mo sites would result in the oxidation of Mo^+ to the 
smaller Mo^+ (0.61 Ä) which would lead to a reduction in 
ao- Further evidence for this mechanism comes from the 
fact that the crystal with the intermediate celí constant 
in Fig. 8 was obtained from the most insulating crystal by 
annealing in a CO/CO2 buffer gas at 1400° C for 1-2 days. 
Thus, insulating single crystals of Gd2Mo207 are oxygen 
deíicient. This analysis is fully consistent with the report 
of Kézsmárki et al (2004) who find that "Gd2Mo207" is 
insulating but can be doped into the metallic state by 10 
percent Ca substitution on the Gd site. Ca-doping will 
of course add holes to the Mo states which compensate 
for the electron doping due to the oxygen deficiency in 
the "as grown" crystals. After considering all the stud- 
ies, we conclude that the MI boundary for stoichiometric 
742M02O7 is indeed between A = Gd and Tb. 



Gdj,Mo20^ crystals 
celi constants ^A) 
poly (metallic) = 10.337(1) 
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FIG. 8 Top: Comparison of electrical resistivity data for three 
single crystals of Gd2Mo2Í 2x^02^^7-0; as a function of unit 
celí constant, ao- Bottom: Variation of the Mo - O dištance 
with increasing celí constant for Gd2Mo2 07-cc crystals. 



Single crystals for the other ferromagnetic metals, A 
= Nd and Sm, are also reported and these are subject 
to the same oxygen deficiency and electron doping but, 
as they are farther from the MI boundary, the effects on 
the physical properties are less drastic. For example the 
transport, optical and magnetic properties, reported by 
Cao et al (1995J for Sm2Mo207 were for a crystal with 
ao = 10.425(1) A compared with the polycrystal with ao 
= 10.391(1) Ä. While the materiál is still metallic (albeit 
a poor metal), the ferromagnetic Curie temperature is 
reduced by 30 K relative to the powder, to 65 K. Note 
also the ZFC/FC divergence below 50 K which suggests 
a spin glassy component. Very similar results were re- 
ported for other A = Sm single crystals, ao = 10.4196(1) 
Ä, for example by Park et al. (2003) and by Moritomo 
et al (2001), which also show - 65 K and a ZFC/FC 
divergence below 20 K. 

The oxide pyrochlores can display a wide range of bulk 
properties including oxygen ion conductivity, supercon- 
ductivity, ferroelectricity and unusual magnetic behavior 
(e. g. spin liquid and spin ice). However these are inti- 
mately related to their stoichiometry, electrical and crys- 
tallographic properties, as shown in the previous para- 
graphs. Whenever possible, we encourage one to mea- 
sure, and report, these basic properties at room temper- 
ature or as a function of temperature. 
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FIG. 9 Properties of an oxygen deficient single crystal of 
Sm2Mo2 07-cc (a) Metallic behavior. (b) Inverse susceptibility 
showing a positive (FM) Curie-Weiss temperature. (c) Sus- 
ceptibility showing Tc = 65 K and a zero field cooled (squares) 
field cooled (triangles) divergence at 50 K (Cao et aL, 1995). 

Part III 

Experimental Results 

VIII. LONG-RANGE ORDERED PHASES 

Many magnetic oxide pyrochlores enter a long range or- 
dered state. These ordered phases are often promoted by 
single-ion anisotropy, anisotropic spin-spin interactions 
and/or further neighbor exchange interactions. Despite 
the fact that they develop long range or der, these sys- 
tems are still of significant interest, providing magnetic 
Systems with axial and planar symmetry, quantum order 
by disorder transition, polarised moments and partially 
ordered systems, for example. 



A. Gd2TÍ207 and Gd2Sn207 

Antiferromagnetically coupled Heisenberg spins on a 
pyrochlore lattice are expected to be highly frustrated 
because the exchange energy can be minimised in an in- 
finite number of ways. Classical (Canals and Lacroix, 
1998; Moessner and Chalker, 1998) and quantum calcula- 
tions (Canals and Lacroix, 1998, 2000) show that the sys- 
tém should remain in a collective paramagnetic state with 
short-range spin-spin correlations at any non-zero tem- 
perature. As discussed in Section IV, Raju et ai (1999) 
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FIG. 10 Left: False colour map of the temperature depen- 
dence of the magnetic diffraction from Gd2TÍ2 07. As the 
temperature is lowered the broad diffuse scattering centered 
at 1.2 sharpens into peaks. Right: Diffraction patterns 
at specific temperatures are shown along with tie marks in- 
dicating where Bragg reflections should be for the proposed 
ground state. At the lowest temperatures the m reflection 
has some intensity and this is highlighted in the inset (Stewart 
et a/., 2004a). 



reported an infinite number of degenerate spin configu- 
rations near the mean field transition temperature when 
dipolar interactions are considered. However, Palmer and 
Chalker (2000) showed that a four-sublattice state with 
the ordering vector /Cord=000 could be stabilized. More 
recent mean field theory calculations (Cépas and Shastry, 
2004; Enjalran and Gingras, 2003) found a higher tem- 
perature second order transition to a phase with three 
ordered sublattices and 1 disordered one. At a slightly 
lower temperature another transition was found to the 
Palmer-Chalker state. In contrast, Monte Carlo simu- 
lations by Cépas et al. (2005) find a single strongly first 
order transition to the Palmer-Chalker state. As we shall 
see below, the experimental situation is even more com- 
plex. 

The S = 7/2 Gd3+ ion, with its half-filled 4f-shell and 
zero orbital momentum is the best example of a Heisen- 
berg spin amongst the rare earths. The bulk susceptibil- 
ity for both Gd2TÍ2 07 and Gd2Sn2 07 follows the Curie- 
Weiss law with an effective moment close to the theo- 
retical limit of 7.94 (ib for the free ion (Bondah-Jagalu 
and Bramwell, 2001; Raju et ai, 1999) and a Curie-Weiss 
temperature of ^cw ^ —10 K. Both compounds however 
enter a magnetically ordered state at ~ 1 K, resulting in 
a frustration index of about 10. 

Early specific heat and AC susceptibility measure- 
ments showed that Gd2TÍ207 entered an ordered phase 
at 0.97 K (Raju et a/., 1999). A subsequent specific heat 
study showed two transitions in zero applied field, at 0.97 
and 0.7 K (Ramirez et al.^ 2002) and additional phase 
transitions, induced in applied fields. These phases were 
later confirmed by single crystal magnetisation work by 
Petrenko et al (2004). 

The náture of these ordered states is hard to deter- 
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mine by neutrón diífraction due to the high absorption 
cross section of gadolinium. However a small, enriched, 
sample of ^^^Gd2TÍ207 has been studied. In the initial 
work at 50 mK, a partiahy ordered noncohinear antifer- 
romagnetic structure, was proposed (Champion et al.^ 
2001). This unusual spin configuration had 3/4 of the 
Gd^+ spins ordered within kagome planes, separated by 
the triangular piane of Gd ions (see Fig. 4), which re- 
main disordered down to the lowest temperature. This 
model was later modified when the diffuse scattering and 
a new magnetic reflection at ^ ^ ^ , was observed and stud- 
ied (Stewart et al.^ 2004). The magnetic scattering from 
^^^Gd2TÍ207 is shown in Fig. 10. Very broad diffuse 
scattering is seen above 1 K. Below this temperature, 
Bragg peaks occur at the reciprocal lattice positions in- 
dicated by the tie marks (vertical bars) in the right panel. 
The lowest angle magnetic reflection, the very weak ^ | ^ 
peak at |Q| = 0.6 only gains intensity below the 
second transition. At all temperatures, diffuse magnetic 
scattering is observed which diminishes as the temper- 
ature decreases and the Bragg peaks grow. The dif- 
fuse magnetic scattering, centred at |Q| ~ 1.2 indi- 
cates that the correlation length of the disordered spins is 
^ 3.5 Ä, the nearest-neighbor dištance and not 7 Ä ap- 
propriate for the model proposed earlier by Champion 
et al. (2001). This picture of two different Gd moments is 
consistent with ^^^Gd Môssbauer experiments (Bonville 
et a/., 2003). Most interestingly, this is not the state 
proposed by Palmer and Chalker (2000). 

The isostructural stannate, Gd2Sn207, has received 
comparatively less attention. As mentioned above, the 
paramagnetic properties are similar to Gd2TÍ2 07 (Mat- 
suhira et al.^ 2002). However, only one transition, at 

I. K, and only one Gd moment has been observed in 
Gd2Sn2 07 (Bonville et a/., 2003). Recent neutrón diffrac- 
tion measurements by Wills et al (2006), shown in Fig. 

II, have confirmed that Gd2Sn207 undergoes only one 
transition down to the lowest temperatures and that 
the ground state is the one predicted theoretically for a 
Heisenberg pyrochlore antiferromagnet with dipolar in- 
teractions by Palmer and Chalker (2000). 

The difference in these two materials highlights the 
extréme sensitivity of the ground state to small changes 
in composition and of the consequential details in the 
spin-spin interactions in these highly frustrated materi- 
als. Wills et al. (2006) proposed that a different relative 
magnitude in the two types of third-neighbour superex- 
change interaction may be responsible for the different 
ground states in these otherwise very similar Gd com- 
pounds (Del Maestro and Gingras, 2007). 

To understand the magnetic ground states further and 
the true cause of the differences between Gd2TÍ2 07 and 
Gd2Sn2 07 one has to investigate the excitation spec- 
trum. Specific heat is a bulk próbe that often provides a 
reliable insight into the low lying excitations. The spe- 
cific heat, C^, of these two compounds just below the 
transition temperature is described well by an anomalous 
power law, C^ ex T^, that had been presumed to be asso- 
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FIG. 11 Neutrón powder diífraction pattern from Gd2Sn207 
at 0.1 K, after a nuclear background (1.4 K) data set was 
subtracted (Wills et a/., 2006). Inset: One of the 3 r7+ basis 
vectors that describe the data (see fit), consistent with the 
Palmer- Chalker model (Palmer and Chalker, 2000). 



ciated with an unusual energy-dependence of the density 
of states (Bonville et a/., 2003, a; Ramirez et a/., 2002; 
Raju et al.^ 1999). Electron spin resonance (ESR) (Has- 
san et al.^ 2003; Sosin et al.^ 2006, a), Môssbauer (Bertin 
et al.., 2002), muon spin relaxation (Bonville et al.., 2003a; 
Dunsiger et al., 2006; Yaouanc, et al., 2005) and neutrón 
spin echo (Ehlers, 2006) have provided many interest- 
ing results. In all these measurements, spin fluctuations 
were observed much below the ordering temperature of 
the compound, seemingly consistent with the unusual 
energy-dependence of the density of states invoked to 
interpret the specific heat measurements. Sosin et al. 
(2006, 2007) failed to confirm the unusually large planar 
anisotropy measured by Hassan et al. (2003), however 
they found a small gap ^ 1 K that coexists with a para- 
magnetic signál in Gd2TÍ207, however no paramagnetic 
spins were observed at low temperatures in Gd2Sn207. 
Recently Quilliam et al. (2007), remeasured the specific 
heat of Gd2Sn2 07 and found that the Cy oc behav- 
ior does not hold below 400 mK (see Fig. 12) and that 
C^ decreases exponentially below 350 mK, providing evi- 
dence for a gapped spin-wave spectrum due to anisotropy 
resulting from single-ion effects and long-range dipolar 
interactions (Del Maestro and Gingras, 2007). 

The current experimental understanding of these 
Heisenberg antiferromagnets is at at this time incom- 
plete. Gd2Sn207 appears to be a good example of the 
model considered by Palmer and Chalker (2000), and 
Gd2TÍ2 07 is perhaps a small perturbation from this. 
However, an explanation for the low lying excitations as 
probed by muon spin relaxation in both compounds still 
eludes theorists (Del Maestro and Gingras, 2007). Under- 
standing these ubiquitous, low temperature spin dy nam- 
ies in a relatively simple model magnet, like the Gd-based 
Heisenberg-like antiferromagnet on the pyrochlore lat- 
tice, would likely help elucidate the náture of the low tem- 
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Temperature (K) 

FIG. 12 Specific heat of Gd2Sn207 from Quilliam et al. 
(2007) (circles) and Bonville et al (2003) (diamonds). The 
power law (dotted line) and a conventional power 
law(dashed line) for a three dimensional antiferromagnet are 
plotted. The upturn seen below 150 mK results from the nu- 
clear electric quadrupole interaction (Quilham et a/., 2007). 

perature spin dynamics in other frustrated magnets. We 
discuss in Section XII.A some recent work on Gd2TÍ2 07 
in a magnet ic íield. 

B. Er2TÍ207 

The Er^+ single ion crystal field ground state in 
Er2TÍ207 is a Kramers doublet with a magnetic mo- 
ment of 3.8 fiB and 0.12 /íb perpendicular and paral- 
lel to the local (111) axis, respectively. Er2TÍ2 07, hke 
Yb2TÍ207, is a realization of an XY-hke systém on the 
pyrochlore lattice. However, while Yb2TÍ2 07 has a small 
ferromagnetic Curie-Weiss temperature (^cw ~ +0.75 
K) (Hodges et al, 2002), Er2TÍ207 has a large nega- 
tive Curie-Weiss temperature (^cw ~ - 22 K) (Champion 
et al.^ 2003), and is therefore a pyrochlore XY antiferro- 
magnet. The lowest-lying crystal field states of Er2TÍ2 07 
are at 6.38 meV (74.1 K) and 7.39 meV (85.8 K) above 
the ground states (Champion et al.^ 2003). 

Specific heat measurements show a sharp feature at 
approximately Tn ~ 1.2 K (Blôte et al, 1969; Champion 
et a/., 2003; Siddharthan et a/., 1999). Neutrón scatter- 
ing below Tn confirms the existence of a Néel long-range 
ordered state with propagation vector /Cord = 000. The 
neutrón scattering data strongly suggest a second order 
phase transition at Tn, with the scattered intensity / 
vanishing at Tn as /(T) oc (Tn — T)^^ with critical ex- 
ponent /3 ~ 0.33, characteristic of the three-dimensional 
XY model. 

As the transition is continuous, the systém is expected 
to order in only one of the irreducible representations of 
the Er^+ site representation. A refinement of the mag- 
netic structure finds that only the two basis vectors 
or í/js of Fs, or a superposition of both, is consistent with 




FIG. 13 Ground state spin configuration determined for 
Er2TÍ2 07 (Poole et a/., 2007). In a macroscopic sample six 
symmetry equivalent spin configurations will co-exist. 



the observed magnetic intensitites of the Bragg peaks. 
Results from a spherical neutrón polarimetry study find 
that the zero magnetic field ordered state is described 
almost entirely by with very little admixing from 7/^3 . 
The spin configuration corresponding to ip2 is shown in 
Fig. 13 (Poole et a/., 2007). The experimental observa- 
tions of (i) an ordering into the noncoplanar '02 structure 
via (ii) a second order phase transition are both surpris- 
ing. 

As discussed in Sections IV and VIII. A, isotropic ex- 
change plus dipolar interactions select by themselves a 
state with spins perpendicular to the local (111) axes. 
Hence, easy piane (X Y) anisotropy should further en- 
hance the stability of this state. This was explicitly 
shown in calculations (Del Maestro and Gingras, 2007) 
pertaining to Gd2Sn207 which does display a Palmer- 
Chalker ground state (Wills et al.^ 2006). However a 
Palmer-Chalker ground state for Er2TÍ207 is seemingly 
robustly rejected by recent neutrón polarimetry experi- 
ments (Poole et al.^ 2007). Since nearest-neighbor anti- 
ferromagnet ic exchange, dipolar interactions and the X Y 
náture of the magnetic moments should conspire together 
to produce a Palmer-Chalker ground state in Er2TÍ207, 
other interactions or effects must be at play to stabi- 
lize the experimentally observed '02 ground state (Poole 
et al, 2007). 

An early numerical study of the classical XY an- 
tiferromagnet on the pyrochlore lattice had found a 
strongly first order transiton into a long-range ordered 
state with propagation vector /Cord=000 (Bramwell et al, 
1994). This finding was subsequently confirmed by 
another study (Champion et a/., 2003; Champion and 
Holdsworth, 2004). In an earlier study, Bramwell et al. 
(1994) had only identified a macroscopic degeneracy of 
discrete ground states, while there is in fact a continu- 
ous degree of freedom for antiferromagnetically coupled 
XY spins on a single tetrahedron, resulting in a con- 
tinuous ground state degeneracy for a macroscopic sam- 
ple (Champion et al, 2003; Champion and Holdsworth, 
2004). Monte Carlo simulations of a classical (111) py- 
rochlore XY antiferromagnet find that the degeneracy is 
lifted by an order-by-disorder transition driven by ther- 
mal fluctuations with, most interestingly, the systém de- 
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veloping long range order in the same tjj2 state as found 
experimentally in Er2TÍ207. The Monte Carlo results 
can be rationalized on the basis of a classical spin wave 
calculation (Champion and Holdsworth, 2004). There, 
it is found that there are zero energy spin wave modes 
over planes in the Brihouin zóne for which the spin wave 
wavevector q obeys gr • a = O and q - (b — c) = O, where a, 
b and c are are the basis vectors of the primitive rhom- 
bohedral unit ceh of the pyrochlore lattice. 

On the basis of the results from Monte Carlo and clas- 
sical spin wave calculations, Champion and Holdsworth 
(2004) proposed that the ip2 in Er2TÍ2 07 is perhaps sta- 
bilized by zero-point quantum fiuctuations, arguing that 
their effects are captured by the classical spin wave ar- 
gument. However, problems remain with this interpret a- 
tion. As mentioned above, the transition at Tn ~ 1.2 K is 
second order in the experiment while simulations show a 
strong first order transition (Champion and Holdsworth, 
2004). Also, the density of zero energy classical spin wave 
States is a constant, a result that is incompatible with the 
experimental specific heat data which shows a tem- 
perature dependence below 1 K. Inelastic neutrón scat- 
tering reveals no evidence for such low lying excitations 
required to generate the temperature dependence of 
the specific heat. The experimental results and mod- 
els presented above suggest a well established long-range 
ordered ground state for Er2TÍ2 07. Yet, muon spin re- 
laxation data reveal that below Tn, the muon polariza- 
tion relaxation rate remains large (2 x 10^ Ms~^) and 
essentially temperature independent down to 20 mK, in- 
dicative of a dynamic ground state (Lago et al.^ 2005). 
This behavior reminds one of low-temperature muon 
data in several pyrochlores including the "ordered" Gd- 
pyrochlores and the "dynamic" Tb2TÍ2 07. 

Unlike Er2TÍ207, Er2(GaSb)07 does not appear to or- 
der down to 50 mK (Blôte et a/., 1969). It is not clear at 
this stage whether the randomness on the B-sites leads 
to a disordered and glassy phase, i. e. an XY-like spin 
glass. Results on Er2Sn2 07 will be discussed in Section 
XI.C. 



C. Tb2Sn207 

Until recently the magnetic and electrical properties 
of the stannates had received little attention. This is 
probably due to the slightly more problematic synthesis 
and the lack of large, high quality, single crystals. 

Like many frustrated magnets, the high temperature 
(T > 40 K) inverse susceptibility is very linear. In 
Tb2Sn2 07, this yields a Curie-Weiss temperature of ap- 
proximately ^cw ~ — 12 K and an effective moment of 
9.8 jiíB^ which agrees well with the value of 9.72 jj^b 
for the ^Fe ground state of Tb^+ (Bondah-Jagalu and 
Bramwell, 2001). Reminiscent of the cooperative para- 
magnet Tb2TÍ207 (see Section XI), the inverse suscepti- 
bility deviates downwards at ~ 30 K upon cooling. How- 
ever, on cooling further, the systém enters a frozen ferro- 
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FIG. 14 Magnetic scattering from Tb2Sn207 at 100 mK. 
Inset: the proposed "ordered" spin ice structure with spins 
canted off the (111) by 13° (Mirebeau et a/., 2005). The broad 
scattering at the base of the Bragg reflections is indicative of 
spin clusters which reach ^ 200 Ä at base temperature. 



magnetic state at 0.87 K. This was first reported by Mat- 
suhira et ai (2002) from low temperature susceptibility 
measurements and then by Mirebeau et al. (2005) who 
determined the magnetic structure from neutrón diffrac- 
tion measurements. In this work by Mirebeau et al. 
(2005), a crossover from antiferromagnetic correlations 
to ferromagnetic correlations was observed at 2 K. This 
was followed, at a lower temperature, by a transition into 
an ordered state and the appearance of Bragg peaks, see 
Fig. 14. Rietveld analysis determined that the spins are 
canted by 13.3° off the local (111) directions with an or- 
dered moment of 5.9 /i^, significantly reduced from the 
free ion moment of 9 /i^, presumably due to crystalline 
electric field effects. 

Ferromagnetically coupled spins on the pyrochlore lat- 
tice fixed along the local (111) directions gives rise to a 
highly degenerate state that has become known as the 
spin ice state (see Section X). By comparing the Tb^+ 
moment value deduced from neutrón diffraction (5.9 /j. b) 
to the smaller value deduced from the nuclear specific 
heat (4.5 jub) Mirebeau et al. (2006) suggested that slow 
collective fiuctuations, on the time scale of 10~^-10~^s, 
coexist with the ordered state at the lowest tempera- 
tures (Mirebeau et al.^ 2006). 

Subsequently, two groups led by Dalmas de Réotier, 
et al. (2006) and Bert et al. (2006) have suggested that 
the ground state is not státie. They claim that the ab- 
sence of any oscillations, even at the shortest times, in the 
muon spin relaxation spectrum precludes the existence of 
a státie moment. Dalmas de Réotier, et al. (2006) argue 
that both neutrón and muon data are consistent with 
a characteristic spin relaxation time of ^ 10~^^s whilst 
Bert et al. (2006) suggest that the dynamics result from 
fiuctuations of clusters of correlated spins with the or- 
dered spin ice structure. On the other hand, the results 
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FIG. 15 Top: The temperature dependence of the suscepti- 
bility for Y and Lu ruthenate. Both data sets show the Ru'^^ 
moments ordering at ~ 75 K. Bottom: Heat capacity mea- 
surements over the entire rare earth ruthenate serieš showing 
how the Ru^^ ordering temperature increases as the size of 
the rare earth ion decreases (Ito et a/., 2001). 



from the neutrón spin echo study of Chapuis, et al (2007) 
suggest no spin dy namies below 0.8 K. 

Mirebeau et al (2006) have suggested that the ^3% 
expansion in the lattice constant between Tb2Sn2 07 and 
Tb2TÍ2 07 results in the dipolar interaction overcoming 
the weakened exchange interaction allowing Tb2Sn2 07 to 
freeze. Clearly more experiment al work is needed on this 
compound to determine the exact náture of the ground 
state. Recently, from inelastic neutrón scattering stud- 
ies, Mirebeau et al. (2007) found that the composition of 
the first two crystal field levels of Tb2Sn2 07 are inverted 
when compared to those of Tb2TÍ2 07. This cannot be 
explained by simple calculations of the local crystal field 
environment and more work is needed. 



D. A2RU2O7 {A = Y, Gd, Dy, Ho, Er and TI) 

All but the largest, A = La and Ce, rare-earth ruthen- 
ates as well as A = Bi and A = TI form as cubic py- 



rochlores (Bertaut et a/., 1959; Subramanian et a/., 1983; 
Yoshii and Sato, 1999). The lattice parameters vary from 
10.07 to 10.355 Ä across the rare earth serieš with the A = 
Bi and TI compounds at the high end of this range. Crys- 
tals of several ruthenates have been made by hydrother- 
mal, vapour transport and floating zóne methods (Sleight 
and Bouchard, 1972; Zhou et al.^ 2007). Historically, the 
ruthenium pyrochlore oxides have been extensively stud- 
ied for their novel conductivity (Carcia et a/., 1982; Pike 
and Seager, 1977) and catalytic activity (Horowitz et al.^ 
1983), but until recently, very little was known about 
their magnet ic properties. Most of the ruthenates are 
semiconductors at room temperature with hopping en- 
ergies of ~ 0.1 e V, but the Bi based sample is metallic 
with Pauli paramagnetism. Ru^+ (low spin 4d^) has an 
S' = 1 moment and the Ru sublattice appears to order 
between 70 and 160 K, depending on the rare earth (see 
Fig. 15). However, Ocw for Y2RU2O7 is reported to be 
-1250 K although it is not clear that this was obtained in 
the paramagnetic regime (Gurgul et a/., 2007). Studies of 
diamagnetic ^4 = Y and Lu provide evidence of the mag- 
netic náture of the Ru^+ ion. These indicate a magnetic 
transition at 75 K, leading to a high frustration index, / 
^17. The magnetic susceptibilities measured under zero 
field cooled (ZFC) and under field cooled (FC) conditions 
show a different temperature dependence (see Fig. 15). 
This is consistent with the observation of a canted an- 
tiferromagnetic structure from neutrón diffraction (Ito 
et a/., 2001; Kmiec et al.^ 2006). Inelastic neutrón scat- 
tering measurements by van Duijn et al. (2007) reveals 
the development of a large gap in the excitation spectrum 
below the ordering temperature. 

Er2Ru2 07 and Gd2Ru207 have been studied by 
powder neutrón diffraction (Taira et al, 2003) and 
Môssbauer (Gurgul et al.^ 2007). In both cases, the Ru 
sublattice orders at the temperature indicated by the 
peak in the specific heat plotted in Fig. 15. This long 
range ordered state is a /Cord=000 type. This is also true 
for Y2RU2O7 (Ito et al, 2001). In the Er and Gd ruthen- 
ates, the rare earth moment is polarised by the ordered 
Ru- sublattice, via the rare earth - Ru^+ exchange inter- 
actions. This reveals itself in a small but finite moment 
on the A-site as soon as the ruthenium orders. This mo- 
ment grows with decreasing temperature before the rare 
earth ion orders in its own right at 10 K (Er) and 40 K 
(Gd). A refinement of the low temperature (3 K) neu- 
trón diffraction data from Er2Ru207 suggests that the 
magnetic moments, both Er and Ru, are in a piane ori- 
ented towards the same [100] direction, but antiparallel 
to each other. The refined moments for Ru^+ saturate 
at the theoretical limit for the S=l ion (2 Jííb) at 90 K, 
but at 3 K the moment on the Er site is only 4.5 /i^, 
50% of the predicted value, and is probably reduced by 
crystalline electric field effects, although the signál has 
not yet saturated at the lowest temperature studied (see 
Fig. 16). 

Since Bansal et al. (2002) reported a possible spin- 
ice ground state (see Section X) in H02RU2O7 and 
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FIG. 16 Left: The temperature dependence of the Gd mo- 
ment from Gd-Môssbauer (Gurgul et aL, 2007) on Gd2Ru207. 
Right: The temperature dependence of the Ru^^ and Er^^ 
moment s in Er2Ru207 as measured by neutrón diffrac- 
tion (Taira et a/., 2003). 
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FIG. 17 Specific heat of polycrystalhne H02RU2O7 as a func- 
tion of temperature around the Ru^^ (inset) and Ho^^ (main) 
ordering temperatures (Gardner et a/., 2005). 



Dy2Ru2 07, there have been several studies of the mag- 
netic behaviour of H02RU2O7 (Bansal et al.^ 2003; Gard- 
ner et al.^ 2005; Wiebe et al.^ 2004). In this compound, 
unHke the other heavy rare earth ruthenates, a small dif- 
ference is seen between FC and ZFC when the rare earth 
ion orders. This happens at 1.4 K when the Ho^+ sub- 
lattice enters a long-range ordered state, unhke a true 
spin ice materiaL Lower temperature (T < 20 K) spe- 
cific heat (see Fig. 17) also highhghts this transition. As 
seen in a number of geometrically frustrated systems, a 
broad feature (centered here at ~3 K) precedes the sharp 
lambda-hke transition at a shghtly lower temperature. 
These features are related to the build-up of short- and 
long-range correlations between Ho spins, respectively. 
The ^ordereď spin ice has the holmium moments canted 
oíf the local (111) axes presumably due to the local fields 
generated by the ordered Ru^+ sublattice. The ordered 
structure of Gd2Ru207 is similar to H02RU2O7. In both 
cases the full moment of the rare earth ion was not mea- 
sured at the lowest temperatures investigated (Gurgul 
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FIG. 18 Neutrón spin echo spectra from H02RU2O7 at 
IQI = 0.4 Ä"^ as measured on NG5NSE (left) and the de- 
pendence of the relaxation time on temperature (right) from 
two different positions in reciprocal space. The Unes through 
the spectra on the left are a simple exponential, where the 
characteristic relaxation time is used in the Arrhenius plot on 
the right (Gardner et al, 2005). 



et al, 2007; Wiebe et al, 2004). 

The magnetic entropy associated with H02RU2O7 at 
these low temperatures reaches neither i? In 2 (expected 
for a completely ordered doublet systém) nor R{hi2 — 
^ In |) (expected for spin ice). In fact, the magnetic en- 
tropy of H02RU2O7 (Gardner et al.^ 2005) resembles that 
of the dipolar spin ices in an applied field of 1 T. Several 
neutrón scattering (Gardner et al.^ 2005; Wiebe et al.^ 
2004) and AC susceptibility (Gardner et al.^ 2005) experi- 
ments have investigated the spin dynamics in H02RU2O7. 
The náture of the Ho-moments is governed by the crys- 
talline electric field levels which are slightly perturbed 
from those measured in H02TÍ2O7 (Rosenkranz et al.^ 
2000) especially after the ordering of the Ru^+ sublat- 
tice (Gardner et al.^ 2005). This is most apparent in 
neutrón spin echo measurements where a 35 K change 
in the gap to the first excited state is observed once the 
Ru^+-ions order (Gardner et al, 2005) (see Fig. 18). 

While TI2RU2O7 is cubic with a room temperature lat- 
tice parameter of 10.18 Ä, when it is cooled below 120 K, 
the structure transforms to an orthorhombic structure. 
Concomitantly, a metal-insulator transition occurs and 
the resistivity changes by 5 orders of magnitude. A sud- 
den drop in the susceptibility also occurs at this tempera- 
ture (Lee et al.^ 2006). These data have been interpreted 
as the three-dimensional analogue to the more commonly 
known one-dimensional S = 1 Haldane gap systém (Hal- 
dane, 1983). Specifically, it is believed that the structural 
phase transition creates effective 5 = 1, chains of Ru-ions 
which do not dimerize. These chains are arguably created 
due to the orbital ordering of the Ru 4d electrons. 

Three (2-h,5+) pyrochlores of interest are 
Cd2Ru207 (Wang and Sleight, 1998), Ca2Ru207 (Mune- 
naka and Sato, 2006) and Hg2Ru207 (Yamamoto et al.^ 
2007), with an S = 3/2 Ru^+ ion and non-metallic 
transport properties. Hg2Ru2 07, like TI2RU2O7, has 
a metal-insulator transition that is concomitant with a 
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Reference (1988) prepared the phases A = Sc, Y, In, TI and Tb- 
(Troyanchuk et al iggé")*^ with pressures between 5-8 GPa and temperatures 
;JOqp - 1500°C. Subramanian et al (1988) were able to 



1500^ 
rdpare A = 



C. Subramanian et 

Y and Dy-Lu using much lower pressures 



(Greedan et al, 1996a' 

(Fujinaka et a/., 1979^^ j ^^gj^pg^^g^i^^^^gg with a hydrothermal method in sealed 
(Troyanchuk et al, 198§Xn tubes including NaClOa, NaOH and H2O at 3 kbar 
(Subramanian et al, 19^%]^ GPa) and 500°C. Greedan et al (1996a) synthesized 



(Reimers et al, 1991)A =:= Sc using a tetrahedral anvil press at 850°C and 60 
(Shimakawa et al, 199^bar (6 GPa). Subsequently, Shimakawa et al (1999) 
(Troyanchuk et al 198&f^^^ ^ isostatic press apparatus at 1000°C - 1300°C 
ÍSubramanian et al 19^\^^ ^^^^ ^'^ ^^^^ (0.04 GPa) for the serieš A = In, Y 

anťl Lu while the ^4 = TI materiál required 2.5 GPa and 

(Shimakawa et al, 199^]3^^o^ ^ piston-cylinder apparatus (Shimakawa et al, 



TABLE II Collected values for unit celí constants, ao, mea- 
sured Curie temperatures, Tc, and Curie-Weiss temperatures, 
^cw, for A2Mn207 phases with A — Sc, Y, Lu, In and TI. 



structural phase transition. Single crystals of Ca2Ru2 07 
have been grown by the hydrothermal method (Mune- 
naka and Sato, 2006) and the effective paramagnetic 
moment is only 0.36 /i^, about one or der of magnitude 
smaller that the theoretical value for the S = 3/2 ion. 
At 23 K, there is clear irreversibility between FC and 
ZFC susceptibilities indicative of a spin glass transition 
and reminiscent of Y2M02O7 (Greedan et al, 1986). 



E. ^sMtisOr (A = Sc, Y, Tb - Lu and TI) 

The rare earth manganese (IV) pyrochlores are not sta- 
ble at ambient pressure at any temperature and must be 
prepared using high pressure methods. The earliest pub- 
lished report is by Fujinaka et al (1979) who synthe- 
sized the pyrochlores A2B2O7 with A = Y, TI and B = 
Cr, Mn using temperatures in the range 1000 - 1100°C 
and pressures from 3-6 GPa. Reports of the high pres- 
sure synthesis of several A2Mn207 phases appeared from 
the groups of Troyanchuk et al (1988) and Subramanian 
et al (1988) about ten years later. Troyanchuk et al 
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LU2V2O7 


83 


74 


1.1 


YTÍO3 


33 


29 


1.1 



(Troyanchuk et al, 198§^c 9)^ 

(Raju et al, 1994) ^ 'j-^^ ^j^it celí constants from all known preparations of 
(Shimakawa et al, 1999^2^^^2 07 phases along with preliminary magnetic char- 
(Fujinaka et al, 1979]^ct[írization data are collected in Table II for A = Sc, Y, 
(Raju et al 1994) "^^^ where only the Mn-sublattice is mag- 

ŕ Shimakawa et al 199§f^ ^' ^^^^^ variety of preparative conditions, the 

agreement among the various groups for these key pa- 

rameters is excellent, suggesting that sample to sample 
compositional variation is small. As seen in the case of 
the molybdate pyrochlores (Section VII), both the unit 
celí constant, ao, and Tc are very sensitive to compo- 
sition. All of these materials appear to be ferromag- 
nets, consistent with the observation of a positive Curie- 
Weiss temperature and an apparent Curie temperature 
as shown for example in Fig. 19 for Tl2Mn207 (Raju 
et al, 1994). Note, however, that there appear to be two 
separáte classes of ferromagnetic materials, if sorting is 
done by the ratio Oqw/Tc. The compounds A = T\ and 
In have Tc values of 120 K and Ocw/Tc ratios of about 
1.25-1.3. Of course mean field theory gives Tc = Oqw for 
a ferromagnet. Empirically, for most ferromagnets this 
ratio is larger than 1, with typical values between 1.1 and 
1.2 as shown in Table III for some selected ferromagnetic 
insulators. As well, one can estimate a theoretical ra- 
tio calculated using the Rushbrooke-Wood relationship 
for Tc and the mean field theory for ^cw, assuming the 
same exchange constant J (nearest-neighbors only) and 
S = 3/2, appropriate for Mn^+, which gives Tc/Oqw = 
1.4 (Rushbrooke and Wood, 1958). On the other hand 
the A = Sc, Y and Lu compounds show corresponding 
ratios within the range 3-4, values well outside those 
seen for typical ferromagnets and well outside the ex- 
pected theoretical limits. This suggests that these phases 
present a much more complex situation than that of a 
simple ferromagnet. Interestingly, when the A element 
has a magnetic moment {A = Tb - Yb), the properties, 
especially the Ocw/Tc ratio, return to typical simple fer- 
romagnetic levels just above 1 as shown in Table IV. 



TABLE III Comparison of the Curie-Weiss and ordering tem- 
peratures for selected ferromagnetic insulators. 



1. Tl2Mn207 

Tl2Mn207 
lowing the 



has attracted 
discovery of 



considerable attention fol- 
giant magneto-resistance 



(GMR) (Cheong et al, 1996; Shimakawa et al, 1996; 
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(Greedan et a/., 1996) 



TABLE IV Collected values for unit celí constants, ao Curie 
temperatures, Tc, and Curie- Weiss temperatures, ^cw, for 
A2Mn207 phases with A = Th - Yb. 



Subramanian et al, 1996). While this topic is outside 
the main area of this review, having little to do with ge- 
ometrie frustration, some eomments are in order. Fig- 
úre 19 shows the evidence for this effect (Shimakawa 
et al, 1996). It was soon realized that the mechanism 
for this GMR must be different from that for the well- 
studied perovskite manganates (Schiffer et al.^ 1995). In 
the latter, double exchange between Mn^+ and Mn^+ is 
in part involved in both the metallic behavior and the 
ferromagnetism. From a combination of accurate mea- 
surements of Mn - 01 distances and core level spectro- 
scopies, it was determined that only Mn^+ is present in 
Tl2Mn2 07 (Kwei et al.^ 1997; Rosenfeld and Subrama- 
nian, 1996; Subramanian et al.^ 1996). The ferromag- 
netism arises from superexchange, not double exchange, 
and the metallic properties are due to accidental over- 
lap between the TI 6s band and the Mn 3d band. The 
coupling between magnetic and transport properties re- 
sults from abnormally strong incoherent scattering of the 
conduction electrons due to spin fluctuations which ac- 
company the FM orderin GMR thus arises from the large 
field dependence of Tc. Interestingly, doping of the TI site 
with In or Sc, the parent compounds of which are both 
insulating, greatly enhances the GMR effect as does dop- 
ing with Tl-vacancies or even Cd (Alonso et al.^ 2000; 
Cheong et al.^ 1996; Ramirez and Subramanian, 1997; 
Velasco et al, 2002). 

Neutrón diffraction confirmed the ferromagnetic long 
range order for this materiál, finding an ordered mo- 
ment of 2.91(4) jiiB per Mn^+ ion, in excellent agree- 
ment with the 3.0 ja b expected for an S = 3/2 state. 
Small angle neutrón scattering is fully consistent with 
the formation of domain walls and spin waves below Tc 
= 123.2(3) K (Lynn et al, 1998; Raju et al, 1994). These 
results are of interest when compared to data from the 
insulating materials A = Y, Ho and Yb, below. 
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FIG. 19 Top: Magnetic properties of Tl2Mn207 showing ev- 
idence for ferromagnetic order in 0.005 T (Raju et aL, 1994). 
Bottom: Evidence for a giant magneto resistance effect in 
Tl2Mn207 (Shimakawa et a/., 1996). 



2. Y2Mn207 

As shown by Reimers et al. (1991), this materiál ex- 
hibits many of the features of a ferromagnet, including an 
apparent Tc ~ 16(1) K from DC magnetization studies 
and saturation with applied magnetic field at low temper- 
atures (see Fig. 20). Note that saturation is not reached 
even at /íq^ = 4 T and that the moment is close to 2 /i^ 
rather than the 3 /íb expected for an 5 = 3/2 ion. 

Behavior quite atypical of ferromagnetism is found in 
the AC susceptibility and the heat capacity is shown 
in Figs. 20 and 21. First, while x' does indeed show 
a sharp increase below 16 K, the dominánt feature is 
a broad maximum centered at about 7 K which shows 
a strong frequency dependence, typical of spin glassi- 
ness. Controversy exists concerning measurements of the 
specific heat. Reimers et al (1991) report no lambda- 
type anomaly near 16 K and a linear temperature depen- 
dence below 7 K, behavior again typical of spin glasses. 
They also reported that about 60 percent of the total 
entropy of i? In 4 (attained by 100 K) is removed above 
16 K. On the other hand Shimakawa et al. (1999) find 
a lambda anomaly at 15 K which they have interpreted 
as evidence for long range ferromagnetic order. Figúre 
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FIG. 20 Bulkpropertiesof Y2Mn207. Top left: DC suscepti- 
bility showing Tc = 16 K in 0.15 mT (circle), 0.56 mT (square) 
and 10 mT (triangle). Bottom left: Magnetization versus ap- 
plied field at various temperatures (K) from top to bottom: 
1.8, 5, 7.5, 10, 15, 20, 25, 30, 35, 40, 45 and 50 K. Right: AC 
susceptibility. (C) x — 20 Hz. The inset shows the x maxi- 
mum for various frequencies, top to bottom: 20, 100, 200 and 
1000 Hz. (D) x'' cit various frequencies, from top to bottom: 
20, 40, 80, 100, 200, 500 and 1000 Hz. The curves are each 
shifted by -0.1 from the preceding one (Reimers et a/., 1991). 



21 compares both results and those for a non-magnetic, 
lattice matched materiál, Y2Sn207. Although the agree- 
ment in the data between the two groups is good, apart 
from the weak lambda anomaly, Fig. 21 shows consider- 
able magnetic heat capacity both above and below the 
lambda feature when compared to Y2Sn2 07. Clearly, 
this is not consistent with a simple ferromagnetic transi- 
tion. As well, even the observation of a lambda anomaly 
in the specific heat is no guarantee of a phase transition 
to a conventional long-range ordered state as has been 
demonstrated for other pyrochlores like Yb2TÍ207 and 
the garnet, Yb3Ga50i2 (Blôte et al.^ 1969; Dalmas de 
Réotier, et al, 2003; Hodges et al, 2002). While the re- 
sults of Fig. 21 cannot be easily reconciled, the claim by 
Shimakawa et al. (1999) that the heat capacity demon- 
strates, unequivocally, an ordered state is not substanti- 
ated and should be investigated further. 

Neutrón scattering data from Reimers et al. (1991) 
and Greedan et al. (1991) also provide evidence for a 
quite complex ground state. Elastic scattering data show 
mainly broad features below 50 K (see Fig. 22). En- 
hanced magnetic scattering forms below 15 K at nu- 
clear Bragg positions corresponding to the 111 and 222 
at 13.9° and 28.0° respectively. These data were fur- 
ther analyzed via Fourier transformation to give the real 
space spin-spin correlation function which suggested that 
the nearest-neighbor exchange was weakly antiferromag- 
netic while further neighbor exchange is strongly ferro- 
magnetic. This is consistent with the observed positive 



Curie-Weiss temperatures and the apparent failure to or- 
ier in a simple ferromagnetic ground state. 
Small angle neutrón scattering (SANS) data from 
reedan et al. (1991, 1996) are also unusual and in fact 
inprecedented. These are compared with correspond- 
ng data for ferromagnetic Tl2Mn2 07 in Fig. 23 (Lynn 
et al, 1998; Raju et al, 1994). At |Q| = 0.030 Ä"^ the 
■erromagnet (Tl2Mn207) shows a decrease in scattering 
ntensity with decreasing temperature below Tc which is 
iscribed to scattering from spin waves which diminish in 
ntensity as the long range order is established and the 
^luster size moves out of the SANS window. On the other 
land, the behavior of Y2Mn207 is the opposite, showing 
m increase below the apparent ordering temperature of 
15 K, indicating that the population of subcritical clus- 
ters does not diminish significantly to low temperatures. 
In addition the full range of data could not be fitted to 
a simple Lorentzian as is normally the rule for ferromag- 
nets, as for example in Tl2Mn207 (Lynn et al, 1998), but 
only to a Lorentzian plus Lorentzian-squared (L+L^) law 
which involves two correlation lengths, and ^2, 



S{Q) = 



A 



B 



(Q + 1/6) (Q + 1/6)2 



(5) 



The L + L^ law is often found for systems in which 
there is competition between ferromagnetic order and 
random-field disorder and where ^1 is associated with 
the ferromagnetic correlations and ^2 with the random 
fields (Aharony and Pytte, 1983; Arai et al, 1985; Rhyne, 
1985). As shown in Fig. 24, for Y2Mn207 the two corre- 
lation lengths differ by nearly an order of magnitude over 
the temperature range investigated. ^1 shows a buildup 
through a rounded maximum to a value near 20 Ä. Sim- 
ilar broad maxima are found in nearly ferromagnetic, 
metallic glasses (Rhyne, 1985; Rhyne and Fish, 1985). 
^2 quickly diverges to a resolution limited value below 
15 K. The most straightforward interpret ation is that ^1 
measures the temperature evolution of ferromagnetic cor- 
relations that never realize a long-range ordered state due 
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FIG. 21 Left: Comparison of the total heat capacity for 
Y2Mn2 07 reported by Reimers et al. (1991) (black) and Shi- 
makawa et al. (1996) (red) and the diamagnetic materiál, 
Y2Sn207 (blue). Right: The magnetic component of the heat 
capacity. Note the absence of a sharp peak and the linear de- 
pendence at low temperatures (Reimers et al, 1991). 
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FIG. 22 Neutrón diffraction data for Y2Mn207 showing dif- 
ference plots (data at 200K subtracted) at various tempera- 
tures. Note the buildup of broad peaks at the 111 and 222 
Bragg positions (Reimers et aL, 1991). 
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to the intervention of random fields as monitored by ^2- 
In this case, both tendencies have onsets near the same 
temperature. Note that there is no anomaly near 7 K, 
the maximum in the real part of the AC susceptibihty. 
While a detailed understanding of the SANS results for 
this materiál is not yet available, there is httle in these 
data to indicate that Y2Mn207 behaves as a conventional 
ferromagnet. 



3. Ho2Mn207 and Yb2Mn207 

Among the remaining 742Mn207 pyrochlores only 
those with A = Ho and Yb have been studied in de- 
tail (Greedan et al.^ 1996). As seen from Table III, the 
6cw/Tc ratio is within the range typically seen for fer- 
romagnets and both show positive Curie-Weiss temper- 
atures which indicates ferromagnetic rather than ferri- 
magnetic correlations between the two sublattices. The 
apparent Tc values are significantly higher than those for 
A = Y, Lu or Sc indicating that the A - Mn exchange has 
a strong influence on the ordering temperature. Unlike 
Y2Mn207, both materials show magnetic saturation at 
5 K for modest applied fields /íqH > 2 T with values of 
12.4 jiíB and 9.2 jj^b per formula unit for A = Ho and Yb, 
respectively. Assuming ferromagnetic A - Mn coupling 
and that the Mn sublattice moment saturates with the 
full spin-only value of 3.0 iib per Mn^+, the A sublattice 
saturation moment s are 3.2 jn b and 1.6 jj. b per Ho^+ and 
Yb^+ ion, respectively. These are considerably smaller 
than the full saturation moments, g^J^ for free Ho^+ 



FIG. 23 Comparison of the temperature dependence of the 
total small angle neutrón scattering (SANS) intensity data 
for the ferromagnet Tl2Mn207 and Y2Mn207 at |Q| = 0.03 
(Greedan et al, 1996; Raju et al, 1994). Similar results 
for Tl2Mn207 were obtained by Lynn et al (1998). 



(10.0 11b) and Yb^+ (4.0 (Ib)-, and indicate the strong 
infiuence of crystal fields on the ground state. A simi- 
lar Yb moment is found in the isostructural compound, 
Yb2V2 07 (Soderholm et al, 1980). 

In addition to DC susceptibihty and magnetization, 
AC susceptibihty, heat capacity and neutrón scattering 
data are available but these do not lead straightforwardly 
to a consistent interpretation. For example, the real 
part of the AC susceptibihty data shown in Fig. 25 in- 
dicates divergent behavior for the two materials. While 
there is an apparent Tc at 38(1) K with a broad maxi- 
mum just above at 30 K for both materials, the data for 
Ho2Mn207 show a strong frequency dependence while 
those for Yb2Mn207 do not. Thus, the Ho2Mn207 phase 
appears to show a so-called re-entrant spin glass state 
where glassy behavior develops at a temperature below 

Heat capacity (see Fig. 25), and neutrón scattering (see 
Fig. 26), data seem to suggest a very complex magnetic 
ground state for both materials. No lambda anomaly 
consistent with true long range order has been observed 
but resolution limited magnetic Bragg peaks are seen for 
both materials, the temperature dependence of which is 
order parameter like and consistent with Tc = 38(1) K 
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FIG. 24 Temperature dependence of the correlation lengths 
^1 and ^2 for Y2Mn207, Ho2Mn207 and Yb2Mn207 (Greedan 
et a/., 1996). 
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FIG. 26 Temperature dependence of magnetic Bragg peaks 
for Yb2Mn2 07 and Ho2Mn2 07 showing order parameter like 
behavior and an apparent T c of 38(1) K for both (Greedan 
ei al, 1996). 



for both compounds. However, the magnetic structure 
for Yb2Mn2 07 and Ho2Mn2 07 is stih unknown. SANS 
data are similar to those for Y2Mn2 07 in that the L + 
law holds and the two correlation lengths have diíferent 
temperature dependences (see Fig. 24). Interestingly, the 
B coefficient in Equ. 5, which measures the contribution 
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FIG. 25 Bulk properties of Yb2Mn207 and Ho2Mn207. Left 
column is the real part of the AC susceptibility for and 
right column is the magnetic component of the heat capac- 
ity (Greedan et aL, 1996). 



from the term associated with the ^2 correlation length 
shows a temperature dependence akin to that of an order 
parameter for all three mater ials. 

More work on the ordered manganese pyrochlores is 
needed to understand fully their magnetic ground states. 
The role of random fields arising from the geometrically 
frustrated Mn sublattice, but whose microscopic origin 
is not understood, play an important role in the deter- 
mination of the magnetic ground state which appears to 
be rather inhomogeneous. The A = Y, Ho and Yb com- 
pounds do not behave like simple ferromagnets and stud- 
ies of their spin dynamics are warranted. In addition, 
diffraction studies using modern Instruments should be 
performed on A = Ho and Yb to determine the ordered 
component of the ground state. 



F. A2lr207 {A = Nd - Yb) 

In the rare earth iridium(IV) pyrochlores, iridium has 
an electronic configuration 5d^ which is expected to be 
low spin, resulting in an 5 = 1/2 systém and this serieš 
is t hus of considerable interest. The pyrochlore struc- 
ture is reported to be stable for A = Fľ - Lu, one of the 
widest stability ranges for any pyrochlore systém (Subra- 
manian and Sleight, 1993). This serieš of compounds not 
only shows a metal-insulator (MI) transition across the 
rare earth serieš at room temperature, like that discussed 
in Section VII for the molybdates, but sever al individual 
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FIG. 27 Top: Resistivity data for the A2IT2O7 serieš show- 
ing the metal-insulator crossover between A = Eu and A = 
Gd (Yanagishima and Maeno, 2001). Bottom: Evidence for 
metal/insulator transitions with decreasing temperature for 
A = Eu, Srn and Nd but not for Pr (Matsuhira et a/., 2007). 



compounds do show a MI transition as a function of tem- 
perature. These materials were first studied in the early 
70's by Sleight and Bouchard (1972), but have not been 
studied systematically until fairly recently. Early reports 
on the electrical transport behavior were contradictory. 
One group, Lazarev and Shaplygin (1978), reported room 
temperature resistivities in the range of poor metals while 
Subramanian and Sleight (1993) argued that the entire 
serieš of compounds were low activation energy semicon- 
ductors. Also, among the earliest measurements were 
specific heat studies for A = Eľ and Lu which showed no 
magnetic anomalies up to 20 K but a rather large gamma 
coefficient which is not expected for materials reported 
to be insulating (Blacklock and White, 1980). Interest 
in these materials was re-kindled in 2001 with the report 
of magnetic anomalies above 100 K for the A = Y, Lu, 
Sm and Eu serieš members by Taira et al. (2001) and 
the study of Yanagishima and Maeno (2001) who found 
a crossover from metallic to insulating behavior with de- 
creasing A rádius and with the metal-insulator bound- 
ary between A = Eu and Gd. Representative results are 
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FIG. 28 Magnetic susceptibility for Sm2lr207 showing evi- 
dence for a magnetic anomaly near 120 K (Taira et aL, 2001). 



shown in Fig. 27. 

There exists a controversy concerning the connection 
between the magnetic and metallic properties of these 
materials. Initial reports by Yanagishima and Maeno 
(2001) indicated that the metallic serieš members, A = 
Pr, Nd, Sm and Eu were not magnetic. This is contra- 
dicted by Taira et al (2001), as can be seen in Fig. 28, 
where the metallic Sm2lr207 also clearly shows a mag- 
netic transition. 

In subsequent work, it was found that hole doping 
of Y2lr207 by substitution of Ca^+ for Y^+ induces an 
insulator to metal transition at doping levels of about 
15 atomic percent (Fukazawa and Maeno, 2002). The 
magnetic transition disappears with the onset of metallic 
behavior. One other important issue addressed here is 
the apparent large gamma coefficient of the specific heat 
which persists at low temperatures, first noted by Black- 
lock and White (1980). Fukazawa and Maeno (2002) 
found gamma to tend toward zero, 0.0(5) mJ/K^ mol Ir, 
at 0.4 K and thus concluded that Y2lr207 is indeed a 
Mott insulator. 

Nonetheless, the issue of the persistence of magnetic 
order into the metallic samples remains an important 
question. Another attempt to address this problém has 
appeared very recently by Matsuhira et al. (2007). These 
authors have investigated the correlation between sam- 
ple quality and their transport and magnetic properties. 
These studies have revealed MI transitions as the tem- 
perature is lowered for A = Eu, Sm and Nd. However, 
Pr2lr207 appears to remain metallic down to the lowest 
temperatures investigated (see Fig. 27). We return to 
Pr2lr207 in Section XI. D. For one sample, Sm2lr207, 
the onset temperatures for both the MI and magnetic 
transition were shown to be identical, suggesting that 
the two phenomena are intimately connected. The sam- 
ples showing these effects were prepared by a different 
route than those from Yanagishima and Maeno (2001) 
and Fukazawa and Maeno (2002) which involved firing 
of the starting materials in air. Matsuhira et al. (2007) 
used P t tubes in evacuated silica with a 10 percent excess 
of Ir02 and many regrinding sequences. X-ray powder 
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diífraction showed single phase samples with well resolved 
Kai/Ka2 splittings at high angles, indicative of good 
crystallinity. No evidence for changes in the diffraction 
pattern were observed below the apparent MI transitions 
from which the authors suggest that these are continu- 
ous rather than first order. As weh, the náture of the 
magnetic transition is stih unclear. The observation of a 
lambda-type anomaly in the specific heat at Tm indicated 
by the susceptibihty is evidence for a long range ordered 
antiferromagnetic ground state but more work is needed. 
For 5 = 1/2 Systems, neutrón diffraction can be a chal- 
lenge but studies on a single crystal should be definitive. 
Clearly, studies of this very interesting pyrochlore serieš 
are at an early stage and more work is needed to resolve 
the discrepancies among the various groups. 
Related 5d Pyrochlores. 

Very httle has been reported about other pyrochlore 
oxides based on 5d transition elements. The synthesis 
of 5 = Os (5ď^) and Pt (5d^) pyrochlores with A = 
Pr - Lu (including Sc and In for B = Pt) have been 
reported (Hoekstra and Gallagher, 1968; Lazarev and 
Shaplygin, 1978). The B = Ft serieš can only be pre- 
pared using high pressures, 4 GPa, and high tempera- 
tures, 1200°C (Hoekstra and Gallagher, 1968). Room 
temperature electrical resistivity values are reported for 
the B = Os serieš and all are in the range for poor met- 
als, very similar to the initial reports for the B = li 
compounds described above. Cd20s207 was first stud- 
ied by Sleight et al (1974) and later by others for its 
unusual MI transition at 226 K (Mandrus et a/., 2001). 
This suggests that a closer study of this serieš could be 
very interesting. 



G. ^2Mo207 (A = Gd, Nd and Srn) 

The earliest report of the existence of rare earth molyb- 
denum pyrochlores appears to be that of McCarthy 
(1971) from a study of the Eu - Mo - O and Sm - Mo 
- O phase diagrams in which the pyrochlore compounds 
were observed. Hubert (1974) was the first to report 
magnetic susceptibihty for Y2M02O7 in the range 300 to 
1000 K. Ranganathan et al. (1983) synthesized the serieš 
A = Gd - Lu along with solid solutions in which the A- 
site contained various ratios of Nd/Yb and Nd/Er and 
reported magnetic susceptibihty and limited electrical 
transport data for 77 K - 300 K. The observation of pos- 
itive Curie-Weiss {9cw) temperatures led the authors to 
suggest that some of these materials might be ferromag- 
nets, for example Sm2Mo207 (Mandiram and Gopalakr- 
ishnan, 1980). As well, samples rich in Nd appeared to 
be metallic (Ranganathan et al, 1983). This situation 
was clarified by Greedan et al (1987) who showed that 
Mo(IV) pyrochlores for A = Nd, Sm and Gd were indeed 
ferromagnets with Tc = 97, 93 and 83 K, respectively, 
which was attributed to the ordering of the Mo(IV) mo- 
ments, an unprecedented observation. As mentioned in 
Section VII, there is a link between the ionic rádius of 



the A-ion and the electrical transport properties of the 
molybdenum pyrochlores which in turn are strongly cor- 
related to the magnetic properties. 

The existing evidence appears to support the view that 
the properties of ferromagnetism and metallic behavior 
are fundament ally linked in these materials. A proposal 
for the origin of the ferromagnetism (and antiferromag- 
netism for Y2M02O7) has been advanced by Solovyev 
(2003). While the argument is complex and has a num- 
ber of elements, the key idea is that the Mo t2^ states are 
split by the axial crystal field component into ai^ and e^ 
states. These states are effected differently by systematic 
changes in crystal structure in proceeding from A = Nd 
to Y. For A = Nd and Gd, the e^ band is found to be rela- 
tively broad and can support intinerant spin up electrons, 
rather in analogy to the double exchange mechanism in 
the manganate perovskites which selects a ferromagnetic 
ground state. However, for Y2M02O7, the e^ states be- 
come more localized, while the ai^ band broadens. This 
is found to favor an antiferromagnetic ground state. The 
marked differences in magnetism and transport proper- 
ties are very surprising given that the structural changes 
in moving from A = Nd Gd^ Y are very subtle. 
For example the Mo - O - Mo bond angle changes from 
131.4° [Nd] to 130.4° [Gd] to 127.0° [Y] (Moritomo et al, 
2001; Reimers et al, 1988). 



1. Gd2Mo207 

Gd2Mo207 has attracted much interest due to its posi- 
tion near the metal/insulator (MI) boundary (see Fig. 7) 
in this molybdate serieš. As already mentioned, the elec- 
trical transport properties are dependent on the stoi- 
chiometry of the sample (see Fig. 8). Stoichiometric sam- 
ples are ferromagnetic and metallic with Tc near 80 K. 
These show a giant negatíve magneto-resistance exceed- 
ing 40 percent below 15 K (Troyanchuk et al.^ 1998). In- 
sulating samples are oxygen deficient and electron doped. 
It has been shown by Hanasaki et al. (2006) that insu- 
lating samples of Gd2Mo207 can be driven metallic by 
application of high pressure with the MI boundary occur- 
ring between 3 and 4 GPa (see Section XII. B). Extensive 
specific heat studies have been carried out on polycrys- 
talline samples for both A = Sm and Gd by Schnelle and 
Kremer (2004). Unlike previous studies by Raju et al. 
(1992), clear maxima were observed for both materials 
near 75 K, indicative of long range magnetic order on 
the Mo sublattice. For Gd2Mo207, a step-like anomaly 
at 11.3 K superimposed on a dominánt Schottky peak 
suggests partial ordering of the Gd spins, which was sup- 
ported by measurement of the entropy. 

Until quite recently, there has been little Information 
regarding the magnetic structure of either the A = Sm 
or Gd molybdate pyrochlores, due to the high neutrón 
absorption cross sections for both elements. Mirebeau 
et al. (2006a) have solved this problém with the use of 
^^^Gd substituted materials. The room temperature unit 
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FIG. 29 Neutrón powder diffraction results for Gd2Mo207 
by Mirebeau et al. (2006a). Top: The magnetic scattering at 
1.7 K after a 90 K data set has been subtracted. The sohd 
line is a fit to a model with colhnear ferromagnetic couphng 
between Gd and Mo moments. Bottom: The temperature 
dependence of the ordered moments for the Gd and Mo sub- 
lattices. The inset shows the temperature development of the 
correlation lengths (L) obtained from the widths of the Bragg 
peaks and the Bragg R for the magnetic structure. 



celí constant for this sample is 10.3481(2) Ä, which in- 
dicates only a very slight oxygen doping. The magnetic 
diffraction pattern at 1.7 K is shown in Fig. 29 along with 
the temperature dependence of the sublattice moments 
and correlation lengths. The absence of a (200) refiec- 
tion is consistent with a collinear magnetic structure and 
the best fit occurs for a ferromagnetic coupling of the Gd 
and Mo sublattices. This agrees with the earliest mag- 
netization studies of Sato et al. (1986) who found that 
the bulk saturation moment could only be understood in 
terms of a ferromagnetic Gd - Mo coupling. 

From the right panel of Fig. 29 it is clear that the 
Mo sublattice orders above 80 K and that the Gd mo- 
ments are polarized due to the Gd-Mo coupling. The 
ordered moments on both sublattices are significantly 
smaller than the spin only values of 2.0 jiB and 7.0 jiB for 
Mo^+ and Gd^+ , respectively, although the Mo moments 
are similar to those found for the A = Nd phase. Muon 
spin relaxation data (Apetrei et al, 2007), which are 
most sensitive to the Gd spins, indicate that strong spin 
fluctuations persist below Tc and as low as 6.6 K. This 
is consistent with the ^^^Gd Môssbauer data of Hodges 
et al (2003) and also, perhaps, with the specific heat 
data of Schnelle and Kremer (2004). The conclusion is 
that Gd2Mo207 is perhaps an unconventional ferromag- 
net with strong spin fluctuations. This maybe analogous 



(a) Mo (c) Mo 




(b) Nd ^ (d) Nd 



3. r 




FIG. 30 Proposed spin configurations [(a), (b) and (c), (d)] for 
the Mo and Nd moments from neutrón diffraction on a single 
crystal of Nd2Mo207 at 4 K (Yasui et al, 2001). 



to the situation encountered for manganate pyrochlores 
in Section VIITE. The remarkable properties of this ma- 
teriál under applied pressure will be described in Section 

xn.B.2 



2. Nd2Mo207 

Magnetism. It is well established that the Nd2Mo207 
phase is a metallic ferromagnet with Tc = 95 K for sto- 
ichiometric, polycrystalline samples. That the ferromag- 
netism is due to the Mo sublattice was first demonstrated 
by Greedan et al (1991) from powder neutrón diffrac- 
tion studies. A canted ferromagnetic structure on the 
Mo sites with an ordered moment of 1.1 ja b was pro- 
posed from analysis of data at 53 K which are dominated 
by magnetic scattering from the Mo. Subsequent stud- 
ies by Yasui et al. (2001) using a single crystal with Tc 
= 93 K have provided a better defined magnetic struc- 
ture for both the Nd and Mo sublattices. The magnetic 
structure at 4 K is shown in Fig. 30. The Mo spins are 
non-collinear making an angle of 9° with the local z axis 
which is parallel to a four fold axis of the crystal. The 
Mo moment is 1.2 /íb, very close to that found earlier 
and slightly greater than half the value expected for a 
spin only S = 1 ion. The Nd moments are aligned nearly 
along the principál 3-fold rotation axes of the tetrahedra, 
in a 2-in/2-out pattern. Nd moments at 4 K are 1.5 /hb 
and the relative orientations of the Nd and Mo sublattice 
moment directions are antiferromagnetic. 
Anomalous Hall Effect. Much attention has been fo- 
cussed on Nd2Mo207 since the observation of the un- 
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FIG. 31 (a) Temperature dependence of the bull: mag- 
netic moment and the (200) magnetic reflection which traclís 
the development of the chiral order on the Nd sites in 
Nd2Mo2 07. (b) Temperature dependence of Hall conductiv- 
ities for Nd2Mo2 07 and Gd2Mo2 07. Note the conventional 
behavior for A — Gd while the conductivity remains large and 
hnite at low temperature for A — Nd (Taguchi et a/., 2004). 



precedented behavior of the so-called anomalous Hall ef- 
fect (AHE) independently by Katsufuji et al (2000) and 
Taguchi et al. (2001). In ferromagnets, the transverse or 
Hall resistivity has two contributions, the "ordinary" Hall 
coefficient, Rq^ which is proportional to the applied mag- 
netic field B and the "anomalous" coefficient, Rs which 
is proportional to the sample magnetization, M, as in 
Eq. (6). 



pH = RoB + AttRsM. 



(6) 



The usual behavior is for the AHE contribution to van- 
ish at low temperatures in ferromagnetic metals with 
collinear spin configurations. However, in the case of 
Nd2Mo207, the AHE actually increases and remains 
large and finite at T = O (see Fig. 31). In one inter- 
pretation, this behavior has been attributed to the spin 
chirality (the two-in, two-out state is 6-fold degenerate) 
of the Nd^+ moment configuration which acts as an ef- 
fective magnetic field and effects the carrier dynamics in 
the same way as a real magnetic field (Taguchi et a/., 
2001, 2004). A strong point of this argument is the con- 
trast between this AHE and that in the Heisenberg, A 
= Gd phase where spin chirality is not present. For the 
A = Gd materiál the AHE vanishes at low temperatures 
as with a normál ferromagnetic metal. However, this in- 
terpretation has been questioned by Yasui et al. (2003a) 
and Kageyama et al. (2001) who argue, based on exten- 
sive neutrón scattering and specific heat studies, that the 



origins of the AHE for Nd2Mo207 are mu ch more com- 
plex than thought originally and that the spin chirality 
mechanism alone cannot provide a quantitative explana- 
tion and may at best play only a minor role (Kageyama 
et al, 2001; Yasui et al, 2003, 2006). On the other hand, 
Kézsmárki et al. (2005) interpret magneto-optical data 
on both Nd2Mo207 and Gd2Mo207 in favor of an im- 
portant contribution from spin chirality. At present this 
controversy appears to be unresolved. Further comment 
on this problém is presented in Section VIII. G. 4 below. 



3. Sm2Mo207 

Stoichiometric samples of this materiál generally show 
Tc = 93 K, although as already mentioned in Section VII, 
for oxygen deficient samples this value can be signifi- 
cantly reduced (see Fig. 9). A crystal of Sm2Mo207 
with a Tc = 73 K was reported by Taguchi and Tokura 
(1999) to show also a finite anomalous Hall Effect coef- 
ficient at T = O K. A giant negatíve magnetoresistance 
of 13 percent near Tc and 18 percent at 4 K has been 
reported (Taguchi and Tokura, 2000), but little more is 
known about Sm2Mo207. 



4. A2{W\oB)207 

Studies of B-site substituted compounds are very 
sparse. Troyanchuk et al. (1998) have reported the prepa- 
ration of V - doped Gd2Mo207 with B-site composition 
M01.2V0.8- Gd2Mo207 is a ferromagnetic metal while 
the corresponding Gd - V pyrochlore does not exist even 
under high pressures. Nonetheless, the A2V2O7 serieš 
members, which include only A = Lu, Yb and Tm at am- 
bient pressure are all ferromagnetic semiconductors with 
Tc in the range 70 K - 75 K (Bazuev et al, 1976; Shin-ike 
et al, 1977; Soderholm et al, 1980). At this particular 
composition, the solid solution is insulating with resistiv- 
ities in the range of íl-cm but remains ferromagnetic with 
a slightly enhanced Tc relative to the undoped materiál 
and a GMR effect of about 40 percent seen for the pure 
materiál has been destroyed (Troyanchuk et al, 1998). 

In another case, the AHE has been studied for 
Nd2Mo207 with Ti substitutions for Mo up to the B- 
site composition M01.7TÍ0.3 (Kageyama et al, 2001). At 
this doping level the compound is still metallic. This 
group had shown earlier that the AHE has two contri- 
butions, proportional to the magnetizations of Mo and 
Nd separately, that is, (likubo et al, 2001; Yoshii et al, 
2000) 



PH = RoB + AnRsMiMo) + AnR'^MCNá) (7) 

Of these two components, the doping experiments 
showed that Rg changes sign from positive to negatíve 
with increasing Ti content while R'^ does not. This ob- 
servation was taken as evidence that the mechanism for 
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AHE in these materials is more complex than the pure 
chirality model advanced initially (Taguchi et al.^ 2001). 



IX. SPlN GLASS PHASES 

The spin glass state is one where the combination of 
(i) randomness and (ii) frustration prevents the develop- 
ment of conventional long range magnetic order charac- 
terized by delta-function magnetic Bragg peaks (Binder 
and Young, 1986). There are a number of experimental 
pathways to the spin glass state. The classic case involves 
dilute concentrations of magnetic atoms in a metallic, 
diamagnetic host such as Aui-a^Fca^ where x ^ 0.05, al- 
ready mentioned in Section I. Here the effective coupling 
between Fe moments is mediated by the Au conduction 
electrons, giving rise to an RKKY interaction whose sign 
depends on the dištance between two Fe moments. This 
leads to competing random ferro- and antiferromagnetic 
constraints at each magnetic site. Another approach is to 
introduce disorder into a magnetically ordered but geo- 
metrically frustrated materiál by dilution of the magnetic 
sites with diamagnetic ions to levels near the percola- 
tion threshold for the particular lattice. An example is 
Eui-ícSra^S where Sr^+ dilution of the face centered cubic 
Eu^+ sites destroys the ferromagnetic ground state and a 
spin glass state emerges near x = 0.5. Similar to the lat- 
ter example is LiRoxYi-x^4^ where Ho^+ is an effective 
Ising spin and where the predominant interactions are 
dipolar. Dipolar interactions are intrinsically frustrated 
since they can be either ferromagnetic or antiferromag- 
netic depending on the orientation of the vector r j — Vi 
between spins at positions r i and r j. As a result, substi- 
tution of Ho by Y generates random frustration. Indeed, 
LiHo^^Yi 

-x^A is an Ising ferromagnet for x down to about 
x ~ 0.2. For x < 0.2, a dipolar spin glass state devel- 
ops (Reich et al.^ 1990; Wu et al, 1993), although this 
conclusion has recently been questioned (Jônsson et al.^ 
2007). 

One question that has attracted much attention in the 
field of spin glasses is whether or not a thermodynamic 
spin freezing transition occurs at nonzero temperature 
for physical dimensions. While it was determined some 
time ago that there is no spin glass transition in two 
dimensions, a seemingly definite conclusion for three- 
dimensional (3D) systems has only recently become avail- 
able, thanks to the work of Ballesteros et al. (2000) who 
pioneered the use of the scaled correlation length method 
in finite-size scaling analysis of spin glass models. A 
large majority of real magnetic systems are better de- 
scribed by Heisenberg spins rather than Ising spins and 
it was believed that the 3D XY and Heisenberg spin glass 
models do not have a transition. The explanation for 
the experimentally observed transition in real systems 
therefore had to invoke random anisotropy as the mech- 
anism responsible for driving the systém into the Ising 
spin glass universality class (Binder and Young, 1986). 
However, recent extensive Monte Carlo simulations on 
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FIG. 32 Left: D C magnetic susceptibility for Y2M02O7 at an 
applied field of 100 Oe (Gingras et al, 1997). Right: Heat 
capacity (a) and entropy (b) for Y2M02O7. Note the broad 
maximum below 20 K and the linear dependence on temper- 
ature below 7 K (Raju et ah, 1992). 



the 3D XY and Heisenberg models employing the Balles- 
teros et al. (2000) method find some compelling evidence 
for a phase transition at nonzero temperature in these 
two systems (Lee and Young, 2003), although the lower- 
critical dimension for these models for a nonzero tran- 
sition temperature seems to be very close to three (Lee 
and Young, 2007). 

From a theoretical perspective, models with couplings 
of random signs, such as the so-called Edwards-Anderson 
model, have attracted by far the most attention. One 
generally believes that as long as there is a transition at 
nonzero temperature, the universality class should be the 
same irrespective of the details of the model, i. e., either 
continuous or discrete distributions of random bonds or 
a randomly diluted frustration. 

Experimentally, a number of signatures are taken as 
indicative of the spin glass state, such as the frequency 
dependence of x' the AC susceptibility, the linear tem- 
perature dependence of the low temperature heat capac- 
ity or the /iSR line shape, but the definitíve approach is 
the measurement of the temperature dependence of the 
nonlinear magnetic susceptibility, xs (Binder and Young, 
1986). The DC magnetization, Mz{T^Bz)j can be ex- 
panded as a Taylor serieš of the applied magnetic field 
B z as 

M,(T, Bz) ^ Xi{T)Bz - Xs{T)Bz^ , (8) 

where T is the temperature. At a second order spin 
glass transition, one expects X3(^) show a power-law 
divergence as Xsi^) ~ {T — Tf)~^ , with 7 a critical expo- 
nent characterizing the spin glass transition at the freez- 
ing temperature Tf (Binder and Young, 1986; Gingras 
et al.^ 1997). Other critical exponents can also be de- 
termined by measuring the full nonlinear magnetic field 
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dependence of Mz{Bz^T) (Binder and Young, 1986; Gin- 
gras et al.^ 1997). 

With this very minimal background materiál in hand, 
we can now discuss the spin glass behaviors observed in 
some of the pyrochlore oxides. 



A. Y2M02O7 and Tb2Mo207 

1. Y2M02O7 

The initial report by Greedan et al. (1986), showing a 
canonical spin glass behavior, (see Fig. 32a) for Y2M02O7 
has sparked considerable interest which continues to the 
present time. Neutrón powder diffraction data, albeit 
of moderate resolution, could be analyzed in terms of 
a fully ordered pyrochlore model (Reimers et al.^ 1988). 
Specific heat studies by Raju et al (1992) showed only 
a broad maximum near the apparent Tj = 22 K and 
a linear dependence on temperature at low T, another 
feature typical of spin glasses (see Fig. 32b). AC suscep- 
tibility (Miyoshi et al.^ 2000), thermo-remanent magne- 
tization by Dupuis et al (2002) and Ladieu et al. (2004) 
provide more evidence for the canonical spin glass behav- 
ior of this materiál. The AC susceptibility in nearly zero 
field shows the classic frequency dependence (see Fig. 33). 

The best evidence for the nearly canonical spin glass 
character came from measurements of the non-linear sus- 
ceptibility, by Gingras et al (1996, 1997). The non- 
linear magnetization could be analyzed, (see Fig. 34) ac- 
cording to a scaling model for phase transitions, yielding 
critical exponents 7, p and ô which satisfied the scaling 
relationship 



(9) 



with values T/ = 22 K, 7 = 2.9(5), (3 = 0.8(2) and ó ^ 
4.7. These agree well with those found for conventional 
spin glasses with dilute magnetic centers and positional 
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FIG. 33 Left: Neutrón spin echo results for Y2M02O7 at 
temperatures spanning T/ = 22.5 K determined from státie 
magnetization data. Note that within this time window, 
spin freezing is not fully established until 10 K (Gardner 
et al, 2001, 2004). Right: Frequency dependent AC suscep- 
tibility, x^ for Y2M02O7 showing classical spin glass behav- 
ior (Miyoshi et al, 2000). 




FIG. 34 Nonlinear magnetization analyzed according to a 
scaling model for a T/ = 22 K and 7 = 2.8 and (3 = 0.75. 
Inset shows a log-log plot which allows determination of the 
critical parameters with ô = 4.73 which follows from the scal- 
ing law with the above 7 and (3 values (Gingras et a/., 1997). 



disorder (Binder and Young, 1986), making Y2M02O7 
indistinguishable from such systems. 

Strong evidence for spin freezing comes also from stud- 
ies of spin dynamics by muon spin relaxation (Dunsiger 
et al.^ 1996). Fig. 35 shows 1/Ti, the muon spin relax- 
ation rate, as a function of temperature. These data show 
features typical of disordered spin-frozen systems such as 
a critical slowing down as Tf is approached from above 
followed by a sharp decrease. A finite 1/Ti persists down 
to the lowest temperatures studied, but this is at the res- 
olution limit of this technique and more needs to be done 
to conclusively determine if the spins are still dynamic in 
the mK temperature range. 

Neutrón scattering experiments, Gardner et al 
(1999a), provide more insight into the possible local or- 
dering and the spin dynamics over a wide temperature 
range. Figúre 36 shows the |Q| dependence of the mag- 
netic elastic scattering at 1.8 K. The obvious feature is 
a broad peak centered at |Q| = 0.44 Ä~^. In terms 
of unit celí dimensions, this value of |Q| corresponds to 
27r/d(110), which involves the cubic face diagonál. From 
the half-width at half maximum (HWHM) a correlation 
length can be estimated as (f = 1/HWHM = 5 Ä. These 
two observations are consistent with a four sublattice 
structure as depicted in the inset of Fig. 36. 

Inelastic scattering data give a detailed picture of the 
evolution of the spin dynamics in terms of four distinct 
regimes: (1) For T > 200 K there are no discernable 
correlations in either space or time. (2) Between 200 K 
and T = 2Tf spatial correlations build up, peaking at 
IQI = 0.44 Ä-\ but the energy distribution is broad. 
(3) Within the interval T/ < T < 2T/, the spatial cor- 
relations are no longer evolving but the spin fluctuation 
rate decreases proportionally to {T — Tf) in a manner 
appropriate to the approach to a phase transition. (4) 
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FIG. 35 The muon spin relaxation rate, 1/Ti versus tempera- 
ture. Note the critical slowing down near T f and the presence 
of a small but finite relaxation rate persisting to very low tem- 
peratures (Dunsiger et al, 1996). 



Below Tf the spin fluctuation rate is small and changes 
very little, consistent with a highly spin frozen state. 

The neutrón spin echo technique provides Information 
on spin dynamics over a shorter time scale than /iSR, 
from about 1 to 10~^ ns. Within this time window, a 
diíferent view of the spin dynamics emerges (Gardner 
et al, 2001, 2004). The left panel of Fig. 33 shows the 
time dependence of the normalized intermediate scatter- 
ing function, S{Q^t)/ S{Q^O) for temperatures spanning 
T f (22.5 K) determined from státie magnetization stud- 
ies. While the data for 10 K show spin- frozen behavior, 
there are still significant spin dynamics below Tf. Fits 
to these data indicate that spin relaxation times increase 
by a factor of 10^ between 25 K and 10 K, from 0.07 ns 
to 8 x 10^ ns. 

The evidence presented above raises a fundamental 
question: why is Y2M02O7 such a typical spin glass? In 
generál, insulating antiferromagnetic spin glasses involve 
dilution of the magnetic centers by diamagnetic ions at 
a concentration below the percolation limit which intro- 
duces both positional disorder and frustration simultane- 
ously, two conditions considered necessary for the estab- 
lishment of the spin glass ground state. In Y2M02O7 the 
magnetic frustration is provided by the magnetic lattice 
topology but is there also disorder of some subtle type? 
While powder neutrón diffraction data show that the av- 
erage structure is well described by the fully ordered py- 



rochlore model (Lozano et a/., 2007; Reimers et al.^ 1988), 
local structure probes such as EXAFS, and NMR have 
suggested the presence of disorder at some level in this 
materiál. Mo K edge EXAFS data were interpreted by 
Booth et al. (2000) to show that the variance in the Mo - 
Mo nearest-neighbor dištance was about ten times larger 
that that for the Mo - O and Mo - Y distances. From 
this and a number of assumptions it was estimated that 
the level of disorder introduced in a pair-wise exchange 
constant would be of the order of 5 percent. This was 
judged to be about a factor of íive too small when com- 
pared with predictions from the Sherrington and South- 
ern (1975) mean field theory for spin glasses. Keren 
and Gardner (2001) carried out ^^Y NMR experiments 
over a limited temperature range. Data obtained above 
200 K (the Curie-Weiss temperature) showed a smooth 
and broad resonance but at 200 K and 92.4 K, a large 
number of peaks of small amplitúde become superim- 
posed on the broad feature (see Fig. 37). This was inter- 
preted in terms of a distribution of Y environments due 
to localized "lattice" distortions and led to speculation 
that these distortions are frustration driven. A similar 
experiment was performed using muon spin relaxation 
by Sagi et al. (2005) down to 20 K. The width of the 
internal field distribution as detected by the muon grows 
upon cooling at a rate which cannot be explained by the 
increasing susceptibility alone. Therefore, they conclude 
that the width of the distribution of coupling constants 
also grows upon cooling, or in other words, the lattice 
and spin degrees of freedom are involved in the freezing 
process in Y2M02O7. 

A model of Heisenberg spins with magnetoelastic cou- 
pling was invoked to account for the high freezing temper- 
ature. The derivative of the exchange strength relative 
to bond length was found to be 0.01 eV/Ä and the elas- 
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FIG. 36 Elastic neutrón scattering for Y2M02O7 shown as 
the difference between data at 1.8 K and 50 K. The broad 
peak suggests the short range order model shown in the in- 
set (Gardner et a/., 1999a). 
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FIG. 37 Y-NMR data for Y2M02O7 showing multiple peaks 
arising below 200K, indicating the presence of multiple Y sites 
(Keren and Gardner, 2001) 
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tie constant to be 0.1 eV/Ä^. In ZnCr2 04, where mag- 
netoelastic coupling appears to be import ant, these val- 
ues are 0.04 eV/Ä and 6.5 eV/Ä^ respectively (Sushkov 
et al.^ 2006). Other theoretical studies that have con- 
sidered classical pyrochlore Heisenberg antiferromagnets 
with random variations of the exchange couplings in- 
clude Bellier-Castella et al. (2001) and Saunders and 
Chalker (2007). 

At the time of writing, the microscopic mechanism be- 
hind the spin glass transition in Y2M02O7 is unresolved. 
Ideally, one would like to find a means to reconcile the 
diífraction studies with the local próbe results. One 
possibility is the application of the neutrón pair distri- 
bution function method (Billinge, 2004; Proffen et al.^ 
2003). Studies using this approach are currently on- 
going (Lozano et al.^ 2007). 

Dunsiger et al. (1996a) showed that a 20% dilution of 
the magnetic Mo^+ site with nonmagnetic Ti^+ reduces 
the freezing temperature but increases the residual muon 
relaxation rate, indicating an increased density of states 
for magnetic excitations near zero energy. While there 
are no other published reports on 5-site substitution in 
Y2M02O7, there may be much to learn from such efforts. 
As already noted, the microscopic origin of the spin glass 
behavior in this materiál is not understood. What is 
known suggests that Y2M02O7 falls into the category of 
a bond disordered spin glass rather than a site disordered 
glass such as Eui-a^Sr^^S, for example. In the context of 
a picture in which Y2M02O7 is an "intrinsic" random 
bond spin glass, it would be of interest to track Tf as a 
function of the doping of the Mo site with a non-magnetic 
ion such as Ti^+. Interesting physics may emerge at low 
temperature close to the percolation limit which is 39 
percent of magnetic site occupancy for the pyrochlore 
lattice. 



FIG. 38 Diffuse magnetic elastic scattering for Tb2Mo207 at 
various temperatures (Gaulin et al, 1992). 



2. Tb2Mo207 

Although this materiál has received much less atten- 
tion than the Y-based pyrochlore, there are strong sim- 
ilarities. Tb2Mo207 is an inherently more complex sys- 
tém with two geometrically frustrated magnetic sublat- 
tices. Also, the A = Tb phase is just on the insulating 
side of the MI boundary for this serieš and is not a bulk 
ferromagnet although the Curie-Weiss temperature is re- 
ported to be 6cw = +17K (Sato et al.^ 1986). Among the 
few detailed studies of transport properties, the materiál 
is shown to exhibit an unusual magnetoresistance with 
both positive (10 percent) and negatíve (30 percent) MR 
at low and high fields, respectively (Troyanchuk et al, 
1988). This behavior lacks an explanation to date. 

Magnetization and neutrón scattering studies indi- 
cated spin glass like behavior with Tf = 25K and intense 
diffuse magnetic scattering (Gaulin et a/., 1992; Greedan 
et al.^ 1990). The diffuse magnetic scattering sets in be- 
low 50 K (see Fig. 38), showing two very broad peaks 
near |Q|= 1.0 and 2.0 Ä~^. This pattern is quite dif- 
ferent from that for Y2M02O7 (see Section IX. A. 1) and 
indicates that the Tb sublattice scattering is dominánt 
here. Note also the upturn at low Q. 

SANS data were also reported over the Q-range 0.019 
Ä-^ to 0.140 Ä-^ Subtraction of the 300 K data from 
low temperature data shows that non-zero SANS ap- 
peared only above |Q| = 0.14 Ä~^, so it is unlikely that 
this SANS tail is of ferromagnetic origin (Greedan et al.^ 
1990). However, these results might be in conflict with 
a recent report from Apetrei et al. (2007) who observe 
SANS above |Q| = 0.25 Ä-^ This issue should be clar- 
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FIG. 39 The muon spin relaxation rate vs temperature for 
Tb2Mo2 07 in an applied field of 5 mT. The inset shows the 
Tb moment fluctuation rate above the spin freezing point 
compared with neutrón data indicated by the inverted tri- 
angles (Dunsiger et a/., 1996). 
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iíied. From inelastic neutrón scattering, the Tb spins 
have been shown to be fluctuating at about 0.02 THz 
above Tj, but spin freezing within the experimental time 
window was seen below 25 K (GauHn et al.^ 1992). Spin 
relaxation studies were extended into the /iSR time win- 
dow (see Fig. 39). Comparing with Fig. 35, which shows 
corresponding data for Y2M02O7, one notes some simi- 
larities and differences. First, the data above T f are ver y 
similar, indicating a critical slowing down behavior, fol- 
lowed by a clear maximum and a subsequent decrease. 
Note, however, that the relaxation time below about 1 K 
remains very large (5 /is), relative to that for A = Y (0.02 
/is). This difference is attributed to the larger moment 
of Tb^+ relative to Mo^+ being roughly proportional to 
the ratio of the squares of the moments. These results 
show that the Tb spins are involved in the freezing and 
that there exists an appreciable density of states for low 
energy magnetic excitations for these materials which are 
accessible even at very low temperatures. 



3. Other A2M02O7 

In fact there is relatively little published Information 
on other members of this serieš for A = Dy, Ho, Er, Tm, 
Yb or Lu. Most show DC and AC susceptibility anoma- 
lies near T = 22 K, as in the A = Y materiál (Raju and 
Gougeon, 1995). Miyoshi et al (2000) show AC suscep- 
tibility data for A = Ho and Tb. The A = Yb materiál 
has been studied using ^^^Yb Môssbauer measurements. 
The major results are that evidence was found for a lower 
than axial symmetry at the Yb site, suggesting some lo- 
cal disorder and that the local magnetic fields acting on 
the Yb nucleus show a random distribution, consistent 
with the bulk spin glass properties (Hodges et aL, 2003). 

There have been a few studies involving mixed occu- 



FIG. 40 A phase diagram for {Ai-xAx)Mo207 pyrochlores. 
(a) Tc or Tf. (b) Resistivity at 5K. (c) Magneto-resistance at 
5K and Tc (d) Unit celí constant (Katsufuji et a/., 2000). 



pation of the A site. Katsufuji et al (2000) investigated 
materials of the type Ai-xA'^ with A = Nd, Sm and 
Gd and A' = Tb and Dy. That is, solid solutions be- 
tween the ferromagnetic phases, A = Nd, Sm and Gd 
and the spin glass materials A^ = Tb and Dy. The re- 
sults are summarized in Fig. 40 where the FM metal to 
spin glass insulator transition occurs at an average A-site 
rádius between that for Gd^+ [1.053 Ä] and Tb^+ [1.04 
Ä] . This is fully consistent with the results from the pure 
samples. Two studies have been reported in which the 
A-site is substituted with non magnetic ions and the re- 
sults are markedly different. In some very early work by 
Sato and Greedan (1987), solid solutions of the type A 
= Yi-x'^dix were investigated. Compositions were chosen 
to span the average A site rádius from Nd to Y, intending 
to cover the ferromagnet to spin glass transition range. 
Surprisingly, perhaps, all samples showed only spin glass 
behavior with Tf values not far from that for A = Y, 
Tf ^ 22 K. Even the compositions A = Yo.sLao.s and 
Yo.6Lao.4, which have an average A site rádius equiva- 
lent to Nd and Sm, respectively, showed no ferromag- 
netic transition, although the Curie temperatures were 
positive, +41 K and +31 K respectively. Unfortunately, 
no electrical transport data have been reported. 

Hanasaki et al. (2007) chose solid solutions with A = 
Eu (Eu^+is, technically, a non-magnetic ion as J = 0) 
and A' = Y or La. The pure A = Eu is reported 
to be metallic and ferromagnetic with T = 50 K or 
so (Kézsmárki et al.^ 2006). The data are from an im- 
age furnace grown crystal and unfortunately no unit celí 
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constant is reported, so some electron doping is likely. 
Hanasaki et al. (2007) report ferromagnetic behavior for 
the La-doped materials but with a re-entrant spin glass 
transition at 22 K. Tc for the Euo.85Lao.15 materiál ap- 
pears to be about 65 K. For this A-site composition the 
average A-site rádius is equivalent to that of Sm^+ and 
one might expect a higher value. Perhaps, electron dop- 
ing is playing a role again. The phases with A = (Eu, 
Y) are spin glasses. These results are quite new and in- 
teresting but the systém should be investigated further. 
Finally, the discrepancy between the results of Sato and 
Greedan (1987) where the A-site involves La and Y and 
those of Hanasaki et al. (2007), just described, is difRcult 
to understand. One obvious difference is that the vari- 
ance of the average rádius is of course greater for the A 
= La, Y combination than for Eu, La or Eu, Y. Thus, 
the influence of the A-site composition on the magnetic 
and transport properties is apparently subtle and merits 
closer investigation. 



X. SPIN ICE PHASES 

The spin ice phenomenology in the pyrochlore oxides 
was discovered in 1997 in H02TÍ2O7 (Harris et al.^ 1997, 
1998a). This systém possesses a ferromagnetic Curie- 
Weiss temperature 6cw ~ +2 K. It was thus quite sur- 
prising that this materiál does not develop long range 
order down to 50 mK, with neutrón scattering revealing 
only broad diffuse scattering features for temperatures 
mu ch lower than ^cw- 

In H02TÍ2O7 the strong axial crystal electric field act- 
ing on Ho^+ gives rise to an almost ideál, classical Ising 
spin. Because of symmetry, the local Ising (quantiza- 
tion) axis is along the local cubic (111) directions, such 
that on a tetrahedron, a spin can only point "in", to- 
wards the middle of the opposing triangular face, or op- 
positely, hence "out" of the tetrahedron. While geomet- 
rie frustration is generally associated with antiferromag- 
netic interactions, here frustration arises for ferromag- 
netic ones. Indeeed, there are six "two in/two out" spin 
configurations that minimize the ferromagnetic exchange 
energy on an individual tetrahedron, and thus six ground 
States. There is an infinity of such ground states for a 
macroscopically large sample, and there is therefore an 
extensive ground state entropy. This entropy, 5*0, can be 
estimated by borrowing the argument used by Pauling 
for estimating the residual protón configuration disor- 
der in common hexagonal water ice (Pauling, 1935). The 
main point is to consider the difference between the num- 
ber of constraints necessary to determine a ground state 
and the number of degrees of freedom that the systém 
possesses. Consider Anderson's Ising pyrochlore antifer- 
romagnet, onto which the local (111) pyrochlore Ising 
model maps, as discussed above (Anderson, 1956). The 
ground state condition is "underconstrained" , demand- 
ing only that the total magnetization of the four Ising 
spins on each tetrahedron be zero. Six of the 2^ = 16 pos- 



sible spin configurations satisfy this condition. Counting 
2^ configurations for each tetrahedron gives, for a sys- 
tém of N spins and N/2 tetrahedra, a total number of 
microstates, Q = (2^^)^/^ = 4^. This number drastically 
overestimates the exact total, Q = 2^ . The reason is 
that each spin is shared between two tetrahedra, hence 
the above 16 configurations on each tetrahedron are not 
independent. Following Pauling's argument, we allocate 
2^ = 4 states per tetrahedron and, assuming that 6/16 of 
them satisfy the constraint, this leads to a ground state 
degeneracy í^o = {22(6/16)}^/^ = (3/2)^/^. The corre- 
sponding entropy 5*0 = /cb log(í]o) = (^^b)/2 log(3/2) is 
of course just Pauling's originál result. 

Not only is the residual entropy of ferromagnetic (111) 
spins on the pyrochlore lattice the same as Pauling's en- 
tropy for water ice but, as shown in Fig. 41, there is also 
a rather direct connection between the spin configura- 
tions in the pyrochlore problém and that of the protón 
positions in water ice. For this reason, the term spin 
ice was coined. In anticipation of the forthcoming dis- 
cussion of the physics at play in Tb2TÍ207, it is worth- 
while to comment on a case where the nearest-neighbor 
exchange interactions are antiferromagnetic for a situ- 
ation with local (111) spins. In that case, the ground 
state consists of all spins pointing in or out of a refer- 
ence tetrahedron. Hence, there are in that case only two 
ground states related by a global spin inversion symme- 
try, and a second order transition in the universality class 
of the (unfrustrated) three-dimensional Ising model is ex- 
pected. Earlier, Anderson had noticed the connection be- 
tween the statistical mechanics of antiferromagnetically 
coupled Ising spins on the pyrochlore lattice and Paul- 
ing's model of protón disorder in water ice. However, 
in Anderson's model, the Ising spins share a common 
(global) z— axis direction, and frustration arises as usual 
for antiferromagnetic interactions with spins on triangu- 




FIG. 41 Illustration of the equivalence of the water ice rule 
"two protons near, two protons far" and the spin ice rule 
"two spins in, two spins out". The diagram (left) illustrates 
a water molecule in the tetrahedral coordination of the ice 
structure with protón positions located by displacement vec- 
tors that occupy a lattice of linked tetrahedra. In spin ice 
(right) the displacement vectors are replaced by rare earth 
moments ("spins") occupying the pyrochlore lattice, which is 
the dual lattice (i. e. the lattice formed by the mid-points of 
the bonds) of the oxide lattice in cubic ice. 
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lar or tetrahedral units. However, since the pyrochlore 
lattice has cubic symmetry, the x, y and z directions are 
equivalent, and this renders Anderson's global antiferro- 
magnetic Ising model unreahstic. It is the local náture 
of the quantization direction that is crucial, and which is 
the origin of the frustration for ferromagnetic interactions 
and for the "ehmination" of the frustration for antifer- 
romagnetic exchange. To see this, consider the fohowing 
toy-model Hamiltonian 



H = —J ^ ^ Si ■ S j 

(i,3) 



(10) 



with classical Heisenberg spins Si on the sites i of 
the pyrochlore lattice, interacting via nearest-neighbor 
exchange coupling J. The second term is a single ion 
anisotropy interaction with A the anisotropy parame- 
ter and Zi a unit vector in the local (111) direction 
at site i. For J = O and A > O, the energy is 
lower if Si points along i^, and one refers to this as 
an Ising anisotropy. The case A < O would be referred 
to as an XY model (Bramwell et a/., 1994; Champion 
and Holdsworth, 2004). However, as discussed in Sec- 
tion V. B, the real microscopic crystal field Hamiltonian 
is more complicated than that considered by Bramwell 
et ai (1994) and Champion and Holdsworth (2004) where 
the limit A ^ oo was taken. 

Hence, we assume an extréme Ising limit. A/ J oo. 
Si is then confined to be either parallel or antiparallel to 
Zi. To implement this energetic single- ion constraint, we 
write Si = (Ji'Zi. Injecting this back into H above, we 
obtain 



(11) 



where we have taken Zi ■ Zj = —1/3 for two distinct cu- 
bic [111] directions. One sees that for ferromagnetic J 
(J > 0), the now "global" a^' Ising variables map onto 
an equivalent Ising antiferromagnet (Anderson's model, 
Anderson (1956))with coupling J/3, and is therefore frus- 
trated. Conversely, for antiferromagnetic J (J < 0), the 
minimum energy state consists of all spins pointing in 
(say crf = +1 V i) or ah out (af = -1 V i). 

As will be discussed in Section X. A, Dy2TÍ207 is 
also identified as a spin ice materiál. It turns out 
that both Ho^+ and Dy^+ carry a sizeable magnetic 
moment, /i, of approximately 10/íb in the crystal field 
ground state of H02TÍ2O7 and Dy2TÍ207. Thus, these 
Systems have a magnetostatic dipole-dipole interaction 
D = /ioM^/(47rrnn^) of approximately D ^ lÄ K at 
the nearest-neighbor dištance rnn- Since D is approxi- 
mately the same as ^cw, it is surprising that the long- 
range and complex náture of the dipolar interactions do 
not lift the spin ice degeneracy and drive the systém to 
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FIG. 42 Dipolar spin ice model: the scaled maximum eigen- 
values, A'^"^(g)/D, in the (hhl) piane. The dipole-dipole in- 
teractions are treated with the Ewald approach. Here, the 
exchange coupling J7nn was set to zero. While the spectrum is 
very fiat, suggesting a weak propensity towards long-range or- 
der, a maximum at a unique ordering wave vector Qmax = 001. 



long range order at a critical temperature Tc ^ D. In 
fact, numerical and theoretical studies have compellingly 
demonstrated that it is precisely the mathematical form 
of the dipole-dipole interactions that is at the origin of 
the spin ice phenomenology in rare-earth pyrochlores. In 
particular, it is the ferromagnetic character of the dipo- 
lar interactions at the nearest-neighbor dištance, and not 
a ferromagnetic nearest-neighbor exchange (which turns 
out to be antiferromagnetic, as anály sis of experimen- 
tal data on H02TÍ2O7 and Dy2TÍ2 07 have revealed) that 
is primarily the source of frustration. Indeed, as dis- 
cussed above, if it were not for the dipolar interactions, 
the nearest-neighbor antiferromagnetic interactions alone 
in H02TÍ2O7 and Dy2TÍ207 would drive the systém into 
a long-range ordered phase. 

The minimal model, called the dipolar spin ice model 
(DSM), that is needed to investigate these questions in- 
cludes nearest-neighbor exchange (first term) and long 
range magnetic dipóle interactions (second term) in. 
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For the open pyrochlore lattice structure, we expect 
further neighbor exchange to be very small, so these can 
be neglected as a first approximation. Recent Monte 
Carlo simulations seem to confirm this expectation (Ruff 
et al.^ 2005; Yavors'kii et al.^ 2007). Here the spin vec- 
tor S^' labels the Ising moment of magnitude \S^'\ = 1 
at lattice site ž, oriented along the local Ising (111) axis 
Zi. The dištance Ir^J is measured in units of the near- 
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est neighbor dištance, rnn- J represents the exchange 
energy and D = {(i^jŕ / {Aiľr^^. Because of the local 
Ising axes, the effective nearest-neighbor energy scale is 

-^eff = -^nn + ^nn where Jnn = J/3 and I^nn = 51^/3, 

since Zi ■ Zj = —1/3 and {zi ■ rij){rij • Zj) = —2/3. If 
I^nn = O and J > O one obtains the spin ice model orig- 
inahy proposed by Harris et al. (1997, 1998) henceforth 
referred to as the "nearest-neighbour spin ice model" . 

The condition Jeff > O is a simple criterion to assess 
whether a systém displays a spin ice state. Mean field 
theory (Gingras and den Hertog, 2001) and Monte Carlo 
simulations (Melko et al.^ 2004) íind a critical valne for 
the transition between "all-in/all-out" Néel order and 
spin ice at Jnn/^nn^ -0.901, hence quite close to the 
naive nearest-neighbor estimate Jnn/^nn = —1- The 
dipolar interactions beyond nearest-neighbors provide a 
weak, extra stabilization of the Néel phase over the spin 
ice state. 

The success of the nearest-neighbor criterion in assess- 
ing whether a systém should display spin ice phenomenol- 
ogy or not indicates that, to a large extent, dipolar in- 
teractions beyond nearest-neighbors are self-screened, as 
originally suggested by a mean field calculation (Gingras 
and den Hertog, 2001) and Monte Carlo simulations (den 
Hertog and Gingras, 2000). The reason why dipolar spin 
ice Systems obey the ice rules is not immediately appar- 
ent. Isakov et al. (2005) have proposed that the ice rules 
result from the fact that the dipolar interactions are, up 
to a perturbatively small and rapidly decaying function, a 
real-space projector onto the manifold of ice-rule obeying 
States. A possibly simpler explanation has recently been 
proposed which involves separating each point dipóle into 
its constituent magnetic (monopole) charges and requir- 
ing that each tetrahedron is, magnetic chargé wise, neu- 
tral, which automatically leads to the conclusion that all 
ice-rule obeying states are, again up to a small correc- 
tion, ground states of the dipolar interactions (Castel- 
novo et al.^ 2007). Such a real-space argument is pos- 
sibly not unrelated to the Ewald energy calculations in 
which the dipolar lattice sums are regularized by effective 
Gaussian charges either in reciprocal space (Enjalran and 
Gingras, 2003) or in direct space (Melko et al.^ 2004). 



A. Dy2TÍ207 and H02TÍ2O7 

The Dy^+ ion in Dy2TÍ207 was first identified as hav- 
ing strongly anisotropic properties from the maximum in 
the DC susceptibility (Blôte et al, 1969) at T - 0.9 K 
followed by a precipitous drop to essentially O at T ~ 0.4 
K. This suggested a strongly anisotropic effective g ten- 
sor with di g± ^ O component perpendicular to the local 
(111) axis, making Dy^+ in effect a (111) Ising model. 
This was confirmed by Flood (1974) from direct measure- 
ment of the magnetic moment and from analysis of inelas- 
tic neutrón scattering data by Rosenkranz (Rosenkranz 
et a/., 2000). Calculations from crystal field theory, com- 
plemented by fitting crystal field parameters (CFPs) to 
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FIG. 43 Dy2TÍ2 07: (a) Specific heat and (b) entropy versus 
temperature, measured by Ramirez et al (1999). The recov- 
ered entropy at 10 K agrees reasonably well with Pauling's en- 
tropy, i?/2 ln(3/3). The data are compared to that calculated 
by Monte Carlo simulations of den Hertog and Gingras (2000) 
for the dipolar spin ice model with exchange Jnn = —1.24 K 
and dipolar coupling Dnn = 2.35 K. (c) H02TÍ2O7: The total 
specific heat is shown by the empty squares and the expected 
nuclear contribution by the solid line. The electronic contri- 
bution has been estimated by subtracting these two curves 
(open circles). Near to 0.7 K, this subtraction is prone to a 
large error (see text). Dipolar spin ice simulation results are 
indicated by the filled circles (from Bramwell et al (2001a)). 



susceptibility data also confirmed the strong Ising náture 
of Dy^+ in Dy2TÍ207 (Jana et al, 2002). Specifically, the 
single ion electronic ground state of Dy^+ is a Kramers 
doublet of almost pure |J = 15/2, mj = ±15/2) sepa- 
rated from the first excited state by a gap A 33 meV 
(~ 380 K) while the rescaled CFPs of Rosenkranz et al. 
(2000) and Jana et al. (2002) finds a much smaller exci- 
tation gap A ^ 100 cm-^ (- 140 K). 

Blôte et al. (1969) measured the specific heat C {T) 
of Dy2TÍ207 between 0.37 K and 1.3 K and found only 
a broad maximum around 1.2 K and no sign of a sharp 
specific heat feature indicating a phase transition to long- 
range magnetic order. Most interestingly, a total mag- 
netic entropy of only |Rln(2) was found. Hence, a no- 
ticeable fraction of Rln(2) expected for a Kramers dou- 
blet was missing. 

Soon after the proposal of spin ice in H02TÍ2O7 (Harris 
et al, 1997, 1998a), Ramirez et al. (1999) showed that the 
missing entropy in Dy2TÍ207 could be determined rather 
precisely by a measurement of the magnetic specific heat 
C{T) between 0.4 K and and 12 K. Figúre 43 shows the 
temperature dependence of C{T)/T. The magnetic en- 
tropy change, ASi^2^ between temperatures Ti and T2 
can be found by integrating C (T) /T between these two 
temperatures: 
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A^i,2 = S{T2) - S{T,) = / ' . (13) 

Figúre 43b shows that the magnetic entropy recovered 
is about 3.9 J mol~^ a value that falls considerably 

short of Rln(2) ^ 5.76 J mo\-^ K'^. The difference, 1.86 
J mol~^ is close to Pauhng's estimate for the residual 
extensive entropy of water ice: (R/2) In (3/2) ^ 1.68 J 
mol~^ thus providing compehing thermodynamic 

evidence for the existence of an ice-rules obeying state in 
Dy2TÍ207. 

As mentioned above, the large 10/íb moments of both 
Dy^+ and Ho^+ lead to a critical role for magnetic dipole- 
dipole interactions in spin ices. This energy, I^nn, at 
nearest-neighbor distances can be estimated from the 
effective Curie constant or from the single-ion crystal 
field doublet wavefunctions. The earher theoretical stud- 
ies (den Hertog and Gingras, 2000; Siddharthan et aL, 
1999) of the DSM estimated 

for both H02TÍ2O7 and Dy2TÍ207. Here r^^ = 
(ao/4)\/2 is the nearest-neighbor dištance and a is the 
size of the conventional cubic unit ceh. As discussed 
above, the factor 5/3 originates from the orientation of 
the Ising quantization axes relative to the vector direc- 
tion rnn that connects nearest-neighnor magnetic mo- 
ments. As we briefly discuss below, the current esti- 
mate on Dnn for Ľ)y2TÍ2 07 ís probably accurate to within 
8%— 10%. Hence, the nearest-neighbor exchange, Jnn, is 
the main unknown. It can be estimated from the high- 
temperature (paramagnetic) regime of the magnetic sus- 
ceptibihty, (Siddharthan et al.^ 1999) or of the specific 
heat, C{T)^ (Jana et al.^ 2002) data. A different ap- 
proach, fohowed by den Hertog and Gingras (2000) and 
Bramweh et al (2001a) has been to determine Jnn by 
íitting either the height of the specific heat, Cp, peaks 
near 1 K or the temperature at which the peaks occurs, 
Tp, against Monte Carlo simulations of the DSM. Inter- 
estingly, fits of Tp or Cp allow for a consistent deter- 
mination of Jnn ~ -1.24 K for Dy2TÍ207. Figúre 43 
shows the good agreement between Monte Carlo results 
and experimental data (den Hertog and Gingras, 2000). 
Note here that the parameter D^n sets the scale for the 
dipolar interactions, at /ioM^/(47rrnn^); the simulations 
themselves use true long-range dipole-dipole interactions 
implemented via the Ewald method (Melko et a/., 2004). 
These results show convincingly the spin ice phenomenol- 
ogy in Dy2TÍ2 07, and also in H02TÍ2O7, as we now dis- 
cuss. 

While H02TÍ2O7 was the first compound to be 
proposed as a spin ice, specific heat measurements 
proved initially less straightforward to interpret than in 
Dy2TÍ2 07, and this led to some confusion. Specifically, 



the rapid increase of the specific heat below 1 K was 
originally interpreted as an indication of a phase tran- 
sition to a partially ordered state around a tempera- 
ture of 0.6 K (Siddharthan et a/., 1999). Instead, it 
turns out that the anomalous low-temperature behav- 
ior of the specific heat in H02TÍ2O7 is of nuclear origin. 
An anomalously large hyperfine interaction between the 
electronic and nuclear spins for Ho commonly leads to a 
nuclear specific heat Schottky anomaly at 0.3 K. A sub- 
traction of the nuclear contribution from the total low- 
temperature specific heat reveals the purely electronic 
specific heat (Bramwell et al, 2001a), C {T). The in- 
tegration of C (T) /T from 300 mK up to 10 K gave a 
magnetic entropy deficit of an amount close to Pauling's 
i?/21n(3/2) zero-point entropy, hence confirming, ther- 
modynamically, that H02TÍ2O7 is indeed a spin ice (Cor- 
nelius and Gardner, 2001). Following the same proce- 
dúre as the one used for Dy2TÍ207 (den Hertog and Gin- 
gras, 2000), a comparison of C {T) with Monte Carlo 
simulations allows to estimate the exchange constant 
in H02TÍ2O7 as Jnn ^ —0.55 K, an antiferromagnetic 
value (Bramwell et al.^ 2001a). 

While specific heat measurements on H02TÍ2O7 are 
problematic, this is emphatically not so for neutrón scat- 
tering experiments. Unlike dysprosium, holmium has 
only one stable isotope whose neutrón absorption cross 
section is negligible. A comparison of the experimental 
neutrón scattering intensity with that calculated from 
Monte Carlo simulations of the dipolar spin ice model 
with an exchange constant Jnn ^ —0.55 K determined 
as above shows excellent agreement. Interestingly, both 
the experiment and Monte Carlo data differ substantially 
from that calculated for the nearest-neighbor spin ice 
model (see Fig. 44). This clearly shows that non-trivial 
spin correlations develop in the materiál as it progres- 
sively freezes within the low-temperature spin ice regime. 
Indeed, those correlations are the precursors of those that 
would ultimately lead to long-range or der, if not pre- 
cluded by spin freezing (Melko et a/., 2001, 2004). 

A puzzling question raised by the good agreement be- 
tween Monte Carlo simulations of the DSM and the ex- 
perimental results illustrated in Fig. 43 is: why don't the 
dipolar interactions drive a transition to long-range order 
at a critical temperature Tc ~ I^nn ^ 2 K? A partial an- 
swer can be found in the mean field theory calculations 
of Gingras and den Hertog (2001). It was found there 
that the q dependence of the softest branch of critical 
modes in the dipolar spin ice model is very weakly dis- 
persive, reaching a global maximum eigenvalue A(qniax) 
at qmax = 001. At the mean field level, this indicates 
that a transition to long range order should occur at a 
critical temperature Tc = A (qmax), with the development 
of delta- function Bragg peaks below Tc. 

First of all, this transition does not occur at Tc ~ I^nn 
because the soft (critical) mode at qmax is in "entropic" 
competition with all the other quasi-soft modes at q 7^ 
qmax- The correlated spin ice regime, which one may 
dub a "collective paramagnet" in Villain's sense (Vil- 
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FIG. 44 H02TÍ2O7: Neutrón scattering in the hhl piane show- 
ing experimental data (upper panel; the sharp spots are nu- 
clear Bragg scattering with no magnetic component), com- 
pared with Monte Carlo simulations of the near neighbour 
spin ice model (middle panel) and the dipolar spin ice model 
(lower panel) (Bramwell et aL, 2001a). Blue indicates the 
weakest and red-brown the strongest intensity. 



lain, 1979), albeit with utterly sluggish spin dynamics 
for T < 1 K in Dy2TÍ207, has so much entropic disorder 
that it is not energetically favorable to localize the sys- 
tém in phase space to a long-range ordered state. In the 
Standard single spin-flip Monte Carlo simulations, this 
transition is not observed because the probability to flip 
a spin once the systém enters a state where each tetra- 
hedron obeys on average the "two-in/two-out" ice rule is 
very small and decreases experimentally very fast with 
decreasing temperature (Melko et al.^ 2004). The low- 
energy excitations deep in the spin ice state correspond 
to nonlocal closed loops of spins flipping from "in" to 
"out" and vice verša so that the systém, as it experiences 
those excitations, remains in a spin ice state (Barkema 
and Newman, 1998). Using such loop excitations, Monte 
Carlo simulations (Melko et al.^ 2001, 2004) found a tran- 
sition to the long-range order predicted by mean field 
theory (Gingras and den Hertog, 2001). In Monte Carlo 



simulations of the DSM, a strongly íirst order transition 
occurs at a critical temperature Tc ^ 0.07I)nn where 
all the residual Pauling entropy is recovered through 
the pre-transitional build-up of correlations and, mostly, 
via a large latent heat (Melko et al, 2001, 2004). For 
^nn ^ 2.35 K believed appropriate for Dy2TÍ207 (and 
H02TÍ2O7), this Tc amounts to 160 mK. However, to this 
date, no experimental work has observed a transition to 
long-range order in spin ice materials down to 60 mK (see 
for example, Fukazawa et al. (2002a)). A common expla- 
nation for this absence of a transition in real spin ice ma- 
terials is that equilibration is lost in the spin ice state (e. g. 
T < 0.4 K in Dy2TÍ207 and T < 0.6K in H02TÍ2O7), 
with the real materials not "beneíitting" from nonlocal 
dynamics as employed in the simulations (Melko et al.^ 
2001, 2004). However, this explanation is certainly some- 
what incomplete since, as we discuss further below, a 
number of experiments report spin dynamics down to 20 
mK. However, before we discuss experiments investigat- 
ing the dynamics of spin ices, we briefly comment on the 
spin-spin correlations in the spin ice regime of Dy2TÍ207. 

Since the orignal work of Ramirez et al. (1999), other 
measurements of the magnetic specific heat data of 
Dy2TÍ2 07 have been reported (Higashinaka et al.^ 2002, 
2003; Hiroi et al, 2003; Ke et al, 2007). Measurements 
of the magnetic specific heat of magnetic insulators are 
notoriously difficult experiments owing to the poor sam- 
ple thermal conductivity. It is perhaps for this reason 
that specific heat data of Dy2TÍ207 from different mea- 
surements show some disparity. For example, the height 
of the magnetic specific heat at its maximum, Cp, for a 
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FIG. 45 Dipolar spin ice model: the phase diagram. The 
antiferromagnetic ground state is an all-spins-in or all-spins- 
out configuration for each tetrahedron. The spin ice con- 
figuration, which includes the q = 001 ground state, is a 
two spins in-two spins out configuration for each tetrahedron. 
The región encompassed between the quasi vertical dotted 
Unes displays hysteresis in the long-range ordered state se- 
lected (q = 000 vs. q = 001) as Jnn/^nn is varied at fixed 
temperature T (Melko et al (2004)). 
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FIG. 46 Dy2TÍ207: ac-susceptibility (upper panel: x^ lower 
panel x") ^ function of temperature at several frequencies, 
illustrating the "15 K peak" from Matsuhira et al (2001). 



powder sample (Ramirez et al.^ 1999) differs by as much 
as 10% compared to that reported in Higashinaka et al 
(2003) for a single crystal. Ultimately, precise measure- 
ments of C (T) are needed if one wants to make an ac- 
curate determination of the residual entropy in the sys- 
tém (Ke et al, 2007). 

The Monte Carlo prediction of long range order devel- 
opping at low temperature in the dipolar spin ice model 
and the experimentally observed collapse of the magnetic 
specific heat below 0.4 in Dy2TÍ2 07 raise the question of 
low temperature spin dy namies in spin ice mater ials. In 
fact experimental studies of this question have led to the 
observation of a much richer phenomenology than that 
which would have been naively expected. 

Measurements of the AC magnetic susceptibility x(^) 
of Dy2TÍ207 down to 60 mK (Fukazawa et al, 2002a) 
and 100 mK (Matsuhira et al, 2001) find that the real 
part of X^ drops precipitously below a temperature 
of roughly 1 K for a frequency of the order of 10 Hz 
(see Fig. 46). At the same time, the imaginary part, x'^ 
shows a rounded maximum. Both x' and x" remain es- 
sentially zero below 0.5 K down to the lowest temperature 
considered, hence signalling an essentially complete spin 
freezing of the systém. Thus, no signatúre of the transi- 
tion to long range order predicted by numerical simula- 
tions (Melko et al, 2001, 2004) is observed. 

The behavior of x' ^ind x" near 1 K signals a spin 
freezing process analogous to what is observed in spin 
glasses. Indeed, measurements of the magnetization of 
Dy2TÍ2 07 show irreversibilities between zero field cool- 
ing (ZFC) and field cooling (FC) below 0.7 K (Matsuhira 
et al, 2001; Snyder et al, 2004). Similar FC-ZFC irre- 



versibilities occur in the Ho2Sn2 07 spin ice (Matsuhira 
et al, 2000). An analysis of the temperature dependence 
of the frequency f m at which x" displays a peak reveals a 
sort of thermally activated freezing behavior which was 
originally parametrized by an Arrhenius form with an 
activation energy of approximately 10 K. However, Sny- 
der et al (2004) questioned the application of the Ar- 
rhenius law to these data. On the other hand, Monte 
Carlo simulations that employ a štandard single spin- 
flip Metropolis algorithm find that the fraction of ac- 
cepted spin flips decreases with decreasing temperature 
according to a Vogel-Fulcher form exp[— A/ (T— T*)] with 
T* - 0.3 K (Melko et al, 2004). The dynamical freezing 
behavior seen in Dy2TÍ207 differs from the critical slow- 
ing down observed in conventional disordered spin glass 
materials (Binder and Young, 1986). Interestingly, for 
the lowest temperature considered, x'' is cut-off on the 
low- frequency regime (Snyder et al, 2004), reminding one 
of what is observed in the LiRoxYi-x^4 Ising systém — 
a phenomenon that has been referred to as "antiglass" 
behavior (Ghosh et al, 2002; Reich et al, 1987). An- 
other difference between the spin freezing in Dy2TÍ207 
and that in conventional disordered spin glasses is the 
magnetic field dependence of Tf. In spin glasses, T f de- 
creases with increasing magnetic field strength, while in 
Dy2TÍ207, the opposite is seen for applied magnetic fields 
up to 5 kOe (Snyder et al, 2004). Measurements on sin- 
gle crystals (Shi et al, 2007) find that the spin freezing 
has a stronger frequency and magnetic field dependence 
for a field along the [111] axis compared to [100], and the 
starting freezing frequency of the single crystal is higher 
than that of the powder sample (Snyder et al, 2001). 
While a quantitative understanding does not exists, the 
behavior at T < 4 K is qualitatively interpreted as a 
signatúre of the collective freezing as the systém enters 
the low-temperature state where the "two-in/two-out" 
ice rules become fulfilled. 

In addition to the low-temperature freezing in the ice- 
rule obeying state, two experimental studies coinciden- 
tally reported results from AC susceptibility x(^) mea- 
surements in Dy2TÍ2 07 above 4K (Matsuhira et al, 
2001; Snyder et al, 2001), finding another freezing pro- 
cess in a temperature range around 15 K, referred to as 
the '15 K feature'. The signatúre of this freezing is only 
seen in the AC susceptibility at finite frequency and not 
in the DC susceptibility data (see Fig. 46). The maxi- 
mum in the imaginary part, x'^ is at 12 K for a frequency 
of 10 Hz increasing to 17 K for a frequency of 10 kHz. 
While the raw data of Snyder et al (2001) and Matsuhira 
et al (2001) are similar, they were analyzed somewhat 
differently. Snyder et al (2001) characterized the freez- 
ing via a single exponential relaxation while Matsuhira 
et al (2001) characterized the AC susceptibility using the 
so-called Davidson-Cole framework, based on an under- 
lying distribution of time scales. Notwithstanding these 
differences, both analyses agreed that there exists a typ- 
ical time scale, r (T), for this freezing phenomenon that 
is parametrized by an Arrhenius form, r ~ ex.p{Ea/T), 
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with an activation barrier energy Ea of the order of 200 
K. A recent /iSR study also finds that the relaxation rate 
of the muon spin polarization in a temperature range 
of 70 K to 280 K can be described by a typical relax- 
ation rate A(T) - exp{-Ea/T) with K - 220 K (Lago 
et al.^ 2007). The exp(£^a/^) dependence of the typical 
time scale characterizing the dynamics around 15 K and 
the energy scale Ea ^ 220 K indicates that the relax- 
ation involves transitions to and from the first excited 
doublet which constitute the main contribution to the 
spin dynamics in the temperature range above 20 K or 
so. While this freezing phenomenon and energy scales 
characterizing it suggest an Orbach process (Finn et al.^ 
1961; Orbach, 1961) that involves both the lattice degrees 
of freedom and the excited crystal field states, a concrete 
microscopic calculation has not yet been done. 

One interesting aspect of the 15 K feature is its be- 
havior when Dy^+ is substituted by diamagnetic in 
Dy2-ccYa;TÍ207. In particular, initial studies found that 
the 15 K peak in x"{^) disappears hy x = 0.4 (Snyder 
et a/., 2001). This was originally interpreted as a sign 
that the 15 K feature is of collective origin. However, in 
a subsequent study (Snyder et al.^ 2004a), it was found 
to re-emerge as x is further increased and is almost as 
strong for x = 1.98 as it is for x = O, but, interest- 
ingly, repositioned at a higher temperature of 22 K for 
a frequency of 1 kHz. This high temperature freezing 
feature is essentially a single-ion phenomenon, akin to 
superparamagnetic spin blocking. While this is not seen 
in AC measurements on H02TÍ2O7, it can be revealed by 
the application of a magnetic field (Ehlers et al.^ 2003). 
Neutrón spin echo experiments on the same materiál con- 
firms the single-ion náture of the 15 K feature (Ehlers 
et al, 2003). 

Another noteworthy aspect of the spin dynamics in 
Dy2TÍ2 07 is the temperature independence of the re- 
laxation time r between 5 K and 10 K (Snyder et al.^ 
2003). Below 5 K, the relaxation time becomes again 
sharply dependent on temperature upon approaching the 
spin ice freezing discussed above. This temperature in- 
dependence of r has been interpreted as a quantum tun- 
neling effect between the up and down Ising spin states. 
This was first observed in H02TÍ2O7 using neutrón spin 
echo (Ehlers et al.^ 2003). Such temperature-independent 
relaxation in Dy2TÍ207 has also been seen in muon spin 
relaxation (/iSR) (Lago et al.^ 2007). However, there is a 
three-orders of magnitude difference between the relax- 
ation rate measured in /iSR and AC susceptibility. Per- 
haps this is because /iSR relies on a local measurement 
that probes all wavevectors of the spin susceptibility. At 
this time, this discrepancy is unresolved. 

A further important topic pertaining to the dynamics 
of spin ice materials is that of the low-temperature spin 
dynamics deep in the frozen spin ice state. As discussed 
above, there is evidence that the electron spin flip dynam- 
ics are exponentially frozen out below 1 K in Dy2TÍ207, 
as well as in other spin ice materials. There is, however, 
some indication of residual spin dynamics in these sys- 
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FIG. 47 The magnetic entropy for Ho2+ccTÍ2-cc07-(5 (inte- 
grated from below T 1 to 22 K) as a function of stuffing at 
H=0 and 1 T. The dashed lines represent the predicted ice en- 
tropy and total spin entropy values. Inset: The temperature 
dependence of the integrated entropy for few compositions at 
H=0 and 1 T (Lau et al, 2006). 



tems that survives down to the lowest temperatures. For 
example, /iSR experiments find a relaxation rate of 0.2 
/is~^ of the muon spin polarization at a temperature of 20 
mK (Lago et al, 2007). This relaxation has been ascribed 
to hyperfine coupling of the electronic and nuclear spins 
which induce a "wobble" around the local (111) Ising di- 
rections. Another work has suggested that the absence 
of a low-temperature nuclear specific heat anomaly in 
Ľ)y2TÍ2 07 may indicate that the electron spin dynamics 
persist to the lowest temperature (Bertin et al.^ 2002). 
This argument would suggest that H02TÍ2O7, which has 
a fully developed nuclear specific heat (Bramwell et al.^ 
2001a), is completely státie. More experimental and the- 
oretical work is required to fully understand the low- 
temperature dynamics of the spin ices. 

Above, we touched on what role magnetic dilution 
plays on the spin dynamics and freezing phenomenon in 
Ľ)y2TÍ207. A recent magnetocaloric study of Dy2TÍ207 
has found a crossover at a temperature of about 0.3 K to a 
low-temperature regime characterized by extremely slow 
relaxation (Orendác et al.^ 2007). In addition, a dilution 
of Dy by 50% of non-magnetic Y, giving DyYTÍ2 07, leads 
to an increase of the relaxation time compared to pure 
Ľ)y2TÍ2 07. This is in contrast with the behavior seen 
above in the formation of the ice state, say above 2 K, 
where a nontrivial dependence of the relaxation time as 
a function of Dy concentration is observed. In particular, 
a level of magnetic dilution less than 12% was found to 
accelerate the relaxation rate while a dilution level higher 
than 12% was found to slow it down again, such that the 
relaxation rates are nearly the same for DyYTÍ207 and 
Dy2TÍ2 07 (Snyder et al, 2004a). 

This discussion illustrates that the náture of the dy- 
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namies in spin ices and the role magnetic dilution is still 
poorly understood. In the context of magnetic dilution, 
we note that a recent paper reports a non-monotonic 
dependence of the residual entropy in Dy2-a:Ya^TÍ207 
and Ho2-a;Ya;TÍ2 07. The data are qualitatively ex- 
plained by a generalization of Pauling's theory for the 
entropy of ice that incorporates site-dilution (Ke et al.^ 
2007). The topic of residual entropy in spin ice when 
the magnetic species are diluted is an interesting prob- 
lém. One would naively expect that the extensive de- 
generacy would ultimately be lifted by disorder (Villain, 
1979). There are two obvious ways that one can en- 
visage this happening. Perhaps the most interesting is 
where dilution would lead to conventional long range or- 
der. For example, in water ice, KOH doping (which ef- 
fectively removes protons) leads to a long-range ordered 
phase of ice called ice XI. The protón vacancies cre- 
ated in the protón structure enhance the dynamics so 
that the systém can develop long range order. Studies 
of diluted spin ice, such as Dy2-x^ (Ke et al.^ 
2007) and Ho2-:cY^TÍ207 (Ehlers et a/., 2006a), have 
so far not reported any signs that diamagnetic dilution 
leads to long range order. Most likely, the extensive 
degeneracy is lifted at high x with the systém leaving 
the ice regime to become a dipolar Ising spin glass akin 
to LiHoa;Yi_a;F4 (Reich et al.^ 1987). Indeed, the sys- 
tematic study of the development of dipolar spin glass 
in both Dy2-a:Ya^TÍ207 and Ho2-a^YajTÍ207 might help 
shed light on the paradoxical antiglass phenomenon in 
LiHoa:Yi_^F4 (Ghosh et al, 2002; Reich et al, 1987), or 
even the possible absence of a spin glass phase altogether 
in that materiál (Jônsson et al.^ 2007). 

On the theoretical front, the DSM of Eq. (12) 
with nearest-neighbor exchange and long range dipole- 
dipole interactions has been fairly successful in describ- 
ing semi-quantitatively the thermodynamic properties 
of Dy2TÍ2 07 both in zero (den Hertog and Gingras, 
2000; Fennell et al, 2004) and in nonzero magnetic 
field (Fukazawa et al.^ 2002a; Ruff et al.^ 2005; Tabata 
et a/., 2006). A detailed comparison between experimen- 
tal results from measurements and Monte Carlo simula- 
tions (Ruff et al.^ 2005) for a field along the [112] direc- 
tion provides strong evidence that exchange interactions 
beyond nearest-neighbors are required to describe quan- 
titatively the experimental data (see also Tabata et al. 
(2006)). This may come as a surprise given the already 
good agreement that was first reported between Monte 
Carlo simulations of the DSM (den Hertog and Gingras, 
2000) and the specific heat measurements of Ramirez 
et al. (1999). However, since that work, as mentioned 
above, sever al other zero field specific heat data sets have 
been reported (Higashinaka et al, 2002, 2003; Hiroi et a/., 
2003; Ke et a/., 2007) and these are no longer so-well de- 
scribed by the originál DSM. A systematic study of more 
recent specific heat measurements in zero and nonzero 
field and magnetization measurements for [110] and [111] 
fields have estimated first (Ji), second ( J2) and third ( J3) 
nearest-neighbor exchange parameters in Dy2TÍ207 (Ya- 



vors'kii et al.^ 2008). Perhaps most interestingly, such a 
refinement of the spin Hamiltonian seemingly allows one 
to explain from a microscopic basis the diffuse scatter- 
ing on the Brillouin zóne boundary (Fennell et al.^ 2004), 
akin to the highly structured inelastic features found in 
the ZnCr204 spinel (Lee et al, 2002). 



B. ^2811207 {A = Pr, Dy and Ho) 

Other pyrochlore oxides with similar properties include 
IIo2Sn207, Dy2Sn207 and Pr2Sn207 (Matsuhira et al.^ 
2000, 2002, 2004). These have not been studied in single 
crystal form, as far as we know, but bulk magnetic prop- 
erties would suggest they are ferromagnetic with large 
Ising anisotropy and very slow dynamics. 

Of these three compounds, only Ho2Sn207 has been in- 
vestigated using neutrón scattering techniques on powder 
samples (Kadowaki et a/., 2002). In the inelastic spec- 
trum they found very slow dynamics below 40 K and 
crystal field levels at 22 and 26 meV, suggesting the lo- 
cal Ho environment is very similar to that of II02TÍ2O7. 
Analysis of the diffuse magnetic scattering also led the 
authors to conclude that Ho2Sn207 was a dipolar-spin- 
ice materiál. 



C. Ho2+ccTÍ2-cc07-(5: stuffed spin ice 

Besides thermodynamic variables such as temperature 
and magnetic field, the chemical composition of the sys- 
tems may be altered to study the relationship of struc- 
tural and magnetic properties. To this end, 'stuffeď spin 
ices with generál formula A2-\-x^2-x07-s have recently 
been synthesized (Lau et al.^ 2006; Ueland, 2008), in 
which additional magnetic ions replace the non-magnetic 
Ti^+ ions ((5 > O, as the oxygen content needs to be ad- 
justed for chargé balance). The entire rare earth serieš 
of titanates have been formed (Lau et a/., 2006a) but 
to date, the only published data is on the 'stuffeď spin 
ice (SSI) compounds. Lau et al (2007) have reported 
phase separation in these compounds, however the mag- 
netic properties are very similar between samples. Single 
crystals have been produced by Zhou et al. (2007). 

The additional magnetic exchange pathways represent 
a major disturbance of the systém, introducing positional 
disorder and, naively, an increased level of energetic con- 
straints to the formation of the spin ice manifold. Sur- 
prisingly, it has been found that the 'stuffeď spin ice sys- 
tems do not freeze, have short range correlations down 
to the lowest temperatures and, most interestingly, have 
the same entropy per spin at low temperature as the 'un- 
stuffeď spin ice (Lau et al, 2006; Zhou et al, 2007). The 
origins of this residual entropy per spin in these systems 
is still debatable. 

Another interesting feature of the residual entropy in 
the stuffed spin ice materials is its robustness to an ap- 
plied field. As seen in Fig. 47 the entropy per spin with 
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FIG. 48 Inelastic neutrón scattering from 
H02.3TÍ1.7O7-5 (Zhou et aL, 2007). Data at 55 mK, 
700 mK and 10 K, noting the appearance of a gapped spin 
excitation at low temperatures. Neutrón spin echo results 
on stuffed (open symbols) and unstuffed (closed symbols) 
spin ice are depicted in the bottom right panel. Within the 
neutrón time window, the spin ice appears státie by 10 K 
(90% of the spin are státie), but the stuffed spin ice has 
persistent dynamics. 



x = 0.67 in 1 T and zero field are identical. This ro- 
bust residual entropy in a very disordered sample like 
H02.67Ti1.33O6. 67 {x = 0.67) needs investigating further. 

Inelastic and quasielastic neutrón scattering on poly- 
crystalline and single crystalline H02. 3X11.707-5 have re- 
vealed subtle changes vis á vis the parent compound. Dif- 
fuse scattering is centered at 0.9 Ä~^, but it is broader 
than the diífuse scattering seen in H02TÍ2O7. Also, as 
seen in Fig. 48, the quasielastic scattering at tempera- 
tures above 10 K opens a gap resulting in a low lying 
excitation, centered at 0.8 meV (Zhou et a/., 2007). 

More research is currently being done by several groups 
on 'stuffed' spin ice. It appears that the stuffed Dy-spin 
ices do not have residual entropy (Ueland, 2008). 



XI. SPIN LIQUID PHASES 

An array of Heisenberg spins forming a three- 
dimensional pyrochlore lattice that interact among them- 
selves via nearest-neighbor antiferromagnetic exchange 
interactions is theoretically predicted to remain disor- 
dered at finite temperature for either classical (Moess- 
ner and Chalker, 1998, a; Moessner, 2001) or quantum 
spins (Canals and Lacroix, 1998, 2000). Ising spins that 
are coupled antiferromagnetically on a pyrochlore lat- 
tice also possess a large ground state degeneracy char- 
acteristic of the spin liquid or cooperative paramagnetic 
state (Anderson, 1956). 



A. Tb2TÍ207 

Tb2TÍ207 is a well studied systém where Tb^+ ions 
are magnet ic and the Ti^+ ions are nonmagnetic. The 
temperature dependence of the magnetic susceptibility 
of Tb2TÍ207 is described well by the Curie-Weiss law 
down to 50 K with (9cw = -19 K and 9.6 /i^/Tb- 
ion. This effective moment is appropriate for the ^F6 
Tb^+ ion. By studying a magnetically dilute sample, 
(Tbo.o2Yo.98)2TÍ2 07, a value of -6 K was established as 
the crystal field contribution to Ocw- Assuming a max- 
imum of -2 K for the dipolar interactions, ^ -11 K was 
proposed as a good estimate for the contribution of the 
exchange interactions to ^cw • Figúre 49a shows the tem- 
perature dependence of the bulk susceptibility down to 
0.5 K. 

Han et ai (2004) found, using neutrón powder- 
diffraction and x-ray absorption fine-structure, that the 
chemical structure of Tb2TÍ2 07 is well ordered with no 
structural transitions between 4.5 and 600 K. More re- 
cently however, Ruff et al (2007) found a broadening of 
structural peaks in a high resolution, single crystal, x-ray 
diffraction experiment below 20 K, as expected just above 
a cubic-to-tetragonal transition, however this transition 
never fully develops even at 300 mK. 

Powder neutrón diffraction measurements on 
Tb2TÍ207 clearly revealed two contributions to the 
scattering at 2.5 K (Gardner et ai, 1999). One of 
these consists of sharp nuclear Bragg peaks due to the 
crystalline order in the materiál and results in a cubic 
lattice parameter of 10.133 Ä. A second diffuse, liquid- 
like, background is also present and was attributed to 
magnetic short range correlations. Gardner et al (1999) 
used an isotropic spin model, correlated over nearest 
spins only, to describe this scattering which captured the 
minimum in the forward direction, |Q| ~ O Ä~^, as well 
as the approximate location of the peaks and valleys in 
the data. 

Inelastic scattering experiments were also carried out 
on polycrystalline Tb2TÍ207 (Gardner et al.^ 2001) and 
a typical scan at |Q| =2 is shown in Fig. 50. Three 
bands of excitations are clearly observed near 0.37, 2.53, 




o 50 100 150 200 250 300 O 0.1 0.2 0.3 0.4 0.5 



T (K) T (K) 

FIG. 49 Left: The temperature dependance of the inverse 
molar susceptibility and the resulting fit to the Curie-Weiss 
law at high temperature (T > 250 K). Right: The low temper- 
ature (T < 0.5 K) dependence of the magnetic susceptibility 
as a function of frequency (Gardner et a/., 2003). 
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FIG. 50 Left: A typical constant |Q| scan at 2 revealing 
low energy inelastic modes at 12 K on a powder sample of 
Tb2TÍ207. Right: A schematic of the crystalline electric field 
levels (Gingras et al, 2000). 



and 3.50 THz and these were shown to be dispersionless 
above 30 K, a characteristic of a single ion effect asso- 
ciated with the rare earth site. A weak, but very inter- 
esting dispersion develops in the lowest energy band at 
temperatures below ~ 20 K, which we wih discuss later 
in this section. 

A crystahine electric field (CEF) level schéme appro- 
priate to the ^Fq configuration of Tb^+ in the A-site 
environment of Tb2TÍ207 was determined on the basis 
of these neutrón and other bulk property data (Gingras 
et al, 2000). The ground state and first excited states 
are both doublets, with two singlets at much higher en- 
ergies, as shown schematically in Fig. 50. The eigen- 
states I i 4) and | ± 5) are believed to make up most of 
the weight of the ground state doublet and of the first 
excited state, respectively. Such a CEF schéme means 
the moment can be considered extremely Ising in náture 
and pointing in the local (111) directions at tempera- 
tures T < 20 K. More recently Mirebeau et al (2007) re- 
ported a serieš of neutrón scattering measurements that 
"refined" the CEF schéme of Tb2TÍ2 07 and compared it 
to Tb2Sn207, which orders just below 1 K, and which 
was discussed in Section VIII. C. 

Large single crystals of Tb2TÍ207 were first grown us- 
ing íioating zóne image furnace techniques by Gardner 
et al. (1998). The resulting single crystals enabled a serieš 
of experiments which could próbe the four-dimensional 
dynamic structure factor S{Q^huj). Figúre 51 shows 
S'(Q) = J S'(Q, hjú)duú measured at T=50 mK and at 
9 K within the (hhl) piane. Bragg peaks are clearly seen 
at the allowed positions for the Fd3m space group, that 
is {hhl) being all even, or all odd, integers. Also one 
clearly observes a "checkerboard" pattern to the mag- 
netic diffuse scattering. This covers the entire Brillouin 
zóne indicating the spins are correlated on a length scale 
much smaller than the unit celí ~ 10.1 Ä. However, the 
single crystal data are clearly not isotropic in recipro- 
cal space, explaining the quantitative failure of the sim- 
ple isotropic model used to analýze the earlier powder 
data (Gardner et al.^ 2001). Also, the scattering does not 



FIG. 51 Diffuse scattering from a single crystal of Tb2TÍ2 07 
at 50 mK (left) and 9 K (right) (Gardner et al, 2001, 2003). 
Sharp Bragg peaks can be seen at the appropriate reciprocal 
lattice positions in the 50 mK data set, but at 9 K, these have 
been subtracted out using a high temperature (100 K) data 
set. 



conform to that expected for an Ising model with local 
(111) spin directions. Specifically, an Ising model with 
local (111) spins cannot simultaneously produce a large 
amount of scattering around 002 and null scattering at 
000 (Enjalran and Gingras, 2004). 

Single crystals also allow the study of the low-lying ex- 
citations as a function of Q, (Gardner et al.^ 2001, 2003) 
as opposed to the modulus |Q| derived from powder sam- 
ples. As was observed in powders (Gardner et al.^ 1999), 
the dispersion develops below ~ 25 K. This low tempera- 
ture dispersion is plotted along the three high symmetry 
directions within the (hhl) piane in Fig. 52, on which are 
overlaid cuts of S'(Q) along [001] (top), [hhO] (middle), 
and [hhh] (bottom). The diffuse scattering is well de- 
scribed by simple near-neighbor antiferromagnetic spin 
correlations on the pyrochlore lattice, as shown by the 
solid lines in Fig. 52. We also see that the minima in 
the dispersion of this low lying magnetic mode (closed 
symbols) correspond exactly with peaks in 6'(Q), even 
though S{Q) is anisotropic in Q. Gardner et al. (2003) 
showed the gap between the first excited state and the 
ground state drops from 0.37 THz at 30 K to approxi- 
mately 0.25 THz at 100 mK but does not soften further 
as the temperature is further decreased. 

Neutrón spin echo (Gardner et al., 2003, 2004) and 
/iSR (Gardner et al, 1999; Keren et al, 2004) mea- 
surements have also been carried out at mK tempera- 
tures to próbe the dynamics of the magnetic moment s. 
Both show large, fluctuating moment s down to at least 
17 mK (Gardner et al, 1999), as expected for a coop- 
erative paramagnet. However, below T ^ 0.3 K the 
5'(Q, t)/6'(Q, 0) baseline is no longer zero indicating a 
proportion (^ 20%) of the magnetic moments are frozen 
in the neutrón spin echo time window. This is consistent 
with the frequency dependence of the AC susceptibility 
shown in figúre 49b (Gardner et a/., 2003; Luo et a/., 
2001). However, we attribute this partial spin freezing to 
a small subset of magnetic moments distributed near a 
small number of defects. Evidence for spin freezing was 
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FIG. 52 Cuts through the checkerboard board pattern of dif- 
fuse magnetic scattering (open symbols) along high symmetry 
directions of the cubic lattice. To complement these data, the 
dispersion of the lowest lying magnetic excitation, at 4 K, is 
also plotted (closed symbols). Fits to the diffuse scattering 
are shown and discussed in the text (Gardner et a/., 2001). 



also found independently in neutrón scattering experi- 
ments performed on other single crystal samples (Yasui 
et al.^ 2002). In this case, hysteresis is observed in the 
scattering near 002 at temperatures less than ~ 1.7 K, 
but again most of the scattering remains diffuse in náture 
even below 1.7 K. 

Dilution studies by Keren et al. (2004), perturbed the 
lattice through the percolation threshold but found that 
fluctuating moments prevail at all magnetic concentra- 
tion as shown by the muon relaxation rate in Fig. 53. 
They also noted that the magnetic coverage of the lat- 
tice was important to the cooperative phenomenon but 
the percolation threshold played no critical role. 

At first sight, it would be tempting to ascribe 
the failure of Tb2TÍ207 to develop long range order 
down to such a low temperature compared to ^cw, 
to the collective paramagnetic behavior of the classical 
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FIG. 53 Temperature dependence of the muon relax- 
ation rate in a field of 50 G in various samples of 
(TbpYi_p)2TÍ207 (Keren et a/., 2004). 



pyrochlore Heisenberg antiferromagnet (Moessner and 
Chalker, 1998, a; Villain, 1979). However, consideration 
of the single-ion crystal field ground state doublet of 
Tb^+ in Tb2TÍ207 strongly suggests that, in absence of 
exchange and dipolar interactions, bare Tb^+ should be 
considered as an effective Ising spin (Gingras et al.^ 2000; 
Rosenkranz et a/., 2000). 

Considering only nearest-neighbor exchange and long- 
range dipolar interactions, the DSM of Eq. (12) would 
be the appropriate Hamiltonian to describe Tb2TÍ2 07 as 
an Ising systém. On the basis of the estimated value 
of the nearest-neighbor exchange and dipolar coupling, 
Tb2TÍ207 would be predicted to develop a /Cord = 000 
four sublattice long-range Néel order at a critical tem- 
perature Tc ~ 1 K (den Hertog and Gingras, 2000), in 
dramatic contradition with experiment al results. Fur- 
thermore, the paramagnetic neutrón scattering of such a 
classical Ising model is inconsistent with that observed 
in Tb2TÍ207 at 9 K as shown in Fig. 51. Mean field 
theory calculations on a Heisenberg model with dipolar 
and exchange interactions and with a simple single-ion 
anistropy (Enjalran and Gingras, 2004) as well as a ran- 
dom phase approximation (RPA) calculation (Kao et a/., 
2003) that considers a semi-realistic description of the 
crystal-field levels of Tb^+ show that the strong scatter- 
ing intensity around 002 should be extremely weak for 
Ising spins. Rather, the intensity is caused by the finite 
amount of fluctuations perpendicular to the local (111) 
directions (Enjalran and Gingras, 2004). In the RPA cal- 
culations (Kao et al.^ 2003), the scattering intensity near 
002 originates from the contribution of the excited crys- 
tal field levels at an energy ~ 18 K to the susceptibility, 
which is made dispersive via the exchange and dipolar in- 
teractions. The experimental observation that the diffuse 
scattering intensity is largely unaltered down to a tem- 
perature of 50 mK, as shown in Fig. 51, suggests that the 
"restored" effective isotropicity of the spins is an intrin- 
sic part of the effective low-energy theory that describes 
Tb2TÍ207. In that context, a recent paper by Molavian 
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FIG. 54 Estimate fluctuation rate of the Yb^^ moment de- 
rived from muon (filled circles) and Môssbauer (open cir- 
cles) spectroscopies (Hodges et a/., 2002). Note that the rate 
plateaus at 200 mK and does not go to zero. 



and Môssbauer spectroscopies, they concluded that the 
Yb^+ spin fluctuations slow down by more than 3 or- 
ders of magnitude to sever al megahertz, without freezing 
completely. This was confirmed in their neutrón powder 
diífraction, where no extra Bragg intensity was reported 
below 0.24 K. Single crystal neutrón diffraction by Yasui 
et al. (2003) has revealed extra Bragg scattering below 
0.24 K from a státie ferromagnetic state, albeit with a 
reduced moment of 1.1 fiB- The latter two studies mo- 
tivated a polarised neutrón scattering study by Gardner 
et al. (2004a). That work conclusively ruled out a frozen 
ferromagnetic state and confirmed that the majority of 
the spin systém continues to fluctuate below the 240 mK 
transition while a small amount of magnetic scattering 
was observed at the 111 Bragg position at 90 mK. 



et al. (2007) argues, on the basis of a calculation that 
considers non-interacting tetrahedra, that virtual fluctu- 
ations between the ground state and excited crystal fleld 
levels lead to an important renormalization of the cou- 
pling constants entering the effective low-energy theory. 
In particular, Molavian et al. (2007) propose that, for the 
simple model considered, the systém is pushed to the spin 
ice side of the phase diagram by those quantum fluctua- 
tions. In other words, within that description, Tb2TÍ2 07 
is a sort of quantum spin ice systém where the spin ice 
like correlations remain hidden down to T < 500 mK. 

In this discussion of the lack of order in Tb2TÍ207, 
it is not clear what is the role of the structural/lattice 
fluctuations very recently reported by Ruff et al. (2007). 
At this time, there is really no robust theoretical under- 
standing of the physics at play in Tb2TÍ207 and more 
theoretical and experimental studies are required. 



B. Yb2TÍ207 

Ytterbium titanate is an insulator with lattice parame- 
ter ttŕ? = 10.026(1) Ä at room temperature (Brixner, 1964). 
Early work suggested an ordered magnetic state just be- 
low 0.2 K (Blôte et al.^ 1969), where a sharp anomaly 
in the speciflc heat was observed. In a detailed mag- 
netisation study, Bramwell et al. (2000) found a Curie- 
Weiss temperature of 0.59(1) K, indicative of weak fer- 
romagnetic coupling and a free ion moment of 3.34(1) 
fiB- Hodges et al. (2001) used Môssbauer spectroscopy 
to investigate the crystal fleld schéme and determined 
that the ground state Kramers doublet was separated 
by 620 K from the flrst excited state producing an easy 
piane anisotropy, like in Er2TÍ2 07. They also found the 
effective paramagnetic moment to be 3.05(8) fiB- 

Hodges and co-workers (Hodges et al, 2002) extended 
their study of Yb2TÍ207 and found an abrupt change 
in the fluctuation rate of the Yb^+ spin at 0.24 K but 
not a frozen ground state as expected from the earlier 
studies (Blôte et al, 1969). Using muon spin relaxation 



C. Er2Sn207 

Er2Sn2 07 and Er2TÍ207 are thought to constitute ex- 
perimental realizations of the highly frustrated XY an- 
tiferromagnet on the pyrochlore lattice. Very little is 
known about Er2Sn2 07 compared to the isostructural 
Er2TÍ207 compound which develops a non-collinear Néel 
state at 1.1 K, possibly via an order by disorder transition 
(see Section VIII. B (Champion et a/., 2003)). Bondah- 
Jagalu and Bramwell (2001) reported a pronounced fleld 
cooled - zero fleld cooled divergence in the susceptibil- 
ity at 3.4 K. The temperature and fleld dependence of 
the bulk magnetization of Er2Sn207 was more recently 
published by Matsuhira et al. (2002). At high temper- 
atures l/x{T) is linear and a flt to the Curie-Weiss law 
results in a ground state moment of 9.59 fiB appropriate 
for the H 15 ion. It also gives a Curie-Weiss temperature, 
Ocw ^ —14 K, suggestive of antiferromagnetic correla- 
tions. At 10 K l/x{T) starts bending downwards, pre- 
sumably due to crystal fleld effects. However, no anomaly 
indicative of a freezing into either a long or short range 
ordered systém was seen down to 0.13 K, suggesting that 
geometrical frustration plays a large role in determining 
the magnetic ground state. /iSR studies by Lago et al. 
(2005) revealed persistent spin dynamics at 20 mK in 
both Er based samples. While the temperature depen- 
dence of the muon relaxation rate is similar, the depo- 
larisation curve for Er2Sn2 07 remains exponential down 
to the lowest temperatures suggesting the systém fails 
to enter an ordered state, while the titanate takes on a 
gaussian relaxation at early times. 

Neutrón diffraction results by Bramwell et al. (2004) 
and Shirai and Bramwell (2007) would however suggest 
that Er2Sn2 07 enters a state with long-range order at 
^ 0.1 K. This ordered ground state is different that 
that of Er2TÍ2 07, which appears to be stable in the 
Er2(Sn,Ti)2 07 solid solution until 90% substitution of 
Ti for Sn. 
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FIG. 55 Top: The temperature dependence of the in- 
verse magnetic susceptibility from a powder sample of 
Er2Sn207 (Matsuhira et aL, 2002). The inset shows the 
data down to 0.13 K. Bottom: Relaxation rates A derived 
from muon spin relaxation measurements for Er2Sn207 and 
Er2TÍ2 07 (Lago et aL, 2005). SoUd curves are the fit of the 
high temperature dat to an exponential activation function. 
A II is the longitudinal relaxation rate for Er2TÍ2 07. Inset: 
thermal evolution of the order parameter critical exponent, 
for Er2TÍ2 07 showing the evolution from exponential to 
Gaussian curve shape across Tn. For Er2Sn207, /3 ~ 1 for 
the entire temperature range. 



D. Pr2lr207 

Recently, much attention has been devoted to the 
metallic materiál Pr2lr207 which is reported to remain 
paramagnetic down to 0.3 K by Machida et al. (2005). 
From magnetization, specific heat and inelastic neutrón 
scattering studies the Pr ground state was determined to 
be a well isolated doublet with Ising-like moments ori- 
ented along the local (111) axes. The antiferromagneti- 
cally coupled spins have an energy scale of 20 K and no 
freezing was observed. Subsequent studies carried out on 
crystals grown from a KF flux showed remarkably com- 
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FIG. 56 D C susceptibility at a field of 20mT applied along 
[111] for Pr2lr207. Note the ln(T) dependence expected for 
the Kondo effect and (inset) evidence for spin freezing (Nakat- 
suji et al, 2006). 



plex behavior, including the Kondo eífect and partial spin 
freezing below 120 mK (see Fig. 56) (Millican et al.^ 2007; 
Nakatsuji et al.^ 2006). The authors have concluded that 
a spin liquid state exists for the Pr moments but no stud- 
ies of spin dynamics have yet been presented to support 
this claim. 

Finally, a report of an unconventional anomalous Hall 
effect by Machida et al (2007) has appeared very re- 
cently and a chirality mechanism, similar to that pro- 
posed for Nd2Mo207 (Taguchi et al, 2001, 2004), has 
been invoked. 



XII. EXTERNAL PERTURBATIONS 

Many studies have been performed on geometrically 
frustrated magnets where an external pertubation (apart 
from nonzero temperature) has been applied to the sys- 
tém. These perturbations can take the form of common 
variables like magnetic field and the chemical composi- 
tion of the systems. Less common perturbations include 
applied pressure and irradiation damage and the effects 
of these too have been investigated on pyrochlore ox- 
ides. To date, no magnetic studies of irradiation damaged 
pyrochlores have been published, though several groups 
have studied the structural properties of pyrochlores af- 
ter being bombarded with radiation (Ewing et al.^ 2004; 
Sickafus et al.^ 2000). Where appropriate, the results 
from changes in temperature or chemical composition 
have been mentioned in the previous sections and we 
will restrict ourselves here to those studies where pro- 
found changes to the magnetic náture of the frustrated 
magnet have been observed. 
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FIG. 57 Magnetic phase diagrams of the Gd2TÍ207 for three 
different orientations of an applied magnetic field (Petrenko 
et al, 2004). These phase boundaries roughly approximate 
those seen by Ramirez et al (2002) on a powder sample. 



A. Magnetic Field 

Most magnetic pyrochlore oxides have been subject to 
an applied field and their behavior monitored. In some 
cases, drastic effects have been observed. Indeed, one 
must remember that these systems are delicately bal- 
anced between sever al competing energy scales, and a 
magnetic field, even small fields, can alter the systém sig- 
nificantly. For example, in the spin ices, H02TÍ2O7 and 
Dy2TÍ2 07, in a field as small as 100 Oe applied along 
[001] or [110], a few of the many ground states are fa- 
vored creating metastable states, magnetization plateaus 
and slow spin dynamics but no elementary changes are 
observed (Fennell et al.^ 2005). In a few cases, the appli- 
cation of a field can drive the systém into a new ground 
state. These facts must be remembered when probing 
frustrated systems, possibly even in the remanent field 
of a magnet. 



1. Gd2TÍ207 

The magnetic field vs temperature phase diagram for 
the Heisenberg pyrochlore magnet Gd2TÍ207 was first 
determined on polycrystalline samples using specific heat 
and AC susceptibility (Ramirez et a/., 2002). A mean 
field model was used to justify the five phases found be- 
low 1 K and in fields up to 8 T. Petrenko et al (2004) 
later reproduced these field induced transitions by mea- 
suring the specific heat along 3 different crystallographic 
directions on a single crystal (see Fig. 57). 

To determine the magnetic structure in these phases, 
neutrón powder diffraction measurements have been per- 
formed. Although significant differences can be seen in 
the diffraction patterns in every phase, only the zero field 
structures have been modeled successfully (Shirai and 
Bramwell, 2007; Stewart et al, 2004). 

With the recent knowledge that Gd2Sn207 is a better 



representation of a Heisenberg pyrochlore magnet with 
dipolar interactions and near neighbor exchange (Del 
Maestro and Gingras, 2007; Wills et al.^ 2006) (see Sec- 
tion VIII. A), one would want to study it in single crys- 
talline form and in a magnetic field. 



2. Kagome ice 

By applying a magnetic field along the [111] direction 
of single crystals of Dy2TÍ2 07 or II02TÍ2O7 one can in- 
duce a new microscopically degenerate phase, known as 
kagome ice. In this state the 2-in-2-out spin configura- 
tion, known as the ice-rules state, is still found, how- 
ever one spin is pinned by the field and the number of 
possible ground states is lowered (Moessner and Sondhi, 
2003). As shown in Fig. 4, the pyrochlore structure as 
observed along the [111] is comprised of stacked, alter- 
nating kagome and triangular planes. When the applied 
field exceeds 0.3 kOe the kagome ice phase is stable and 
exhibits a residual entropy of ~ 40% of that seen in zero 
field spin ice. 

Matsuhira et al. (2002a) performed magnetisation and 
specific heat experiments on Dy2TÍ207, proving the ex- 
istence of this phase. Since then, several pieces of work 
have been done, mainly on Dy2TÍ207, with probes like 
AC susceptibility and specific heat (Higashinaka et al.^ 
2004; Sakakibara et al, 2004; Tabata et al, 2006). Very 
recently Fennell et al (2007) have performed neutrón 
diffraction measurements on H02TÍ2O7 in the kagome 
ice phase and speculate that a Kasteleyn transition takes 
plače. Wills et al (2002) investigated via analytical and 
Monte Carlo calculations an ice-like state in a kagome 
lattice with local Ising anisotropy before one was ex- 
perimentally realized. While the physics of kagome ice 
emerging from a spin ice single crystal subject to a [111] 
field is rather interesting, the behavior of spin ice subject 
to a magnetic field along one of the [110], [112] or [100] 
directions is also worthy of a brief comment. 

The pyrochlore lattice can be viewed as two sets of or- 
thogonal chains of spins. One set is parallel to the [110] 
direction and the other along [110]. These two sets are re- 
ferred to as a and [5 chains, respectively. The application 
of a field parallel to [110] and of strength slightly larger 
than 1 T leads to a pinning of the spins on the a chain, ef- 
fectively freezing those spins. Both the net exchange and 
net dipolar field produced by the frozen spins on the a 
chains vanish by symmetry for the spins on the j3 chains. 
Since the j3 chains are perpendicular to the field along 
[110], the spins along those [5 chains are free to interact 
only among themselves. Just as a [111] field leads to a de- 
composition of the pyrochlore lattice into weakly coupled 
kagome planes with their normál along [111], a field along 
[110] induces a "magnetic break up" of an otherwise cu- 
bic pyrochlore systém into a set of quasi one-dimensional 
(5 chains predominantly weakly coupled by dipolar inter- 
actions. On the basis of Monte Carlo simulations (Ruff 
et al, 2005; Yoshida et al, 2004) and chain mean field 
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theory (Ruíf et a/., 2005), one expects that dipolar spin 
ice would have a transition to a long range ordered state, 
referred to as q = X order (Harris et a/., 1997). However, 
while neutrón scattering experiments on H02TÍ2O7 and 
Dy2TÍ2 07 (Fennell et al.^ 2005) and specific heat mea- 
surements (Hiroi et a/., 2003a) on Dy2TÍ2 07 find some 
evidence for a transition with a field along [110], the tran- 
sition is not sharp. It is not known whether the apparent 
failure of the systém to develop long range q = X order 
is due to either an imperfect ahgnment of the field along 
[110] or the inability of the systém to properly equilibrate 
upon approaching the putative transition. 

The [100] field problém was probably the first one stud- 
ied theoretically (Harris et al.^ 1998). It was predicted 
on the basis of Monte Carlo simulations on the nearest- 
neighbor spin ice model that the systém should exhibit a 
field-driven transition between a gas-like weakly magne- 
tized ice-rule obeying state and a strongly polarized state. 
It was originally argued that this transition is first order 
and terminates at a critical point, similarly to the gas- 
liquid critical point (Harris et a/., 1998). However, recent 
work that combines numerical and analytical calculations 
argues for a more exotic topological Kastelyn transition 
driven by the proliferation of defects in the ice-state as 
the temperature or the field strength is varied (Jaubert 
et a/., 2007). The same work speculates that there ex- 
ists preliminary evidence for such a rounded Kastelyn 
transition in the neutrón scattering data on H02TÍ2O7 
subject to a [100] field (Fenneh et al, 2005). On the 
other hand, while the Monte Carlo data of Harris et al. 
(1998) find a sharp and a broad feature in the specific 
heat as a function of temperature for sufficiently low 
(but nonzero) [100] field, specific heat measurements (Hi- 
gashinaka et al.^ 2003a) on Dy2TÍ2 07 only find one peak 
down to a temperature of 0.35 K. Perhaps long range 
dipolar interactions are of some relevance to the phe- 
nomenology at play for the case of a [100] field. Indeed, 
long range dipole-dipole interactions are definitely crucial 
to the qualitative physics at play for fields along [110] and 
near [112] (Ruff et al.^ 2005). Clearly, more experimental 
and theoretical work on the problém of spin ice in a [100] 
magnet ic field is called for. 



3. TbsTisOr 

Neutrón diffraction and susceptibility data have re- 
cently been collected on the cooperative paramagnet 
Tb2TÍ2 07 in high magnetic fields (Rule et al.^ 2006; Ue- 
land et a/., 2006). In the AC susceptibility work by Ue- 
land et al (2006) unusual cooperative effects and slow 
spin dynamics were induced by applying fields as high as 
9 T at temperatures above 1 K. 

Rule et al (2006) performed single crystal neutrón 
scattering on Tb2TÍ207. They were able to show that 
the short range spin-spin correlations, discussed in Sec- 
tion XI (Gardner et al.^ 2001, 2003), are very dynamic 
in origin and that they become of longer range as a field 
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FIG. 58 Neutrón scattering with AE = ±0.5 meV from 
Tb2TÍ2 07 at 1 K (Rule et al, 2006). In zero field (a), 2 T 
(b) and 8.5 T (c). In panel d) constant Q seans at 0.1 K 
are shown, indicating that the low lying magnetic scattering 
( slow spin dynamics) is quenched as the field increases. 



is applied. Applying a fairly small field of 1000 Oe, at 
400 mK, causes the magnetic diffuse scattering to con- 
dense into peaks (see Fig. 51), with characteristics sim- 
ilar to a polarised paramagnet. Figúre 58 shows similar 
data at 1 K. In zero field, Bragg peaks are seen at 113 
and 222 from the underlying nuclear structure and broad 
diffuse scattering, magnetic in origin is centered at 002. 
In a field up to 2 T (not shown) the diffuse scattering 
sharpens into a 002 peak and in higher fields (/íqH > 2 
T) a magnetically ordered phase is induced, with the de- 
velopment of true long range order accompanied by spin 
waves. 



B. High Pressure 

The study of pressure on samples in polycrystalline 
and single crystal form has become more common in 
recent years, with a paper on the pressure effect on 
the magnetoresistive materiál Tl2Mn2 07 being among 
the earliest (Sushko et al.^ 1996) in pyrochlore physics. 
More recently, studies have been performed on the trans- 
port, structural and magnetic properties of other py- 
rochlores (Mirebeau et al, 2002; Saha et a/., 2006; Zhang 
et al, 2007, 2006). 

Sm2Zr207 undergoes a structural distortion at low 
pressures. However, the distortion is not obser ved be- 
tween 13.5 and 18 GPa. Above 18 GPa, the pyrochlore 
structure is unstable and a distorted fluorite structure, in 
which the cations and anión vacancies are disordered, is 
observed (Zhang et al, 2007). Similar things happen to 
Cd2Nb2 07 at 12 GPa (Samara et al, 2006; Zhang et al, 
2006). Gd2TÍ207 undergoes a subtle distortion of the lat- 
tice at 9 GPa at room temperature (Saha et al, 2006). 
Neutrón diffraction at 1.4 K on H02TÍ2O7 at pressures up 
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to 6 GPa saw no change in the spatial correlations of this 
spin ice compound (Mirebeau and Goncharenko, 2004), 
but a recent study of Dy2TÍ207 observed a small change 
in the magnetisation at 13 kbar (Mito et al.^ 2007). 



1. Tb2TÍ207 

The Tb2TÍ207 cooperative paramagnet has been stud- 
ied by means of single-crystal and polycrystahine neu- 
trón diffraction for an unprecedented range of thermo- 
dynamical parameter s combining high pressures, mag- 
netic fields and low temperatures (Mirebeau et al.^ 2002, 
2004a). Mirebeau et al. (2002) found a long range, mag- 
netically ordered state near 2 K for a hydrostatic pressure 
exceeding ~ 1.5 GPa. Powder neutrón diffraction data 
as a function of pressure are shown in Fig. 59. Weak 
magnetic Bragg peaks occur at positions not associated 
with nuclear reflections although indexable to the Fd3m 
space group. They also found two šmah reflections that 
were attributed to a long range modulation of the main 
structure. 

To investigate the transition further, Mirebeau et al. 
(2004a) studied single crystals where the ordered mag- 
netic moment may be tuned by means of the direction of 
the anisotropic pressure component. It was shown that 
an anisotropic pressure component is needed to suppress 
the spin liquid state, but an anisotropic pressure alone 
with no hydrostatice component produced no effect at 
1.4 K. When an ordered phase was induced, the direc- 
tion of the anisotropic pressure was important with a fac- 
tor of 30 difference seen in the strength of the magnetic 
scattering between different directions (Mirebeau et al, 
2004a). One wonders if the rich behavior of Tb2TÍ2 07 
under pressure is related to the recent observations of 
seemingly dynamical lattice effects in zero applied pres- 
sure (Ruff et al, 2007). 



2. yl2Mo207 {A = Gd and Tb) 

The properties of Gd2Mo207_a; change remarkably 
with pressure. Kim et al. (2005) applied moderate pres- 
sures of 1.6 GPa and found a significant decrease in 
the saturation moment at 4 K with increased hystere- 
sis. Mirebeau et al. (2006a) and Miyoshi et al. (2006) ex- 
tended the pressure range and were able to drive the sys- 
tém into a glassy state at 2.7 GPa. In the study by Mire- 
beau et al. (2006a) the melting of the magnetic structure 
was monitored by neutrón diffraction (see Fig. 60). Thus, 
the application of pressure can change the sign of the 
Mo— Mo exchange interaction from ferromagnetic to an- 
tiferromagnetic. In addition, a phase with mixed Tb/La 
substitution on the A-site, namely(Lao.2Tbo.8)2Mo207, 
was studied by Apetrei et al. (2007a). This mate- 
riál, which has nearly the same unit celí constant as 
Gd2Mo207_a^ [10.3787(8) Ä], shows ferromagnetism with 
Tc = 58 K but with a different magnetic structure in 
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FIG. 59 Pressure induced magnetic Bragg scattering from 
polycrystalline Tb2TÍ2 07 at <2 K. The data shows the pres- 
sure evolution of státie correlations in the cooperative para- 
magnet. The inset depicts the pressure dependence of the 
transition temperature (Mirebeau et a/., 2002). 



which the Tb moments form the familiar two-in two-out 
Ising structure, but there is still ferromagnetic coupling 
between the A and B magnetic sublattices. It is also 
remarked that the ordered moments are reduced signif- 
icantly from the free ion values by 50 percent for Tb, 
30 percent for Gd and 70 percent for Mo. The case for 
Tb is not surprising and can be ascribed to crystal field 
effects, while the result for Mo is in line with previous 
reports (Gaulin et al., 1992; Greedan et al, 1991; Yasui 
et al, 2001). 



Part IV 

Conclusions 

We have highlighted most of the advances made over 
the past 20 years in classifying the magnetic ground 
state in many geometrically frustrated magnets with the 
pyrochlore lattice. An extensive review on other frus- 
trated magnetic systems like the spinels and kagome lat- 
tice magnets would be most welcomed. Although signifi- 
cant advances have been made in the understanding and 
characterization of these pyrochlore magnets, much more 
work can be done. As the systems get more complicated 
(i. e. through dilution of the magnetic sublattice or mix- 
ing of magnetic species) we hope and expect to see more 
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FIG. 60 Pressure dependence on the magnet ic scattering from 
Gd2Mo2 07. Note the broadening of the reflection near q — 1 
which is complete by 2.7 GPa (Mirebeau et aL, 2006a) 



new and exciting physics, but experiment alists will have 
to be more diligent in their characterization of their sam- 
ples, since it is now clear that the magnetic ground state 
is very sensitive to small perturbations. Understanding 
the role local disorder and spin-lattice interactions play 
in these materials, as well as the origin of the persistent 
spin dynamics seen well below apparent ordering tem- 
peratures, usually by a local próbe such as /iSR is crucial 
in advancing this field. Another area that researchers 
should invest some time studying is the transport prop- 
erties at metal/insulator or superconducting transitions. 
While a connection to geometrie frustration is often a 
feature of the introduction to papers involving supercon- 
ducting pyrochlores like Cd2Re207 (Hanawa et al.^ 2001; 
Sakai et al.^ 2001) and the serieš MOS2O6, where M = 
K, Rb and Cs (Yonezawa et al.^ 2004, a, b), there is no 
credible evidence yet that magnetism, frustrated or oth- 
erwise, play s any role in the superconducting phase. Fi- 
nally, while the models, created by theorists, for under- 
standing geometrically frustrated systems have become 
increasingly sophisticated, there persist many outstand- 
ing questions. What is the origin of the 15 K anomaly 
in the titanate spin ices? Why does Tb2TÍ207 remain 
dynamic and what is the true náture of the Tb^+ mo- 
ments? What is the true ground state for Yb2TÍ207? Is 
the ground state of Er2TÍ207 truly selected by quantum 
order-by-disorder? These must be addressed in the near 
future. 
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